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Synopsis 


Waveguide fed slot arrays are extensively used in many Radar and microwave 
communication systems. Though there exist many types of slots which could 
possibly be used as radiating elements in an array, the longitudinal broad wall 
slot and the inclined narrow wall slot have received extensive consideration. 
Two other slots, namely the centered longitudinal slot and the vertical narrow 
wall non-wrap-around slot are recognized as ideal candidates to build an array 
due to their properties of having no cross-polarization and the ability to build 
a true linear array. The vertical narrow wall waveguide slot in particular is 
preferred, specifically in all applications where scanning in E-plane is desired. 
There have been a number of attempts at making the narrow wall waveguide 
slot radiate, by attaching metallic wires in the neighborhood of the slot. Specif- 
ically, this technique makes use of either a bent screw or a pair of tilted wires, 
attached to the edges of the slot to achieve the necessary coupling from the 
incident wave. 

This thesis introduces a fundamentally different method of exciting the verti- 
cal narrow wall slot. It makes use of a thin printed strip behind the slot to excite 
it. The strip can either be inclined with respect to the slot or can take up an 
L shape. The slot, along with the exciting strip, can be printed simultaneously 
on a dielectric substrate and can replace the narrow wall of the waveguide. 
This method possesses all the advantages of the printed circuit technique and 
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therefore can achieve good accuracy and square ended slots, which are gener- 
ally hard to achieve by the conventional slot machining process. I lu‘ coupling 
to the slot can be adjusted predominantly by varying the strip inclination with 
respect to the slot. There exist several possible paranu'ters in tlu' structure 
which could be made use of in achieving the desired excitation hwel. 

For the analysis, this structure is modeled as a coupled integral eciuation and 
is solved by the Method of Moments (MoM). The integral eciuations are sed up by 
jjjaking use of the Green's functions of an inhomogoncously filled w<iveguidc. 
The required Green’s functions are derived using the standard tc^chniquG of 
normal mode expansions. The analysis of this structure via. MoM turns out to 
be extremely involved due to the dependence of the waveguide Green’s 

functions. This, together with the inclination of the strip along the axis of the 
waveguide, makes all the source and field domain variables coupled, resulting 
in serious implementation problems in setting up the MoM matrix. This tliesis 
presents a systematic quadrature scheme for making this problem tractable 
and computationally efficient. The thesis also introduces algorithms for solving 
the eigenvalues of LSE and LSM normal modes in an inhomogcneously filled 
waveguide. 

Among the Green’s functions made use of in the analysis, it is found that, 
the Green’s function corresponding to the axial E field duo to an axial olectric 
current is a true distribution, due to the presence of the Dirac delta function. 
Method of computing the field, in cases where one of tho source domain 
quadratures are dropped, leads to divergence in the series representing the E, 
field. This problem is identified as one involving the correct interpretation of 
distributions, and is solved by representing the singularity as an or! lumormat- 
ized LSM series. Since the LSM eigenfunctions are not orthonormal to begin 
with, they are orthonormalized using Gram-Schmidt procedure. Algorithms 
for orthonormalization of functions through Gram-Schmidt procedure has been 
developed. This result will have implications in cases where an axial radiating 
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source in an inhomogeneously filled waveguide is considered. 

A program is developed which computes the S parameters, equivalent circuit, 
slot voltage and higher order mode scattering coefficients etc. Extensive para- 
metric studies are conducted by varying strip offset, strip angle, strip length, 
slot width, substrate thickness and substrate permittivity. 

In order to validate the theoretical model, experiments are carried out at 
X-band. Waveguides with a part of the narrow wall machined off are used to 
fabricate the test pieces. A dielectric substrate, printed with the slot and the 
exciting strip on either side, is fitted into the waveguide in the region where the 
narrow wall is machined off Through Reflect Line (TRL) calibration is used to 
establish the measurement plane at the center of the slot. The necessary TRL 
calibration standards are fabricated in the laboratory. The measurements are 
done on a HP8510C vector network analyzer. Experiments are also conducted 
to ascertain the effect of the ground plane. 

The present study on the strip excited waveguide slot indicates the suitability 
of this structure in building planar arrays working at Ku-band and beyond. It 
is possible to build a planar array, by attaching, a printed dielectric substrate 
consisting of slots and exciting strips to a machined waveguide channel. An 
array of this type has all the slots along the array axis and does not suffer from 
the problem of cross-polarization. 
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Introduction 


1.1 Introduction 

Rectangular waveguide slot array antenna has been quite popular in large high 
power radars. Although a variety of slot configurations on a rectangular waveg- 
uide have been studied and reported in the literature, the broad wall offset slot 
and the narrow wall inclined slot have been the configurations chosen in most : 
of the practical antennas, for their simphcity. For designing a high perfor- 
mance antenna with very low side-lobes, the accuracy of excitation control is of 
paramount importance. Thus, a very accurate evaluation of the characteristics ' 
of the slot has received considerable attention among the researchers. 

There seems to be always a conflict between the slot geometiy that can be ; 
analyzed relatively easily, and the slot geometry that can be fabricated easily. 
For example, a rovmd ended slot is easier to machine than a square ended slot; i 
but the analysis of a roxmd ended slot is very difficult compared to the square : 
ended one. Ultimately, firom the cost considerations, the geometiy that is easier | 
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to fabricate is invariably chosen in practical antennas. 


Assuming that accurate slot characterization is available, the rectangular 
wavegmde slot array has some hmitation, while using broad wall offset slots 
or inclined narrow wall slots in array configuration. In the first case, the array 
is no longer a linear array because of positive and negative offsets as well as 
the need to vary the offsets along the array. This introduces a second order 
perturbation in the pattern which can be troublesome while designing low side 
lobe antenna. In the second case, the slot inclinations produce a cross-polar 
component off axis, which cannot be reduced beyond a certain level. 

The rectangular waveguide slot reported in this thesis is the result of a search 
for an ideal radiating element that can ameliorate most of the drawbacks men- 
tioned earlier. However, this geometry is best suited for antenna arrays at 
higher microwave and millimeter wave frequencies, where the smaller size of 
the anteima allows the application of printed circuit technology in the fabrica- 
tion. 


The present work introduces a narrow wall vertical slot which is excited 
by a short strip, printed and kept behind the slot. The strip can be inclined 
with respect to the slot, or takes an L shape. The strip and the slot can 
he printed simultaneously on a dielectric substrate and can replace the narrow 
waU. Since the strip does not touch the waveguide walls, this excitation scheme 
can be regarded as electromagneticaliy coupled. The resonant frequency and 
the admittance level of the slot can be controlled, over a wide range, by changing 
the parameters of the strip. The use of printed circuit technique gives good 
dimensional accuracy, square ends and no wall thickness— features difficult to 
achieve with the conventional slot machining process. 
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1.2 Background and Motivation 

1.2.1 Conventional Radiating Slots 



D - Centered Rotated 
E - Offset Rotated ( Compound ) 

F - Narrow wall vertical ( Non Radiating ) 

G - Narrow wall inclined 
H - Vertical Non wraparound ( Non Radiating ) 

Figure 1.1: Types of waveguide slots. 

A slot on a waveguide radiates if it is cut across the current flow lines of the 
dominant mode. Of the eight slots shown in Fig. 1.1, all the slots with the 
exception of (A), (F) and <H) are radiating. For example, the slot (A) does not 
radiate, because the TExo wall currents are zero at the broad wall center line. 
The slots (F) and (H), do not radiate since they are cut along the ctirrent flow 
lines of the TJSio mode. 

This principle, namely, a slot radiates if it is cut across the current flow lines 
of the dominant mode, can be used to identify the slots which radiate and the 
ones which do not. The same rule can be used, to decide how to perturb the 
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fields around the slot, so that a non-radiating slot is made radiating. 

The wall current in the neighborhood of the slot decides the magnitude and 
phase of the radiated field. Consequently, the slot voltage can be adjusted in 
magnitude and phase, by changing the offset and/or inclination of the slot. For 
example, in the longitudinal offset slot (B), the radiated power goes up with the 
offset from the center line. For the compound slot (E), both the rotation and the 
location of the slot decides the slot voltage. 

A waveguide with a series of slots is a good approximation to a linear array. 
A planar array is formed by stacking several linear arrays and feeding them 
appropriately. Since the voltage across each slot radiator is adjustable, an array 
of slots can synthesize a variety of radiation patterns. 

The slots shown in Fig. 1.1, with the exception of the non-radiating ones, be- 
long to the conventional slot repertoire. These slots, especially the longitudinal 
offset slot (B) has received wide attention and is well xinderstood. The research 
on slot antennas began with Stevenson [1], and can be traced along the works 
of Ohner [2], Vu Khac [3], Yee [4], Lyon and Sangster [5], Park [6], Elliott [7], 
Joseffson [8], Rengarajan [9] and Joubert [10]. 

The literature on conventional slots is not surveyed here, since the present 
work falls in the class on non-radiating slots, and the attempts at making 
them radiating. A brief survey of longitudinal slots can be found in [7, 11]. 
Certain sigmficant developments which are not mentioned in these references 
are discussed below. 

The compound slot (E) can achieve any complex slot voltage, since the location 
and the angle can be varied. This property makes it attractive as a coupling 
and radiating element in a planar array. The principles of compound slots are 
now sufficiently well developed [12, 13, 14], so that compound slots can be used 
for realizing a variety of patterns. 

The effect of TF 20 mode, on the performance of arrays have been studied by 
Rengarajan [15] and Elliott [16] . The array design methodology of Elliott [16, 
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17, 18] is capable of incorporating the internal and external mutual coupling. 



Figure 1.2: Slots on inhomogeneously filled waveguides [10]. 


For a weakly excited slot, two standard assumptions on the longitudinal 
offset slot analysis were shown to be erroneous by Espino [19] and Me Na- 
mara [20]. It was shown that, the longitudinal component of the slot E-field 
cannot be ignored, and the transverse variation of the transverse E-field is not 
necessarily even about the slot center line. 

Joubert [11] has analyzed slots on inhomogeneously filled wavegmdes. [See 
Fig. 1.2]. This work is relevant to the thesis, because the analysis of the strip 
excited slot also involves an inhomogeneously filled waveguide. 

1.2.2 Broad Wall Slots tliat do not Radiate 

In an array of centered slots, all the slots, line up along the broad wall center 
line. Therefore, the array is strictly linear unlike an array of longitudinal offset 
slots. Since all the slots are identical, the external mutual coupling can be 
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Figure 1.3: Excitation schemes for centered longitudinal slots [21, 22, 23], 

incorporated once for all in the array design. These properties have prompted 
the following attempts for exciting the centered longitudinal slot. 

Fig. 1.3 shows two methods of exciting the longitudinal centered slot. In 
method (A), which is due to Clapp [21], a vertical screw or a post introduces the 
necessary asymmetry in the TExq mode fields. The coupling can be increased 
by increasing the depth of the screw. The polarity of the slot E-fileld can be 
flipped, if the screw is changed to the other side of the slot. 

A. Datta et al. [23] have analyzed a variant of the Clapp method, where they 
consider two fiat strips on either side of the slot. The analysis incorporates the 
waveguide wall thickness and can take care of arbitrary strip locations. 

The method (B) in Fig. 1.3 has been proposed by Tang [22] . In this technique, 
an iris is introduced into the waveguide, such that, the waveguide fields are 
made asymmetrical. The slot radiates if the iris is offset from the broad wall 
center line. 

This method combines various interesting features. The slot is a series 
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R— L—C circuit in shunt, while the iris is a parallel L—C circuit in shunt. Conse- 
quently, the reactances of the slot and the iris tend to cancel each other around 
the resonance. Therefore, the iris excited slot is noticeably broad band. The 
susceptance is relatively insensitive to the the iris offset, so that coupling and 
resonance are somewhat independent. The range of slot voltage control is sig- 
nificant that the author is able to achieve patterns with side lobes levels from 
10-34 dB. 

1.2.3 Narrow Wall Slots that do not Radiate 

The attempts on the centered longitudinal slots remain limited to the iris and 
probe excited geometries or the variants thereof However, the narrow wall 
vertical slots have received better attention. This is partly because narrow 
wall slots are preferred, if the beam is to be scanned in the E-plane. 

Conventional narrow wall slot arrays consist of the inclined slots (G) shown 
in Fig. 1.1. The inclination of the slot leads to a cross-polar component in the 
radiated field. The presence of the cross polarization restricts the maximum 
achievable side lobe level, in a narrow wall inclined slot array. The cross-polar 
effects are more serious, if the slots are made more inclined as in an efficient ! 
array. 

I 

The narrow wall slots are generally cut through the narrow wall, so that a 
part of the slot lies on the broad wall, wrapped arotmd. This happens because | 

I 

the resonant length of the slot usually exceeds the narrow wall dimension of the | 
guide. The wrap arotind is undesirable for making a planar array since the ad- ; 
jacent waveguides cannot be stacked together. It also produces a discontinuous > 
ground plane so that a standard mutual coupling evaluation cannot be carried j 
out. The vertical slots considered up to the year 1989 were of the wrap-arovmd | 

i 

type. 

The vertical non-wrap around slot emerges as the ideal choice for a narrow , 



8 


Introduction 



Figure 1.4: Excitation schemes for narrow wall slots from Clapp (21]. 

wall slot array. The attempts on vertical slot excitation also began with Clapp 
[21]. Apart from certain variations, the ideas on narrow wall slot excitations 
can be traced back to the following ideas of Clapp. 

Fig. 1.4 shows t;he two well known methods of Clapp. In (A) a bent screw is 
used to couple from the TEio mode. In (B) a loop links the two sides of the slot. 
We shall call these as probe fed and loop fed respectively. In the probe method, 
the location and the length of the vertical part of the screw decides the coupling 
to the slot. For the loop fed slot, the effective loop area in the transverse plane 
decides the coupling. Therefore, rotating the screw, or twisting the loop changes 
the radiated power. The polarity of the induced slot voltage can be flipped in 
both cases. 

J. S. Ajioka [24] introduced a method which came to be known as the Ajioka 
element. With reference to Fig. 1.5, the Ajioka element consists of two inclined 
metallic rods attached to the sides of the slot. The other end of the rods are at- 
tached to the upper and lower broad walls respectively. The element is referred 
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Figure 1.5: The Ajioka element [24). 


to as "Polarization pure, rod excited transverse slot element". The inclination 
of the rods make them couple from the TEio mode so that the slot is radiating. 

Ajioka describes the design and implementation of a —40 dB side lobe, 
15% bandwidth, 41 dB gain planar array at S-band. The coupling level was 
adjusted in this array, by altering the position of the rods on the narrow wall. 
The Ajioka element is the first reported use of a vertical narrow wall slot for a 
practical application. 

Hashemi-Yeganeh and Elliott [25] modeled the Ajioka element with thin 
strips instead of rods. Method of moments with pulse and point matching were 
used for the analysis. They showed that, the slot in the Ajioka element need 
not be wrap-around, and can be co n fined to the narrow wall. That is, a slot 
which is shorter than the resonant length, can be made resonant via. loading by 
the rods. On one hand, this observation leads one to consider non-wrap arotmd 
slots made resonant via. suitable loading. But the non- wrap around slots have 
lower bandwidth than full length slots. Thus the advantages of being non-wrap 
.aroTxnd comes at a loss of bandwidth. 
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Figure 1.6; Variants to Ajioka and Clapp methods [26]. 


Two variants of the methods envisaged by Clapp and Ajioka are considered 
in [26]. As shown in Fig. 1.6, the slot is excited by a dipole or a loop kept within 
the waveguide. A pair of horizontal wires attach the dipole/loop to the slot. In 
the method described in this reference, the structure is viewed as a dipole/loop 
within the waveguide, loaded by a slot via. a two wire transmission line. Since 
the dipole cannot interact with the slot except via. the transmission line, the 
problem is decoupled and solved. That is, the impedance of the slot is found, 
which is then transferred to the dipole. The problem is thus reduced to that of 
the scattering by a loaded dipole within a waveguide. 

There have been attempts at modifying the slot geometry itself, to counter 
the problems of wrap-around. The C-slot [28] and the I-slot [27] fall into this 
category. As shown in Fig. 1.7, the C-slot and the I-slot fit entirely into the 
narrow wall. It is possible to meet the resonant length due to the essential 
shape of the slot. Both C-slot and the I-slot has cross polar radiation and C-slot 
in particular, suffers from a poorly defined phase center. 
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Figure 1.7: The I and C type narrow wall slots [27, 28]. 


A modification of Ajioka element has appeared recently as a patent [29]. 
They consider vertical wrap-around slots, excited by flat vertical plates. Each 
plate is right triangular such that the base and the height are kept fused to 
the waveguide walls. The technique appears close to the Ajioka element in the 
sense that the inclined wires are replaced by triangular plates. 


1.3 Strip Excited Waveguide Slot 

In this thesis, a fundamentally different method for exciting the narrow wall 
vertical non-wrap aroxmd slot is considered. Fig. 1.8 to 1.9 shows two views 
of the strip excited slot. The slot is excited by a metallic strip printed and 
kept behind the slot. The element can be viewed as an inhomogeneously filled 
waveguide with an inclined strip at the air-dielectric interface. 

It is possible to fabricate this element by etching the strip and the slot. 
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Figure 1.8: Printed strip excited narrow wall slot. 



Figure 1.9: Cutaway view of printed strip excited slot. 
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simultaneously on a dielectric substrate. The use of printed circioit technique 
gives good dimensional accuracy and square ended slots. These are difficult 
and expensive to achieve with the conventional slot machining process. 

The dominant mode of a guide with an E-plane dielectric slab is the LSEiq 
mode. The presence of a thin, low substrate does not significantly alter the 
mode structure of the waveguide. Therefore, the LSEiq mode is very close to 
the T Eiq mode of the homogeneously filled waveguide. 

The behavior of the strip excited slot can be explained making use of the 
following intuitive ideas. The LSEiq mode electric field is vertical and it drives 
a current along the strip. The strip current scatters on its own and alters the 
current flow pattern in the neighborhood of the slot. Note that, for exciting 
the vertical slot, it is sufficient to generate a vertical magnetic field. Thus, the 
incident wave couples to the vertical slot via. the inclined strip. 



A - Inclined centered 
B - Inclined offset vertical 
C - Inclined offset horizontal 
D - Inclined compound offset 
E- L shaped 


Figure 1.10: Different types of exciting strips. 


The strip excited slot has numerous parameters that can be used to change 
the coupling level and the resonant frequency. The behavior of this element can 
be altered, by changing any of the following parameters; strip dimensions, slot 
dimensions, position of the strip with respect to the slot, substrate thickness or 
substrate permittivity. 

If the slot dimensions and the substrate parameters are kept unaltered, a 




14 


introduction 


substantial range of coupling can be achieved, by changing the strip inclination 
alone. The resonant frequency and the coupling are not independently control- 
lable. However, a specific admittance at a given frequency can be achieved, by 
adjusting the strip length and the strip inclination together. It is also possible 
to flip the slot voltage by altering the strip inclination from positive to negative. 

The slot can also excited by an L shaped strip instead of an inclined strip. 
Fig. 1.10 shows an L shaped strip, exciting the slot. The L-strip couples from 
the incident wave through the vertical part, while the strip to slot coupling 
takes place along the horizontal part. The inclined strip can also be moved 
relative to the slot, vertically, horizontally or both ways, as shown in Fig. 1.10. 
The slot radiates a little even in the extreme when the strip does not overlap 
the slot anywhere. 

Fig. 1.11 shows the details of a suggested planar array schematic. The array 
is composed of a machined waveguide channel and a printed dielectric sub- 
strate. The upper side of the substrate is etched with radiating slots, whereas 
the bottom side consists of the exciting strips. The strips and the slots can be 
etched together as in a standard printed circuit fabrication. The substrate can 
be attached to the waveguide bottom channel, with screws via. the set of plated 
through holes on the substrate. This configuration is ideally suited for higher 
microwave frequencies including the millimeter waves. 


1.4 Organization of the Thesis 

This thesis investigates the behavior of the strip excited slot using both the- 
oretical and experimental techniques. The strip excited slot has numerous 
parameters nine parameters: strip and slot dimensions, strip angle, strip off- 
sets, substrate thickness and permittivity — that affect its behavior. Therefore, 
theoretical results are made use of for an extensive parameter wise characteri- 
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slots ( top Side ) exciting strips ( bottom side ) 



Plated through holes 



Figure 1.11: A suggested planar auray configuration. 


zation. The experimental results are used, primarily to validate the computed 
results. However, experimental results are used wherever computed results 
are unavailable. This includes the effect of the groimd plane, losses etc. 

In chapter 2, the strip excited slot is modeled as a boundary value problem 
in an inhomogeneously filled waveguide. A pair of coupled integral equations 
are used to enforce the boundary conditions, with Green’s functions as kernels. 
Green’s functions, for an inhomogeneously filled waveguide, and the half space 
are derived and elaborated upon in Chapter 3. 

The implementation of this problem via. MoM is considerably involved, since 
the waveguide Green’s functions contain the ' factor. For an inclined strip 
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in the axial direction, all the variables participating in the MoM matrix quadra- 
tures gets tied up. Since a total of 9 distinct Green’s function are necessary 
for the solution, the closed form evaluation of MoM matrix entries becomes 
formidable. 

An algorithm has been developed for implementing the MoM matrix quadra- 
tures. It is based on the identification of a multiple integral as a binary tree, 
and the exponential representation of trigonometric functions. This chapter 
also contains the computation of secondary parameters before dealing with 
coding and runtime details. 

Chapter 3 deals with two kinds of topics: supplementary results necessary 
for the analysis of Chapter 2 and topics of an independent nature. Green’s 
functions in an inhomogeneously filled waveguide are derived using the method 
of normal mode expansions. Since, the Green’s function is a series in terms of 
the LSE and LSM modes, the solution of the LSE and LSM eigenvalues are 
dealt with subsequently. The eigenvalues are found using a combined Newton- 
Raphson-Bisection algorithm with good starting values and tight bounds. 

Section 3.4 deals with a topic which has not been done elsewhere. The Eg{Jz) 
Green’s function is identified as a distribution since it contains a Dirac delta 
term, in addition to the LS E and LS M series. For a thin strip, the source domain 
quadrature along the thickness of the strip is neglected. This makes the Eg field 
to remain as a distribution. This problem manifests itself as a divergent LSM 
summation. It is shown that, for the Eg expressions to make sense (classically), 
it is necessary to represent the Dirac delta in terms of an orthonormalized 
LSM series. An algorithm for implementing the Gram-Schmidt procedure is 
developed in this connection. The numerical subtraction of Dirac delta is also 
shown to be an intrinsically unstable operation. 

Method of moment (MoM) is widely used for solving integral equations in 
electromagnetics (EM). A careful evaluation of MoM as numerical method has 
not been carried out so far. MoM involves the straightforward implementation 
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of a very naive idea: the direct solution of the first kind integral equation. Since 
this problem belongs to the class of ill-posed problems, the solution becomes 
unstable after reaching a certain accuracy. In EM, this problem does not show 
up often (or goes unreported) partly because the right-hand- side is very smooth 
and known analytically while the kernel is not so smooth. 

In Section 3.6, we illustrate the instability in MoM based on our examples 
as well as the works of S. W. Lee, R. Mittra and M. Leroy [30, 31, 32] in 
connection with the relative convergence phenomenon. The results indicate 
that the computed results depend on a parameter k , and that the optimum 
choice of « is a compromise between accuracy and stability. The instability in 
MoM is interpreted from the point of Tikhonov regularization. 

Chapter 4 is on fabrication, measurements and the design of appropriate 
calibration standards. The test pieces are fabricated using a standard X-band 
waveguide with a part of the narrow wall machined off The substrate is kept 
pressed from inside the waveguide, and is attached to the guide wall via. sol- 
dering. The S parameters are measured using a vector network analyzer using 
TRL standards, since the conventional waveguide standards are inconvenient 
for setting up the measurement plane in an inhomogeneously filled waveguide. 
This chapter contains the design and implementation of the appropriate TRL 
calibration standards. The effects of the ground plane and losses in the element 
are also illustrated. 

Chapter 5 summarizes the behavior of the strip excited slot based on the 
computed and measured results of the previous chapters. The computed results 
are validated by comparing it with a select set of measurement. It is shown 
that the strip excited slot is a shunt discontinuity which goes through a series 
resonance. It is also shown that, the strip length and strip inclination can be 
used to realize a wide range of operating points. The effect of all the parameters 
in the strip excited slot are presented and interpreted qualitatively. 

Chapter 6 deals with the conclusions and the ideas for future work. 



Chapter 2 


Analysis of Strip Excited 
Waveguide Slot 


2.1 Introduction 


This chapter contains the complete discussion on the formulation, setting up 
and the solution of the strip excited slot as a boundary value problem. The for- 
mulation begins with the geometrical descriptions and the coordinate systems, 
and proceed to model the strip excited slot as a pair of coupled integral equa- 
tions. The integral equations enforce the boundary conditions, and the Green’s 
functions act as the kernels. The analysis requires 9 waveguide Green’s func- 
tions and one external half space Green’s function. These are derived and 
elaborated upon in the next chapter. Chapter 3 contains various issues related 
to the present analysis that it serves as a companion chapter to the present. 
Method of Moments is used to approximate the integral equations to a matrix 
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equation. Throughout the discussion, the structural aspect of the problem is 
emphasized and the necessary S5mibolism is developed. 

The major problem in implementing the MoM hes in being able to do the MoM 
matrix quadratures analytically. The presence of the factor \^ ~ ^'\ and the 
inclination of the strip along the waveguide axis lead to serious difficulties in 
representing the quadratures. An algorithm is developed for implementing 
these quadratures based on the realization that it is a problem in developing 
the requisite symbolism. The algorithm is based on 2 key ideas: the exponential 
representation of trigonometric functions, and the identification of a multiple 
integral with a binary tree. The major exceptions (division by zero, underflow 
etc.) and the restrictions of the algorithm are discussed. The number of func- 
tion evaluations in the quadrature are estimated which indirectly suggest the 
simplifications possible in the implementation. 

A systematic investigation of the secondary parameters that can be computed 
with the present analysis has been carried out. The set of useful parameters 
include — S matrix, T equivalent network. Slot voltage, LSEiq and LSE 2 Q mode 
coefficients etc. A simplified model is put forward for evaluating the external 
mutual coupling to a first order. A brief discussion on the structure and organi- 
zation of the code is included along with the details of CPU time requirements. 


2.2 Geometry and Coordinate Systems 

Fig. 2.1 shows the notations and the coordinate systems used for the analysis. 
The antenna consists of an infinite inhomogeneously filled wavegmde. The 
dielectric substrate is printed with a slot and a strip as shown in the figure. 
The inhomogeneously filled waveguide is a rectangular wavegmde of size ax b 
with an E-plane dielectric slab of thickness t and relative permittivity e^. The 
(r, y) coordinate system is used to specify quantities related to the waveguide. 




a 



Strip geometry 


Figure 2.1: Geometry and coordinate systems 

The strip of size 210; x 2li makes an angle 9 with the waveguide axis. Since it is 
possible to translate the strip with respect to the slot, we specify the coordinates 
of the center of the strip as {yi,zi). The {u,v) coordinate system is introduced 
at the center of the strip in order to facilitate the evaluation of various MoM 
matrix entries related to the strip. The transformation between (u, v) and (y, z) 
coordinates are given by, 

y = yi otu + 0v 
z = zi — Pu + av 


( 2 . 1 ) 

( 2 . 2 ) 
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u = a{y-y{} - p{z- zi) (2.3) 

= Piv-yi) + a{z- zi) (2.4) 

where /3 = cos 9 and a = sinO. 

The slot is vertical and has dim ensions 2ws x 2/^. It is assumed to be centered 
at z = 0. The slot is generally kept centered on the narrow wall of the waveg- 
uide. However, the slot can be slightly offset on the narrow wall in practice. 
The parameter y, specify the position of the center of the slot on the narrow 
wall, ys is usually 6/2 except when the slot is offset on the narrow wall. 


2.3 Formulation of the Boundary Value Problem 

2.3.1 Introduction 

The strip excited slot can be fed by a propagating LSEio mode from one end 
of the waveguide. It can also be excited by a plane wave incident externally 
on the slot. The first case corresponds to the case when the antenna is used 
for transmission, whereas the latter corresponds to the case when the antenna 
is used for reception. In general, the strip excited slot can be excited by any 
valid mode of the inhomogeneously filled waveguide. It can also be excited by 
another slot kept excited in the neighborhood. This is the approximate model 
used for computing the external mutual coupling in an array. One can thus 
identify various instances of exciting the element depending on the context in 
which the element is used. A formulation of the strip excited slot should be 
able to compute all the parameters of interest in any of these situations. 

While formulating the strip excited slot as a boundary value problem, it is 
not necessary to pay attention to any specific method of excitation. Using the 
MoM, it will be shown that a boxmdary value problem can be reduced to the 
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solution of a matrix equation of the form Ax = b. The specific nature of the 
excitation decides the structure of b, and the matrix A remains independent of 
the excitation scheme. Once the A matrix is formulated it is possible to compute 
any parameter of interest for any specific excitation. For the convenience of 
further discussion, we shall assume that the element is excited by a propagating 
LS Eio mode incident firom z = -oo. This is the primary method of excitation and 
the one that is used while making an S-parameter measurement. The alternate 
methods of excitation and the results therein are discussed in Section 2.5. 

Consider an infinite inhomogeneously filled uniform rectangular waveguide 
containing the strip excited slot. The narrow wall containing the slot is assumed 
to be unboimded, so that the slot is part of an infinite ground plane. We 
shall assume the conductor and' dielectric losses to be zero and no thickness 
whatsoever for the strip and the slot. 

The incident wave causes a current flow on the strip and an E-field distri- 
bution in the slot. These sources can be regarded as induced sources causing 
a scattered field. (The electric field distribution can also be regarded as an 
induced source, by modeling it as a magnetic current, making use of Equiva- 
lence Theorem [33].) The scattered field together with the incident field has to 
satisfy all the necessary boundary conditions. In a boundary value problem, 
the boimdary conditions are enforced to solve for the induces sources. 

The following definitions will be adhered to throughout the discussion — J 
denote the electric current, whereas M denote the magnetic current; E stands 
for the electric field and H stands for the magnetic field; and My denote the 
induced current in the directions indicated by the subscript; and Hy denote 
the fields in the specified directions. 

If the length to width ratio of the strip (slot) is large, it is generally acceptable 
to make certain approximations on the induced currents. It can be assumed that 
the strip current is entirely axial, whereas the slot E-field is entirely transverse. 
Note that, this amounts to considering only and My and neglecting J„ and 
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z o 

Figure 2.2: Induced currents and the incident field 

components. [See Fig. 2.2]. Therefore, it is necessary to consider two induced 
sources, namely the electric current on the strip, and the magnetic current 
My in the slot. 

The total tangential E-field has to be zero on the strip, whereas the total 
tangential magnetic field has to be continuous across the slot. We introduce 
further approximations and enforce the boundary conditions only on the axial 
field components. That is, the total Ey, field has to be zero on the strip, whereas 
the total Hy field has to be continuous at the slot. If the symbol E j m (for example) 
denote the E-field due to the source type M, then the boundary conditions can 
be stated as. 




^y\ju ^y^My ~ 


My 


on the strip 
0 at the slot 


(2.5) 

( 2 . 6 ) 


Equation 2.5 is a condition on Ey at the strip, whereas Equation 2.6 is a 
condition on Hy at the slot. The equations have been arranged such that, 
they suggest the matrix partitioning it will lead to, on application of MoM. 
The excitation is kept on the RHS of Equation 2.5 to 2.6, while the self and 
mutual interactions are on the LHS. Since the incident field does not have an Hy 
component the RHS of Equation 2.6 is zero. Had the slot been excited externally, 
the RHS of Equation 2.6 would be non-zero, while that of Equation 2.5 would 




RSE - stRip SElf interaction 

SOR - Slot On StRip 

ROS - StRip On Slot 

SNT - Slot iNTernal interaction 

SXT - Slot eXTernal interaction 

SSE - Slot SElf interaction 

Figure 2.3: Illustrating the interactions 

be zero. 

Note that Eu\j^ is the component at the strip due to the strip current 
Hy\j^ is the Hy component at the slot due to the strip current Ju and so on. 
Therefore, Ey\j^ is the strip’s self interaction, whereas Hy\j^ is the strip’s effect 
on the slot. Fig. 2.3 illustrates all the interactions between the strip and the slot. 
The abbreviations RSE, ROS etc. explained in Fig, 2.3 will be used to denote 
the submatrices resxilting from MoM formulation. All the field components in 
Equations 2.5-2.6 except the are within the waveguide. is the 

Hy component due to the slot current My via. the external half space. Each 
term in Fig. 2.3 can be associated with the terms in Equations 2.5 to 2.6, so 
that 



RSE 


SOR 


source term 


on the strip 


(2.7) 
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ROS 


+ ^y\My - 


ext 

y 




SNT 


SXT 


SSB 


0 


at the slot 


( 2 . 8 ) 


The scattered field components can be related to their respective sources 
via. Green’s fimctions. Two t 5 T)es of Green’s functions will be used, namely, 
the internal waveguide Green’s functions and the external half space Green’s 
function. Since, Green’s functions relate the E, H fields to the sources in the 
(y, z) coordinate system, it is necessary to express (2.5-2.6) in terms of the (y, z) 
coordinate system. 

With reference to Fig. 2.2, we have 


Jy 

— Ol J>ii 

(2.9) 

J. 

1 

II 

(2.10) 

Eu 

= aEy - pE^ 

(2.11) 


so that 

E.\j, = “ E,\j^ - e S.Ia (2.12) 

Ey\j^ — Ey\jy + Ey\j^ 

+ B.\j^ 

EAm, = “ ■B.Im, - (2-13) 

Hy\j, = ( 2 - 14 ) 


where a ~ sin^ and /? = cos^. By substituting (2.12-2.14) in (2.5-2.6), we get, 
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Section 3.1 contains the derivation and properties of the waveguide Green’s 
functions. These relate the 3 field components Ey, Ex and Hy, to the 3 source 
types namely Jy, and My. Section 3.5 discusses the external half space 
Green’s function. The waveguide Green’s functions Ey(Jy ), Ey{.L),..., Hy{My) 
and the external half space Green’s function Hy^^{My) are defined such that, 


Ey\jy - 

jg, Ey(J,) J, dS' 

(2.17) 

^y\jx = 

J, dS' 

(2.18) 

^y\My = 

1^, Hy{My) My dS' 

(2.19) 

II 

H 

Hf{My) My dS' 

(2.20) 


where stands for the integration in the respective source domain. By sub- 
stituting (2.17-2.20) etc. in (2.15-2.16), we get the coupled integral equations 
in Ju and My. 


^ Ig, Ey{jy)jydS' + j Ey{J,)J,dS' ^ Ex(Jy)JydS' 

^ Is'j, Ey{My)MydS' - ^ Ex{My)MydS' 


-aEl^^ 


(2.21) 


lg,^Hy{Jy)jydS' + HyiJ,)J,dS' + J Hy{My)MydS' 


- Hl^\My)MydS' = 0 ( 2 . 22 ) 

where represent quadrature over the strip, whereas represent quadra- 
ture over the slot. The prime on a domain (say 5^) indicate that it is a soruce 
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domain. 


2.3.2 Method of Moment Solution 

Method of Moments (MoM) is used to solve the coupled integral equations in 
(2.21—2.22). A rigorous discussion of MoM can be found Section 3.6. The 
application of MoM to an integral equation of this type involves two distinct 
approximations: restricting the induced currents to a finite dimensional space, 
and enforcing the boundary conditions on a finite grid. We introduce a basis 
for and My and assume the currents to be approximated by finite linear 
combinations. The choice of basis for strip and slot currents are governed by 
the boundary conditions the currents are to satisfy. The currents are such that, 
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Figure 2.4: Basis and testing functions for the strip 


= 0 at u' = ±^2 


(2.23) 
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Figure 2.5: Basis and testing functions for the slot 


My 

= 0 

at y' 

= 

(2.24) 

U.l 

~ j^l - 


-V2 / , 

as V ±ivt 

(2.25) 

\My\ 

~ - 

- 

/ 

as z ±u‘* . 

(2.26) 


(2.23-2.24) indicate that the currents vanish at the end points of their re- 
spective domains, whereas (2.25-2.26) describe the edge condition. The cur- 
rents shotild have the inverse square root growth at the strip (slot) edges. 
Fig. 2.4 to 2.5 show the basis functions used for the strip and slot currents. 
Entire domain sinusoids are used in the axial direction, whereas Bxi even edge 
current distribution is used along the transverse direction. The subscript "s" is 
used for the basis functions on the slot, while the subscript "1" is used for the 
basis functions on the strip. Therefore, 



Ni 

E 4 fit 4>i 


i=l 


(2.27) 
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2=1 


where 

fti 

(t>i 

k-ii 

fzs 

4^s 


COS k^i u', for 
sin k^l u\ for 

= l/ \/l - {v'/wiY 

= ITT j 2 li 

^ f cosk^siy' - Vs), 
\ sinA:,,(y' - y,), 

= l/ zjl- {z'/WsY 

= iTT l2ls . 


1 = 1,3,... 
i = 2, 4,... 


for z = 1, 3, . . 
for z = 2, 4, . . . 


(2.28) 


A'; and Ns are the number of strip and slot basis functions respectively. Note 
that both even and odd bases (fi for all z) are used along the axial direction. 
This choice of bases accommodates every possible current distribution along 
the axial direction of the strip(slot), but constrain the transverse variation 
to be even symmetric. In general, it is necessary to have flexibility in the 
transverse direction as well. This choice of basis functions is widely used in 
slot discontinuity problems [7, 11]. 


Fig. 2.4 to 2.5 also show the test functions used in reducing the integral 
equations to matrix form. Entire domain sinusoids with constant transverse 
variation are used as test functions. The coupled integral equations are reduced 
to a matrix system Ax = b, by substituting the expansions for and My in 
(2.21-2.22), and t akin g inner product with the test functions. Equations (2.9- 
2.10) are used for expressing Jy and in terms of Equation 2.21 gives 
rise to Ni equations via. inner product with Ni strip test functions, whereas 
Equation 2.22 gives rise to Ns equations via. inner product -with Ns slot test 
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functions. The resulting matrix equation can be arranged as, 

■ RSE SOR 
_ROS SSE 

' ' b ' 

Each of the submatrices derives their name from the interaction diagram of 
Fig. 2.3. In the following expressions, primed quantities denote source level 
variables, whereas unprimed quantities represent field (test) level variables. 
The index i is used on basis functions whereas j is used for testing functions. 
The strip is indicated as S'^ or Sr depending on whether it is a basis/test domain. 
The slot is indicated as S'r or Ss depending on whether it is a basis/test domain. 
A quadratmre like indicates a 4-D quadrature with the source level 

integration on the slot followed by the field level integration on the strip. All 
bold quantities represent intermediate submatrices. The following definitions 
hold good for the matrix partitioning. 

• RSE= stRip SElf interaction matrix of size Ni x iV) 


RSE 

= a RSE 

V - /?RSE, 


RSE, 

= a RSE 

yy - 0RSEy, 


RSEz 

= a RSE 

zy ~ 0 RSE*J 


RSEyy{j, ^) 

= / 
Js%,Sh 

fii 4>i fji 

dS'dS 

RSEy^{j, i) 

= /, 

Byi'^z) fu fjl 

dS'dS 

RSE.y{j, i) 

= / 

Js'^,Sr 

E.{Jy) fji 

dS'dS 

RSE^^{j,i) 

11 

Rz{Jz) fii <I>1 fjl 

dS'dS 


• SOR= Slot On StRip cross interaction matrix of size Ni x jV, 

SOR = a SOR, - /3SOR* 




jginc 



0 


(2.29) 
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SORy{j,i) = / Ey{My)f,s<PsUdSUS 

SOR,{j,i) = [ E,{My)f,s<^sfjidS'dS 


m 


ROS= stRip On Slot cross interaction matrix of size Ng x Ni 


ROS = 

nROSj; - /?ROS, 

ROSy{j,i) = 

f Hy{Jy)fa(i>ifjsdS'dS 


Js'^.Ss 

ROSsiiN) = 

[ Hy{Js)hid^if,sdS'dS 




• SSE= Slot SElf interaction matrix of size Ns x Ns 


SSE = SNT - SXT 

SNT{j,i) = [_ Hy{My)U(f>sfjsdS’dS 

SXT{j,i) = H^\My)U<Psf,sdS'dS 

J s$ 



'/r 



1 

1 

Ii = 

A 

and 

r = 

TS 

^2 


Jki 



=?■ 
1 


• is the incident field vector of size Nixl 

E^%j) = a I El^^fyidS 

If the incident field is LSExq, can be taken as where 7 is the 

propagation constant of the LSEio mode. 

• 0 is the zero vector of size Ns x 1 

In the next subsection, the problem of evaluating the submatrices RSE, 
SOR, ROS and SNT are taken up. The evaluation of SXT is discussed in 
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Section 3.5. The matrices A and b can be computed to solve for the coeffi- 
cient matrix x. Any parameter of interest can be computed once the current 
distributions and My are reconstructed from I* and I®. 


2.4 Quadrature Schemes for MoM Matrix Entries 

2.4.1 Introduction 

Algorithms for implementing the internal MoM element quadratures are de- 
veloped here. For simple MoM problems, it is possible to give closed form 
expressions for the MoM matrix entries. However, the expressions for MoM 
matrix entries for this problem are extremely involved and voluminous. It is 
therefore impractical to quote such expressions, for the reason that, the expres- 
sions can neither be coded efficiently nor does it provide any insight into the 
structure of the integrals. 

Our discussion of quadrature methods will be based on the stRip SElf inter- 
action (RSE) formulation, since this presents the most difficult case. All the 
other matrices, namely SOR, ROS and SNT can be evaluated based on the 
method developed for RSR evaluation. Consider the definition of two typical 
RSE intermediate matrices. 

RSByy{j,i) = f Ey{Jy)fucl>ifjidS'dS (2-30) 

RSEy,{j,i) = f Ey{J,)fu<f>ifjidS'dS (2.31) 

where the Greens functions (See Section 3.1 for derivation and details) are 
given by, 
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^y{Jy) — E ^Ey{jy) sinui(a-x') coskyy' sm.ui{a — x) cos e f I + 
E ^By{jy) smvi{a- x' ) coskyy' smvi{a-x) coskyy ^ 

Ey[Jz) = E smui(a-a:') sin/:^/^' smzii(a-a:) cos /Cj^y e""^' 1^“^' ■ © + 

E C^EyiJ,) smtJi(a-z') sin kyy' sinvi{a-x) cos /c^y e I 0 


These Green’s functions, namely Ey{Jy) and Ey{J 2 ) represent two generic 
cases, so that we can make the following observation regarding the structure of 
Green’s functions in inhomogeneously filled waveguides. Each Green’s function 
consists of 2 double indexed series, one through the LSE modes and the other 
through the LSM modes. Each series has a coefficient part, say and 

a function part. The function part depends on both primed and unprimed 
variables. It is composed of the product of foxir trigonometric functions and an 
exponential factor which depends on \z — z'\. On some Green’s functions, we 
also have the signum function <d{z — z'), which merely changes sign based on 
z — z'. Since the source (field) level integrations are to be carried out over the 
strip (slot), the variables x or x' does not appear in the quadratures. If "tr" 
stands ("tr" denotes the word trigonometric) for sin or cos, then, we are looking 
at integrals of the form, 

I = r, „ tr {ky y) tr {ky y') tr (%z u) tr {ku u') e <i>i 0 dS' dS . (2.32) 

Sr 

In (2.32), we have replaced the basis and testing functions with their def- 
initions in terms of "tr". The evaluation of the entire MoM matrix requires 
9 waveguide Green’s functions which amotmts to the evaluation of 18 double 
series. The integral X or the variants thereof appear within each of these series. 
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Therefore, these integrals gets evaluated once every mode, and 10,000 modes 
or more are generally used to get convergence with respect to modes. This ne- 
cessitates a very efficient method for representing and coding these integrals. 



Figure 2.6: Salient points on the strip 

The integral in (2.32) cannot be evaluated in the present form due to the 
presence of the factors |c - c'| and 0. It is necessary to get rid of the jc -- z'\ 
dependence of the integrand by partitioning the strip into regions of trapezoids, 
where either z > z' or z < z'. Fig. 2.7 illustrates the partitioning scheme used 
on the strip, whereas Fig. 2.6 indicates certain salient points on the strip. It 
will be proved later in the discussion that the integral for I can be computed 
hy evaluating the integrand at these 12 points on the strip. The symbols za, 
^D, Vcr • • stand for the coordinates of the points .4, D, C and so on. Since the 
primed and unprimed variables are on the strip, one can identify 3 cases of z', 
namely zd > z' > za, za > z' > zc and zc > z' > zg. For each of these z' ranges, 
we can subdivide the test domain into a set of 2 or 3 trapezoids. (Triangles 
and rectangles can also be regarded as trapezoids in general). Therefore, the 
quadrature for X can be split over these 8 domains as follows. 



Sr 



case a 



partioning m source domain partioning in lest domain 
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Figure 2.7; Partitioning the strip into trapezoids 



case c 


(2.33) 


The meaning of an integral like /^/ is as follows. Asstiming zd > z* > za, 
evaluate the field point integral This integral depends on both y' and 
z' because of two reasons: the domain of integration depends on z' , and the 
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integrand is a function oft/' and 2 '. Evaluate this over the source domain S[ to 
get In (2.33), the + sign is chosen if 0 is absent and - sign otherwise. 

The evaluation of any integral in (2.33) is similar. Therefore, we choose a 
typical term and continue the discussion. I is redefined as tlie part of the 
integral over Si and 5^. Therefore, 

I = / tr {ky y) tr {ky y') tr {k^i u) tr {ku u') (pi f/.V' dS (2.34) 

Si 

It is possible (in theory) to carry out the integration in (2.34) analytically. 
However it will turn out to be an arduous task for the following reason. A 4-D 
quadrature of this kind is evaluated as a sequence of 4 iterated integrals. Each 
intermediate quadrature will be referred to as a partial quadrature. Therefore, 
there exist an order for the quadrature say u, v, u' and v'. Since the function <pi 
depends only on v', and the partial quadrature with respect to 1 / is done at the 
end, we can ignore (pi until the last integral. For doing a partial quadrature, 
say in u, we have to use trigonometric product to sum formulae and express the 
integral in the form, 

r%““tr(6u)du. (2.35) 

J Ui 

Note that, the application of trigonometric product to sum formulae for n terms 
lead to 2”“^ terms in the integrand. The indefinite integral of (2.35) doubles the 
number of terms and the limit substitution again doubles the terms. Therefore, 
at the completion of just one partial quadrature we are left with an elaborate 
number of integrals. The situation is made worse because the domain of in- 
tegration depends on the variables awaiting for the quadrature. One has to 
keep track of this interaction between the terms of the integrand to get the 
parameters a and b. If we continue this thought experiment, it is obvious that 
the expressions becomes unmanageable. In fact, it can be shown that the in- 
tegrand becomes the sum of 8 x 4® = 512 terms prior to the v' quadrature. 
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This a fundamental problem in evaluating the waveguide Green’s functions on 
inclined domains [34]. 

There is yet another problem which is intrinsic to the integrand in (2.34). 
Since the integrand is a product of 4 functions of the "tr" type, there exist 
2^ = 16 possibilities for the integrand. Even if the integral for one specific case 
is evaluated, there remains the task of evaluating the rest 15 cases. 

Note that, apart from the inability to keep track of this large number of 
terms, this problem does not present any profound mathematical puzzles. An 
attempt at evaluating this quadrature is discussed in [35]. The method makes 
use of an integral representation of the waveguide Green’s function such that 
the factor is absent. This technique however is applicable only to a 

select set of integrands and cannot be applied if the edge current distribution 
is used. The method also necessitates the investigation of the complex plane 
behavior of the integral (poles, branch cut etc. ). The algorithm that is proposed 
here is based on the conviction that, this is a problem in clever book keeping and 
one that eventually needs a coding solution, than the mathematical niceties the 
author suggests. 


2.4.2 Algorithm and Illustration 

The algorithm for the evaluation of integrals of the type (2.34) are presented 
here. Let bold j denote since j is already used for indexing the test 

functions. Consider the definition of X with a few modifications, 

X = [ triky) tr (kju) tr{ky') tT{kiu') e~'^ ~ 4> dS' dS (2.36) 

Js[,Si 


n% 

T’ 


where 


k 


for n= 1, 2, 3, . . . 
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k, = for j= 1.2,3,... 

4, = <t>{v') = 

and 7 is assumed to be real and non-negative. The steps of tlie algorithm 
evolves naturally from the following step: replace the four term trigonometric 
product in the integrand of (2.36) by the exponential representation. That is, 
we regard the term tr [ky) tr (fcju) tr {ky') tr [k,u') as a linear combination of 
gj {ky±kju±ky' ±kxu') ^ If WO define, 


5 = / , ^j{ky ± kjU ± ky' ± k u') ^ - 7 {z ~ z') ( 2 . 37 ) 

JS'i, Si 

then there exist 8 possibilities in the integral for £. Let , . . , denote 

the corresponding integrals such that, 









gj{ky — kju — ky' - kiu') ^ "'i' ,(5' f/5 
gj{ky - kjU - ky' k,u') ^ ^s' dS 

gj{ky + kjU + ky' + ka') 


It then follows that. 


^ = ±167 ± ± ± ± ± £‘ ± £• ± 

± £1 ± £2 ± £3 ± £4 i £5 ± £6 ± £7 ) J (2.38) 


for some choice of signs. The term j is absent in some cases. In this represen- 
tation stands for the complex conjugate of £”. The first implication of this 
result is as follows. Once we evaluate the 8 integrals £®, £^ . . . , £'’' it is possible 
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to evaluate X for any of the 16 choices in the integrand. For example, if the 
integrand is composed of "cos {k y) cos [k^ u) sin {k y') sin {k^ u')", then we have 
the following choice of signs. 


X 


(^0 _ _ ^2 ^ ^3 ^4 _ ^5 _ ^6 ^ ^7 ^ 

(^0 - gl - gl + ^3 ^ ^4 _ 55 _ 56 + 57)1 


Therefore, the type of tr’s involved in X can be ignored, and one can concentrate 
on an integral of the generic form. 


f = Jg, g X ‘I’ ‘‘S'dS 

A'’ = qJ i^y ~k 'k S^kiU') Q—j{z — z') 


(2.39) 

(2.40) 


where si, 52 and 53 are ± 1. 

Two observations can be made regarding the integral in (2.39): one regarding 
the structure of the integrand, and the other related to the domain of integra- 
tion. 

• The integrand is a complex function of 4 real variables: two primed variables 
and 2 unprimed variables. The integrand is the product of an exponential 
function X and the edge current function (j). The function X can be regarded as 
one exponential function of a composite argument. The argument is a complex 
linear combination of the integration variables. 

• The domain of integration is a sequence of 2 trapezoids. Each trapezoid is 
defined by a set of 4 straight lines. These bounding lines are linear functions 
of the variables awaiting for the quadrature. 

For example, the domain of integration for the quadrature is defined 
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Figure 2.8: Descriptions of the domains S'^ and Si 


as, [See Fig. 2.8] 






~li < u < "i" sP 

Vp < V < 

-h < u' < Ct;' + 
-Wi < v' < i-wi 


where 


c 


Vp 


= a/(3 

j 

~ -L-ZJi 

_ 

P 

Hz' ~ Zi) - 0 k 
a 


The quadrature is done after choosing a quadrature sequence. For RSE 
matrix evaluation, it is necessary to do the partial quadratures in the sequence 
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X 


I — 








-X, 




u’ 


Xj 






u’ |i" 


X, 


's I? 






[o 




V ’ 


V ^ 


V ^ 




V ^ 


V’ 


^3 


^2 
V ' 


1 


V' 




Figure 2.9: Quadrature sequence as a binary tree 


71 W; // X 7 
TTW; /o" ^ 

nw, 

Itw, I* X4 

nwi ll X, 

nwi X^ 
Jltv, /„' X, 
it«,, /; 


u, V, u' and v' . Since the quadrature in v' is done at the end, the term (p can 
be ignored from the integrand until the integrations with respect to u, v and 1 
u' are completed. Consider the quadrature with respect to u. If we regard the 
integrand X as e^'“ ^ then the integral consists of 2 terms of similar structure ■ 
evaluated at the upper and lower limits of the integration variable. That is, j 

t 

I 

X Xo] (2.41) i 

ACh 1 

1 

where X^ stands for X evaluated at the upper hmit of the integration variable | 

and Xq stands for X evaluated at the lower limit of the integration variable. 
Note that, the integrand retains this structure after every partial quadrature | 
because of the nature of the integrand as well as the domain of integration. | 
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Thus every integration splits the integrand into the abovt‘ foi-ni. A sequence of 
partial quadratures thus forms a binary tree. Fig. 2.9 shows a tk'tailed view of 
this evolution. 

One begins with the integrand Af, perform the quadratun) in ii leading to 
terms Xi and Xq. We associate the term arising out of the upper limit substi- 
tution as the upper branch of the tree (designated 1 ). The term arising out of 
the lower limit substitution is traced along the lower branch (designated 0). 
Therefore at the end of 3 quadratures there are just 8 term.s A’o, .I'l, . . . , AV 

The state of the integrand is specified by tracing down the branch in binary. 
For example, the term Xq corresponds to the term A' evaluated at the upper 
limit of u (1 ), upper limit of t; (1) and the lower limit of u' (0). Similar numbering 
scheme is used for the coefficients of the quadrature variables. 

The quadrature in v' is different from its predecessors. At the completion 
the 3 previous quadratures, we are left with an integral of the following form, 



e (j) {v') dv' 


(2.42) 


Using the transformation v' = wi cos (p, this integral can be reduced to this 
known form, 

'1 .TT 

Trio 


TTWi 


C08<p 


(2,43) 


kikv' wt) 


where Io(kv' wi) is the modified Bessel function of the complex argument A',,' wi. 
Thus the effect of the quadrature in v' is to replace the integrand A' by k tvi Iq A'. 
X stands for X evaluated at v' = 0. Note that this operation does not increase 
the number of terms. This procedure is sufficiently systematic that it is possible 
to write down the sequence (by inspection) as follows; 




1 

K 


U 
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1 

TT 


1 

kyl 


^(A', - X,) - - X,) 

ito it, ~ 


IT W I 
JZy 


1 

kyO 


tn I ^ ~ ^ 

^ (-^0 ^^3 - .^0 -^ 2 ) — (-^0 ~ ^0 ^ 0 ) 


If all the quadrature limits are distinct, then the case illustrated in Fig. 2.9 
is valid. That is, there are eight coefficients for v', four coefficients for u', two 
coefficients for v and one coefficient for u quadratures. However, if the domains 
are not distinct like this case (Note that the RSE domain, (See Fig. 2.7— 2.8) 
is bounded by lines having only 3 distinct slopes), not all the coefficients are 
distinct. It can be shown that there exist the following (distinct) coefficients in 
the entire RSE quadrature. 


• Distinct coefficients for ku 

1 . ±72 + j (ki ± Sikj) 

• Distinct coefficients for 

1. j {k + sikja)/p 

2. ±71+ i k2 

• Distinct coefficients for ku' 

1 . i[(l + 52)^1 + {sikj ± S 3 ^i)] 

2. j[-~k{P^/a) + S2ki + szK] 

3. T 72 + i (s2 ki ± S3 K) 


• Distinct coefficients for kyi 
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1 . + s^kz] 

2. J ( 1 + S2 ) ^2 

3. ±7i + j sz kz 

4. j{szk + szha)ll3 

where 

ki = koi 
kz = kp 

71 = ja 

72 = 7/5 

In the above expressions, the choice of the sign ± is to be made based on whether 
z > z' or z < z' . 

The terms Aq to can be evaluated using a simple idea. The definition 
of Xk is repeated once again. A/., is the value of X evaluated at the kth. stage 
of the binary tree, with v' = 0. For example, the value of Af? is evaluated as 
follows. We identify the set of points on the strip domain fSee Figs 2.6-2.8] 
which belongs to the intersection of the 4 lines. 

u = Cv + , upper limit of u 

v= +wi , upper limit of 

u' = -li , lower limit of u' 

= 0 

With reference to the Fig. 2.6, this takes the unprimed variables to the point 
F and the primed variables to to the point E. We list the value of /To to Xj for 


AV = 1 HHGG 
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1 FFBE 
1 AAGG 
1 BEEF 
eoG DDGG 
Qde DDEE 
1 AAGG 
1 EEEE 

where 


^DG 

= e~ 

-i\zd 

— 

Zg\ 





^DE 

= e~ 

-1 \zd 

— 

Ze\ 





HHGG 

= ei 

(kyn 

+ 

Si kj 

Ujj 

-1- 

^2 k Hq 

53 k^ '^G 

FFEE 

= ei 


+ 

Si kj 

Up 

-t- 


+ S 3 K u'^) 

AAGG 

= eJ 

(kVA 

+ 

sikj 

Ua 

-f- 

S2 k Pq 

+ S3 kx Uq) 


2.4.3 Exception Handling and Restrictions 

In this section, we discuss certain miscellaneous topics related to the quadra- 
ture. These include some important exceptions, restrictions and limitations of 
the algorithm. 

Exception Handling 

(1) Modes other than LSEmJLSMmn (m,n > 1) 

While developing the quadrature scheme, the propagation constant 7 is as- 
sumed to be real and k to be nonzero. However, for LSEmo modes, we have 


^6 = 
^5 = 

^2 = 
= 

= 
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A: = 0 and, in particular for LSEiq mode 7 is imaginary. Tlio algorithm does 
not need any modification for the case of LSEmo mode with ni > 1, except that 
(2.38) is valid in the following simplified form. 

^ (2.44) 

If 7 is imaginary it is necessary to change (2.38) substantially as follows. If 
7 is imaginary and k = 0,we are dealing with an integral of the form, 

^ Is[ Si 

Define, 



^ = ±57 [(^° ± «‘) ± (f'° ± f'M] (2.46) 

(2) Division by Zero 

A division by zero oconxs if any of the coefficients b„, K or b.. becomes zero 
for some combination of parameters. However all the distinct coefficients have 
eady been listed in page 43 that it is possible to identify such cases and 
avoid them . The 6 coefficients can be thought of as belonging toSdistinct groups 

dependingon the possibihty of becoming zero. For example, the coefficient b. = 
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T Oi + i ^2 can never become zero. However the coefficient = j (k + si kj a) / (3 
becomes zero only if /c = k^ a. Such cases can be ignored because it requires a 
strange coincidence of parameters to make it zero. This case can definitely be 
avoided either by computing with extended precision or using parameters with 
an infinite-small offset. For example, we can use 6> = 30 0001° instead of 30°. 
The coefficient 

ku' — i[(l + ^2)^1 + {sikj + 53^:1)] 


however becomes zero for the condition, 

i = j 
S2 = -1 

Si = -S 3 . 


This case definitely occurs on the principal diagonal of the RSE matrix and 
needs to be taken care of The quadrature scheme illustrated in Fig. 2.9 will go 
through the following changes in this case. If ky,' = 0 at the completion of the 
quadrature in v, then it is necessary to evaluate an integral of the form, [See 
Fig. 2.8] 



g kyf 'IL 


du' = u' 


<v'->rU 

i—ll 




(2.47) 


Therefore the v' quadrature becomes. 


T / 

KV + (?., + /,)] dv' = 

“Wl 

7rwif(wi- r COB (p dip + {U + h)- [ ^ dip^ 

V ^ R •^° J 

/l(fc„/ tui) /o(V Wl) 

Therefore, the integral can be reduced to a closed form for this case. The 
modified Bessel functions /„ has the property that, 

17„(2)1 < Io{x) = 0(e+"), for n>0 


i 

I 

i 
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for a: = 1 Re{z) |. Since the arguments to the Bessel functions grow indefinitely 
with respect to the modes, it is necessary to keep the term /„ bounded in the 
quadrature expressions. This can be done by combining the term arising from 
with the Bessel functions. For example, the timm .t ir, /“ .I 2 in 
consists of 


= 7rm,7o(7l +JS2 1-2 
= TTio, /» DEEE 

V S. / 

Note that, it is necessary to treat e~'^^ Jo together .so as to avoid floating 
point overflow at large 7. All standard mathematical libraries return /„(:) as 
In{z) to avoid overflow [36]. 

Restrictions in the Algorithm 

The restrictions in the algorithm arise from our inability to carry out the v' 
quadrature in terms of the Bessel functions, J„, For the cases that has been 
considered so far, it was possible to partition the domain .such that the if' quadra- 
ture could be cast in one of the following forms, 

( 1 ) 

J 

(2) f ' v' <l>{v') dv' 

Fig. 2.10 shows 2 cases where such a partitioning scheme does not hold good. 
If the strip inclinations exceeds then the partitioning scheme of Fig- 2.7 
cannot be apphed for RSE formulation. Therefore, the quadrature in t*' cannot 
be reduced to one of the standard forms, and we have a limit on the maximum 
allowable strip inclination. However, Omaai is of the order of 80® if the length to 
width ratio of the strip is more than 7. 

Similar results hold good if the strip is partially overlapped over the slot. 
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Figure 2.10: Restrictions in the algorithm 

This is shown in Fig. 2.10b. Therefore, cases having arbitrary horizontal offsets 
between strip and the slot cannot be handled. In fact, the offset should be such 
that, 

> +^5 

Zc < -Ws 

Certain Simplifications in the Algorithm 

It was shown in the discussion of the quadrattire that the integral in (2.39) is 
evaluated as a linear combination of Xq to X^ evaluated at the salient points of 
Fig. 2.6. If we observe the partition scheme of Fig. 2.7 and Equation (2.33), it 
can be seen that the contribution from S 2 and S 3 are added together in case(a). 
Therefore, it so happens that the contribution from the point P' cancels while 
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combining the integrals over S 2 and Sz. That is, one need not retain the term 
arising from P' in v quadrature for Therefore, there need to be only 6 

terms in the final integral. Similar observations can be made n-^rarding point 
Q' while combining /s',Si case(c). 

A substantial reduction in computation is achieved by introducing some 
symmetry assumptions. Note however that the result mentioned above does 
not need any symmetry assumptions. If the strip is assumed to be centered on 
the narrow wall, the following results hold good. 

RSE(;,i) = 0, for \j-i\ = odd 

This brings about 50% reduction in the number of matrix elements to be 
computed. There are further simplifications resulting from the symmetry of 
the partitioning scheme. With respect to Fig. 2.7 and (2.33) it can be shown 
that we need to evaluate only the first 4 quadratures. That is, the quadratures 
fs',S2’ fs',S2> Js',S3 can be shown to be the complex conjugates of 

fs',Ss> fs',,S 2 ’ /si.si respectively. 

The cost of a quadrature scheme is specified in terms of the number of 
function evaluations, and not by the multiplication count. For this quadrature 
algorithm, the number of function evaluations are very small. The algorithm 
requires the evaluation of Xq to ^7 and to for 8 sign factors (.si, .^2 and 
S 3 ) over the 8 partitioned domains. However, not all functions are distinct in 
the sense that one can evaluate them from a canonical set, by forming the 
conjugates, taking product etc. An exact estimate of the number of function 
evaluations in (2.32), for the case of? # j, is given in Table 2 . 1 . 

This is the cost of function evaluations for each mode. However, the evalua- 
tion of e7 ^ y and ^ need not be repeated after it has been evaluated for a 

mode because of the curious structure of k and Note that 

, riTT „ 

k = — , for n= 0 , 1 , ... 
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Table 2.1: Function evaluations in the quadrature method 


Function type: e 

lo 

eJky 


No of evaluations: 2 

8 

12 

3 


kj = 0’ for j = 1, 2, . . . 

Therefore, one can evaluate them recursively using relations of the t 3 rpe, 

Qjimr/b)y _ ^j{n-lTr/b)y ^ji-JTjb)y 

If we regard the integral (2.32) as being evaluated alone, the cost per quadra- 
ture is not substantial. However, since the integrals appear within double in- 
dexed series, the cost of MoM matrix evaluation is substantial, mostly because 
of the large multiplication coimt. 

2.5 Computation of Secondary Parameters 

It is possible to compute various parameters of interest once the boundary 
value problem is formulated in the matrix form. The most important param- 
eter is the S matrix with respect to the propagating LSEio mode because it 
can be measxired and can be used for validating the theory. However, other 
parameters such as slot voltage (I^), normalized radiated power {.Prod), LSEio 
and LSE 20 mode coefficients etc. are of interest for array design purpose. Once 
the induced currents are computed, the evaluation of these parameters follows 
automatically from their definitions. Therefore, the formulas for the evaluation 
of the secondary parameters are not quoted, but the relevant implementation 
aspects are explained. 

••JfcA 125679 
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Figure 2.11: Secondary parameter computation schematic 


Fig. 2.11 shows the complete schematic of computing the secondary parame- 
ters. We regard 2: = -00 as port 1 while = +00 is regarded as part 2 . In order 
to compute the S matrix, an LSEw mode is assumed to be incident from port 1 
and the reflected and transmitted dominant mode amplitudes are computed. 
This gives 5 u and S21. S22 and 5 i 2 are computed by assuming an LSEiq mode 
from port 2 . Since each excitation leads to a column b, it is possible to solve the 
S matrix m one step by solving Ax =: b with multiple right-hand side. Stan- 
dard Ax = b solvers allows multiple right-hand sides 136 J. The ISEio mode 
coefficient B^o is identical with whereas the LSE2, mode coefficient B20 can 
be evaluated by extracting the amplitude of the mode. 

The S matrix can be converted to normahzed Z matrix via. the transforma- 
tion {I + s) {I -S) 1. 'Ifm T network representation is preferred if the discon- 
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tinuity is a shunt admittance. The T network is defined by three impedances, 

= - 2^11 — -212 
■22 = -222 ~ 212 
Zt = 212 

where zi, are the series arms and Zt is the shunt arm. The quantity 
1 — 1‘5'iiP — |5'2 i|^ represents the ratio of the radiated power to the incident 
power from port 1 . We shall refer to it as the normalized radiated power Prad- 
If the internal excitation is via. the LSE 20 mode, it is possible to proceed 
exactly as in the case of LSEiq mode incidence and solve for Bio, ^ 20 - The slot 
voltage is obtained by integrating the slot E-field across the width of the slot. 
The slot voltage is evaluated at the center of the slot with assumption that the 
incident wave has a certain magnitude at some reference point. 

The computation of external mutual coupling between two strip excited slots 
is an exceptionally difficult problem. The analysis is considerably simphfied, 
by considering it as an excitation problem from a slot kept excited at the voltage 
Vk. The extraneous slot is regarded as causing an external magnetic field 
The column b will have zero entries corresponding to the rows which enforce 
conditions on the strip. [See Fig. 2.11]. The lower nonzero part of b can be 
computed using the half space Green’s function. 


2.6 Aspects of Coding 

Referring to Fig. 2.12, the code is organized into various translational units and 
kept in different directories^. The Unix utility make is used to compile and 
Hnk various units. Though subroutines have been used wherever necessary, 
they were avoided while writing the quadrature routines fill_matrix_RSE etc. 


^The source code can be made available upon request. It takes about 174 kB memory 
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pbgram main 
set frequency 
read specks 
do MoM — 
stop 
end 


iubroutine do_MoM 

solve LSE LSM roots 

comp Green’s function 

fill matnx RSE 

fill matrix SOR 

fill matrix ROS 

fill matrix SNT 

fill matrix SXT 

fill column b 

assemble A 

solve Ax= b 

compute results 

return 

^d 


subrootino grenn!; function 

do Gram Schmtdl 

get 8 representation 

return 

end 


subroutine comp suits 
comp S-paramoters 
comp 2-paramotors 
make T network equivalent 
comp slot voltage 
comp power radiated 
comp 10, 20 mode coetl. 
return 
end 


Figure 2.12: Code structure 


because of the following reasons. It is found that the quadrature algorithm 
uses all the variables at some stage or other. Therefore, in the strict sense, 
there are no intermediate variables. If a modular structure were opted for the 
quadrature routines, it would require massive storage via. common locations. 
The variables in common are static and carmot be forced to automatic storage. 
Therefore, this algorithm was implemented as one complete code. The compiler 
utility include is used to make the code structured. I’he compiler option 
automatic can be used to force d3naamic storage of all the variables. 


Table 2.2: Specification of the Computing Eiivii-onim'nt 

Machine DEC 3000/600 

Operating System DEC OSP/1 2.0 
Compiler DEC Fortran 77 v 3.3 


Numeric Precision Double 
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Table 2.2 shows the relevant details of the computing environment. A typical 
run at a frequency takes about 9.0 minutes of CPU time for Mmax ~ 100, 
Nmax = 100, Ni — 9, Ns = 9 and = b/2. If the strip is assumed to be centered 
with respect to the slot, then the execution time comes down to 3.3 minutes 
because of the savings due to the symmetry. 


2.7 Summary and Conclusions 

This chapter contained a complete account of the analysis of the strip excited 
slot as a boundary value problem. The discussion began with the development 
of the necessary S5anbolism as notations, coordinate systems etc. The boundary 
value problem was set up using a pair of coupled integral equations in terms of 
the induced currents. The transverse currents were neglected and the boundary 
conditions were enforced only on the axial field components. This simplified the 
model. Suitable acronyms were coined to model the strip and slot interactions 
which were consistently used up to the stage of code organization. 

An algorithm was developed for implementing the MoM matrix quadratures. 
The development of the algorithm was based on the conviction that the method 
should be systematic enough for programming. The operation count of the 
algorithm was made which indicate a surprisingly small number of function 
evaluations. However, the number of multiplications were found to be large 
enough which tend to make the problem computationally intensive. The major 
exceptions in the algorithm were discussed. The method was found to be appli- 
cable only below a critical angle for the strip. It also could not handle arbitrary 
horizontal offsets between the strip and the slot. These two restrictions do not 
seriously limit the operating range of the analysis. 

A brief discussion on the structure of the code and programming practices 
were made. The run time information was dealt with along with a discussion 
of the computing environment. 



Chapter 3 


On Green’s functions and Method 
of Moment 


3.1 Introduction 

This IS an extraordinarily big chapter which contains nearly all the theoretical 
work done in this thesis. The topics dealt here do not have a continuous theme, 
but they are related to the topic of Green’s functions in aonu* way or other, 
ain sections such as (3.2) and (3.5) are essentially an elaboration or an 

xtension of well known works [33, 7], but all other sections aro by tind large 
original contribution by tbc author. 

Many sections of this chapter are of independent interest though this is 
primarily as a companion to Chapter 2. Here is a section-wise summary 
of the contents. In Section 3.2, the Green’s functions in an inhomogeneously 
aveguide are derived using the method of normal mode expansions from 
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Collin [33]. The structure and properties of waveguide Green’s functions decide 
the subject matter of the next two sections. Section 3.3 discusses the solution 
of LSE/LSM eigenvalues based on a combined Newton-Raphson and Bisection 
algorithm. 

Section 3.4 is an illustration of a fundamental problem while handling a 
Green’s function which is a distribution. The Ez{J^) Green’s function is taken 
as an example to illustrate the subtleties in computing E^, if one of the source 
domain quadratures are dropped. The results are interpreted from the point of 
view of distributions and a solution is proposed based on the LSM representa- 
tion of Dirac delta. Section 3.5 contains a method of computing with the half 
space Green’s function. This work is partly based on the work of R. S. Elliott [7] . 

In Section 3.6 we deal with the theory and practice of method of moment. 
MoM is a numerical method which is based on a very naive idea — the direct 
solution of a first kind integral equation. Since this problem is intrinsically 
illconditioned, the sequence of approximate solutions becomes unstable after 
reaching a certain accuracy. In this section, this issue is brought to the limelight 
along with an attempt to consolidate various MoM thumb rules from the point 
of regularization theory. 


3.2 Waveguide Green’s Functions 

3.2.1 Introduction 

The purpose of this section is to derive the waveguide Green’s functions nec- 
essary in setting up the boundary value problem of Chapter 2. The Green’s 
functions in an inhomogeneously filled rectangular waveguide has already been 
derived by [37, 38]. However, it will be re-derived in a form which suits the 
requirement of our analysis. We shall use a method from "Field Theory of 
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Guided Waves" by R. E. Collin [33] which is popular and simple*. The principle 
of derivation is briefly as follows. 

By a normal mode of awaveguide, we understand a valid solution ofMaxwell’s 
equation inside the waveguide having e~"^‘ dependence. The set of normal 
modes are countable and have distinct^ propagation constants If the wave 
guide is homogeneously fllled, the normal modes are degenerate and there are 
two modes with the same propagation constant, namely the T"!*! and TM mode. 
For the inhomogeneously filled waveguide, one can identify two distinct sets 
to which the normal modes belong to: the LSE and LSM modes. The normal 
modes are significant because they form a complete orthonormal set among the 
valid solutions inside the waveguide. There also exist suitable? orthogonality 
relations between the normal modes which enable us to exprims the Green’s 
function as a series of normal modes. Some Green’s functions such as E.{J.) 
do not fit into this form, and contain terms other than the normal modes. The 
outline of this section is to derive the normal modes, and assume the Green’s 
function as a series in terms of them. The coefficients of the normal mode series 
are evaluated by applying the reciprocity theorem. 


3.2.2 Derivation of Normal Modes 

Consider an inhomogeneously filled waveguide as shown in Big. 3. 1 . We desire 
to derive the LSE and LSM normal modes of this waveguide. This is done by 
making use of the vector potentials F* and A* and the relationship between the 
fields and the potentials as enumerated below. 

= 0 

E =: -ivxF 


^Waveguide is assumed to be inliomogeneously filled unless otherwise stated. 


(3.1) 

(3.2) 
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Region.1 (Air) 



Wi/': 

T 






X 


Region.2 (Dielectnc) 

/£r 





Figure 3.1: Iiihomogeneously filled waveguide. 


H = + (VV-F+Zc^f) 

(3.3) 



r 

11 

pR 



and 


(V^-i-k^)A^ = 0 (3.4) 

H = +-VxA (3.5) 

Ai 

E = +-J— (VV-A+k^A) (3.6) 

A = Ax X. 


(3.1— 3.3) will be used to derive the LSE modes, while (3.4— 3.6) will be used to 
derive the LSM modes. The potentials Fx and Ax are defined suitably, such 
that the boundary conditions at the waveguide walls are satisfied, in addition 
to satisfying (3.1) and (3.4) respectively. This is done by the choice. 


Fx 

where f{x) 
+ kQ — ky — uf 


f{x) cos{kyy) 

A sin ui (a — x) , 
B sin U 2 X , 


7 


+ k^ 


k^ 


ul 


a > X > t 
t>x>0 
= 0 


(3.7) i 
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/c, = ^ forn ^ 0. 1.,.. 

^ b 


where 


■Ax 

a{x) 


a(x.) sinikyi/) r^ '’ 
j C cos (.'t (a .r), 

1 D cos ^2 X , 



= 7^ + k'^ " »■ 

ky = ^ forn = 1,2 — 
^ b 


n • r ' t 

( * X 0 

' 0 


(3.8) 


The constants A, B, C and D and ui, « 2 , !’i and will be rtdatod making 
use of further conditions. The potentials and .4^ should satisfy the conti- 
nuity conditions for and Hz at the interface at x ' t. This leads to the 
characterization of LSE and LSM modes as shown below. 

For LSE modes, we have 


A sin uid 

-Aui cos uid 
-ul - kl + kl -I- 7^ 
-ul - fcy + 4- 7^ 

For LSM modes, we have 

C vi sin vi d 
C cos vi d 
-v! - kl + kl + 7 ^ 
-vl ~ kl+k^ Jr 


B . 


~ sm «2 1 

(3.9) 

fr 


B 


— «2 cos «2 f 

(3.10) 

(f 


0 

(3.11) 

0. 

(3.12) 


D 

i!% sin m t 

(3.13) 

«r 

■D cos W 2 f 

(3.14) 

0 

(3.15) 

0. 

(3.16) 


The equations which characterize the modes can be expressed more concisely, 
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by taking the ratio of the first two equations, while subtracting the last two 
equations from each set. This leads to (3.17—3.18) as the characterization of 
LSE roots, and (3.19—3.20) as the characterization of LSM roots. 


tan( ui d ) ^ tan( U 2 1 ) 

Ui U2 

ul-ul = kl (e. - 1) (3.18) 

€rVi taxi{vid) = -V2 tan{v2t) (3.19) 

^1-4 = kl (cr - 1) (3.20) 


Equations 3.17-3.20 admit a countably infinite number of solutions, which 
grow monotonically. The roots can be ordered by indexing with m. We shall 
refer to the roots, later, as the mth LSEfLSM root. The detailed algorithm for 
solving them will be dealt with in the next section. However, some properties 
of the roots are worth emphasizing. If the waveguide is homogeneously filled, 
then. 


mn 

Ui = U2 = Vi =V2 = . 

a 

For an inhomogeneously filled waveguide, all the roots tend to mrcj a asymptot- 
ically. For low filling factors, the LSE and LSM roots are very close to each 
other and are only sHghtly different from m-Kj a. 

The constants B and D can be eliminated from the definition of normal modes 
using (3.9) and (3.14), and put 4 = C = 1, so that, 


where 


{ sinui(a — x) cos(fcj,y ) , a>x>t 

erRte sin 142 3: cos{ ky y ) , t>x>0 

sm( uid) 
sin( 142 1 ) 


(3.21) 
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Ax 

where Rtm 


cos (a -• J-) sin{ y) t * ' 
Rtm cos v<i X Hin( ky y ) r * ' 
cos(t;},<i) 

COS( V2 t ) 


u 

t 


r ■ t 
r ' (} 


(3.22) 


LSE normal modes propagating in the +: direction wit! dt‘rivo(i from (3.21) 
with —jz dependence using (3. 2-3.3). LSM norma! nuKie.s pntpagating in the 
+z direction will be derived from (3.22) with : dependtuicc* using (3. 5-3.6). 
If the mth mode in the +z direction is represented as, 


Em = {St +£,z)e-'^^ ^ r- (3.23) 

H,n = {% +n^z)e‘^' ^ fi,le (3.24) 

then the mode propagating in the -2 direction is given by. 

Em = {St -€xz) (3.26) 

Hm = +nxz)e^-^^ = . (3.26) 


The caligraphic letters 6t, Hz etc. will be used to denote the normal mode field 
components without the dependence. 


3.2.3 Orthogonality Relations 

The orthogonality relations between the normal modes will b<i introduced at 
this stage. General orthogonality relations between modes in uniform cyiindri- 
cal guides with arbitrary filling is discussed in Bresler (39). For an inhomo- 

geneously filled waveguide, the following (restricted) form of orthogonality is 
available. That is, 
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[ Em XH^-Z dS = 0 . (3.27) 

Js 

Em is the electric field of the mth normal mode and Hn is the magnetic field 
of the nth normal mode. If E and H belong to distinct sets such as LSE and 
LSM, then the above relation is valid for all m and n. However, if both the 
vectors belong to the same set say LSE, the relation is valid only when m ^ n. 
Since the integrand is a scalar triple product involving z, and the integration 
domain is the transverse cross-section, the above relation can be reduced to one 
involving the transverse components only. That is, 

j^£t,mxnt^n-zdS = 0. (3.28) 


3.2.4 Derivation of Green’s Functions 

We shall adopt the notion of Green’s function as the field due to an elementary 
source (Dirac delta). For the analysis of Chapter 2, it is necessary to consider 
3 source types namely Jy , X and My. It is also necessary to consider 3 field 
types caused by the aforesaid sources, namely Ey , E, and Hy. This amormts to 
9 Green’s functions which we shall indicate as Ey{Jy) , Hy{Jf) etc. The rule for 
getting the fields from the soxirces will be (3.29) without any scale factors such 
as —jbJix because of the above definition. Thus, 

Ey = j^ Ey[Jy) JydS' . (3.29) 

Let the field radiated by the elementary source of concerned type at {x' , y', z') 
be represented as a series of LSE and LSM modes, 


E = a„ 


^ti,LSE ^ 


-j(z-z') 


-h 


E 

n 


•p Hr 

^n,LSM ® 


-7(W) 


(3.30) 


f 
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H = O.m'^ m,LSB ® 

m 


n 


(3.31) 


for z> z' and, 


E = XI ^m^m,LSB ‘ ^ + X 

m " 


J.SM 


ISM ' 


* • ( 


(3.32) 

(3.33) 


for z < z'. 



Figure 3.2: Radiating sources inside waveguide. 

The coefficients am, bn, a'^, 6', can be evaluated by applying reciprocity the- 
orem to a suitably defined volume. Fig. 3.2 shows the volume containing the 
radiating sources to which reciprocity theorem is applied. Nolo that this region 
is bounded by the surfaces Si, S 2 and the waveguide' walla. 

Consider the following pair of source-field situations: 

• The source pair J, M producing the fields E, H as defined in (3.30-3.33). 

• The source free solutions En, H„, which corresponds to one of the normal 
modes in ±z direction. 

The reciprocity relation between this pair of fields, (E,„H„) and (E.H) is 
given by, 
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£ (e X H„ - E„ X h) - ndS = ( J ■ E„ - M • Hn) • (3.34) 

Note that (E„,H„) satisfies source free Maxwell’s equations while (E.H) 
satisfies Maxwell’s equations with the concerned source type. That is, 

V X E = -jw/xH - M (3.35) 

VxH = +jujeE + 3 (3.36) 

The reciprocity relation reduces to the following form on applying to this 
volume and for this pair of fields. 



h h 


y^(j-E„ - M-H„) dy (3.37) 

One can choose E„, H„ as LSE/LSM mode fields with 2 choices on the direc- 
tion of propagation. By making use of the orthogonality relation (3.28), and the 
4 choices on En, H„, we have the following results. 

• E„ , H„ are the LSE fields propagating in +z direction. 

h = 0 

h = +2a'^PLSEe-'^^' 

• En , H„ are the LSE fields propagating in -z direction. 

h = 0 

h = -2am Else 

• E„ , H„ are the LSM fields propagating in +z direction. 
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h = 0 

h = 

• E„ , H„ are the LSM fields propagating in -c direction. 


h = 0 
h = 

where 

PlSE = l^[£t,LSE X %.LSt!)’7. (IS 
^LSM = J^(£t,LS.M X ’Pt.LSM)'7 tlS . 

The coefficients a,n,bn,a'^ and 6^, can be evaluated by computing the right- 
hand-side of (3.37) in the above 4 cases, leading to the following n^sulta. 


• Jy as source. 

_ , !/) 

tjjyn mmm mmm, wmmt ■ 

~2 P,,K 

fj _ ^ 'Jf' ^ 

^'m ~ +«m 

K = +hn 

• Jr as source. 

n - ^z,LSE{x\y') 

am — ii—l 

+2 Pise 

bn ~ y') 

, +2 Plsm 

Qm = 

b'n = -bn 


(3.38) 

(3.39) 


(3.40) i 

(3.41) I' 


‘V) 




7('i 



I, ' 


hi\ r 
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• My as source. 


a = '^y^Lssix' , y') 

-2Plse 

h _ Py,LSM{x' ,y') 

Uji — — _ 


(3.42) 

(3.43) 


(3.38-3.43) can be substituted in (3.30-3.31) to obtain the following 9 Green’s 
functions. 

Ey{Jy) = sinMi(a-r’) cos kyy' smui(a — x) cos kyy 

+ 21 sin vi {a -x') cos kyy' sin vi{a - x) cos kyy e I 

P^iPj) = ^£--(Jy) sinui(a-r') cos kyy' smui(a — x) siakyy ^ 0 

+ 22 ^E^iJy) sinvi(a-x') cos kyy' sin^;i(a-r) sin e I I 0 

Pyi'^y) — Y2 sinui(a — r') cos kyy' COSW 2 X sin kyy ^ 0 

■+■ 22 {o- — x') COS kyy' cos v^x sin kyy e I I 0 

Py{Jz) = Y1 sinui(a-r') sinkyy' sintii(a-r) cos e ^ I 0 

+ 22 Eeaj.) sint;i(a-r') sinkyy' sintii(a-a:) cos /c^y e I I 0 

Ez(Jz) = Y2 Ee^{J:,) sinui(a — r') sinkyy' sm.ui{a — x) sinkyy 

+ 22 (Jx ) (a - ) sin kyy' sin vi (a - x) sin kyy e I I 

Hy{Jz) = Y2, sin Ml (a — a:') sinkyy' coSM 2 a; sinkyy 
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+ Yl sinr;i(a - re') sin kyp' cosu 2 -r sin kytjc 


Oml-- •'! 




(-) 


E,{M,) = ECJ, costi2a:' sin^j,y' sinui(a-2;) coskyi/f 

+ Y1 ^EyiMy) COSV 2 X' sin kyp' sinvi{a - x) cosky pc '’'"‘h’ ''I (-) 

E^{My) = C%^{My) cosu 2 a:' sinkyp' sinui(a-a:) sinA:j///(' "•■■I' '"'I 

+ ECS cosz; 2 x' sinkyp' sint;i(a-a:) sin A:,,;/ r I ■ •'! 


Hy{My) = Y ^HyiMy) U 2 x' sIh A'y^' 005^2:1’ sin kype 

+ E ^Ey(My) COSV2X’ sin kyp' cost;2a; sinkype ' 


■V I • ' I 


■Im I-’ “S' ! 


where, 


^Ey{Jy) 


^Ey{Jy) 


^kiJy) 


^B 4 Jy) 


^HyiJy) 


/^m 

^Hy{Jy) 




/^m 

^E,{J,) 


^Hy{Jz) 
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■2 Plse ^^0 
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“2 Plsm 
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■Vi ky'Yr, 
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Of the 9 waveguide Green’s functions listed above, the source point {x\ y\z') 
and the field point {x, y, z) are restricted to be in either the dielectric region 
or the air region of the inhomogeneously filled waveguide. Green’s functions 
of the B{J) t 5 rpe has both source and field point on the air side, whereas the 
type H (M) has both points on the dielectric side. Green’s functions of the H ( J) 
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type has source point on the air side and the field point on the dielectric side, 
whereas the E{M) type has it vice versa. 

It is possible to (partially) verify the Green’s functions, by reducing them to 
the case of homogeneously filled guide. Therefore, substituting Cr = 1 or t = 0 
should lead to Ey{My) = 0 and Hy{Jy) = 0. All the summations on Green’s 
functions are double indexed from m = 1, n = 1 (to m = +oo, n = +oo), except 
for the LSE part of Ey{Jy) Green’s function. For this case, the sum starts from 
m = 1, n = 0 through the dominant mode. In the above expressions 7 e and jrn 
represent the LSE and LSM propagation constants. 

The Green’s function E,{JS} as given above is valid only if the field point 
does not coincide with the source point. In cases where this is not true, it is 
necessary to add a Dirac delta function to the above expansion. Therefore, the 
Ez{Jz) valid at all points is given by [40], 

Ez{Jz) = Ez{Jz) + — (3.44) 

jujeo 

This result has significant implications throughout the present work and 
Section 3.4 is dedicated entirely to this Green’s function. 


3.3 LSE and LSM Propagation Constants 

3.3.1 Introduction 

The LSE/ LSM eigenvalues U 2 , vi and V 2 appear throughout the Green’s 
functions at various places. It is necessary to solve these values at every 
frequency in order to evaluate the Green’s functions. This section deals with 
the construction of robust algorithms for solving these roots. 

The set {ui,U 2 ) which satisfy (3.17-3.18) simultaneously are defined as the 
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LSE roots, whereas the set (dj, v^) which satisfy (3.19—3.20) simultaneously are 
defined as the LSM roots. We make the substitutions x = U 2 t, y = md in (3.17- 

3.18) and x = ^ 2 ^, y = Vid in (3.19—3.20). The LSE/LSM defining equations gets 
transformed to, 


tan X tan y 

X ~ y 

LSE family 

(3.45) 

rtanr = — — ytany 

T 

LSM family 

(3.46) 

y2 

.2 _ „2 

^2 ~ ^ 

Hyperbola 

(3.47) 


where 

_ 

^ ~ 1 

c — /cq t y (-t — 1 • 


The solutions of (3.45) and (3.46) are disjoint family of curves which are 
sketched in Fig. 3.3. (3.47) is a hyperbola whose positive asymptote is the line 
y — rx. "We shall investigate the salient properties of the 3 set of curves, such 
as the trajectory, slopes etc. in order to construct a solution scheme. 

We shall make use of a modified form of the algorithm rtsafe from [41]. This 
uses both Newton-Raphson and bisection so that it is necessary to provide both 
a starting value and an interval for the root. The algorithm makes use of 
Newton-Raphson for a fixed number of times. It then switches to bisection if 
the root wanders off the interval, or if Newton-Raphson is unable to achieve 
the accuracy within a certain niimber of iterations. This scheme is necessary 
for the present problem because the functions involved are multiple valued and 
the solution may fall into another possible region. The idea of the following 
discussion is to provide good starting values and safe bounds for the roots so 
that solution scheme is robust. 
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I Lj U Li 1 L 1 1 

0 TC 271: 371 ^ 

Figure 3.3: LSE/LSM family of curves 


3.3.2 Concept of a Root Square 

With Reference to Fig. 3.3, the region of interest, namely r > 0 and y > 0 can 
be divided into squares of size tt x tt centered at (?Ti7r, utt) where rn. n = 0, 1, . . .. 
Each such square will be defined as a root square. The squares with m, n > 0 
will be called proper root squares, while squares with m = 0 or 7z = 0 will be 
called improper root squares [See Fig. 3.4]. A root square centered at {init. nix) 
will be referred as "(m, n)th root square". 

For each root square we identify 3 points, namely bottom right .4 (.zy. /;(,), 
center B{xc, Uc) and top left C{xi,yt). For improper root squares, either C or ,4 
will be outside the region of interest depending on whether m = 0 or ‘ii = 0. 

The LSE/LSM curves pass diagonally across each root square such that it 
passes through A, B and C. It is this property that prompts us to introduce the 
notion of a root square. Note that in (m, 0) root squares, the LSE curve does 
not pass through B. This is because tanx/x = 1 when y = 0. Therefore, in 
improper root squares LSE curves does not pass through B. 

We can view the root finding problem as one of identifying the sequence of 
root squares through which the hyperbola passes. Each such root square will 
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n% 


A{Xr,y^) 


Figure 3.4: Definition of a root square. 


give two roots, one for LSE and the other for LSM. For example, in Fig. 3.3 we 
can identify the sequence of root squares where the hyperbola intercepts the 
LSEJLSM curves. The root sequence for this example is, 

( 1 . 1 )— ^( 1 . 2 )— 

We have thus completed the first step of the root finding. That is, to move 
along the hyperbola from its apex and identify the root square sequence. The 
remaining two steps of the algorithm are now dealt with — finding good starting 
values and safe bounds, so that a root finding algorithm can start polishing a 
root. 


3.3.3 Choosing a Starting Value 

To find a good starting value in a root square, we proceed as follows. Since 
the LSE/ LSM curves pass through A, B and C, the intersection between the 
diagonal AC and the hyperbola can be taken as a starting value. It is easy to 
find t.hi.c! point xq , since it involves the solution of a quadratic. That is, 
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Xo = 


^Jiyb + Xr)^ + - 1 ) • 


{Vb + Xr)2 + (rc)2 

5-2 _ 1 


(Vb + Xr) 


(3.48) 


There is yet another reason for opting for this starting value. We will show 
that the slope of the LSE curves at A, B and C tends to —1 in the limit. 
Therefore, in the limit, in a root square with the hyperbola, the LSE curves are 
close to the diagonal .4C. 

To prove this, we estimate the slope of LSE curves at .4, B and C using the 
following approximations of tanx. 


1. If X is around rmr, then, tanx ~ x — rmr for m = 0. 1, . . . 

2. If X is around m7r/2, then, tanx ~ ±l/| X — (m7r/2)j for m =■ 1,3, ... 
Using these results in (3.45—3.46), we get 


dx 

dy 

dx 

dx 



(3.49) 

(3.50) 


(3.51) 


Note that yb/xr, yc/xc and yt/xi tend to r, since y = rx is the asymptote of the 
hyperbola. Therefore, the slope at the three points tend to -1. These results 
do not hold good for LSM curves. However, experience with the solution has 
suggested that this is a good enough starting value. 

Therefore LSE! LSM roots can be solved by applying Newton-Raphson method 
to (3.52-3.53) with xo as the starting value. That is, 


f[x) = + 5 -!^^ = 0 with y = r \Jx^ - (3.52) 

^ y > 

/(x) = X tanx + ^ytany = 0 with y = r\/x"^ - 


(3.53) 



3.3 L SE a nd LSM Propagation Constants 


75 


It is also possible to solve an LSA'I root using an LSE root as the starting 
value. 

3.3.4 Choosing a Root Bound 



Figure 3.5: Passage of the hyperbola through a root square. 

We define a root bound [;ri,.r 2 ) as an interval on the .r axis such that the f{x) 
has a simple zero within it. We shall soon see that the bounds which we are 
about to compute has the nice property that /(.tj) < 0 and f{x 2 ) > 0. This will 
be useful in bisection algorithm to decide how to update the root interval. 

We introduce two functions x{y) and y{x) to facilitate the discussion of root 
bounds. Note that these functions give just the coordinates on the hyperbola. 

y{x) — ryx^ — 
x{y) = \l{y/rf + c^ 

Fig. 3.5 shows the 4 possible cases by which a hyperbola can pass through 
a root square. In each case, it is easy to identify a root bound. For example, 
if it is case 1, we can choose the following root bound [.r;, a:(y()]. For each case. 
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we can also identify the next^ root square to which the hyperbola gets in. For 
example, in both case 2 and case 4 the next root square is (j??, + 1, n). Table 3.1 
summarizes the essential information. Note that case 3 and case 4 do not occur 
in (m, 0) root squares while case 1 and case 2 do not occur in (0, 7i) root squares. 

Table 3.1: Look up table for the root finding scheme 


No 

Case 

Root Bound 

Next root square 

1 

yb < y{xi) <yt & xi < x{yt) < Xr 


( 771. , 76 + 1 ) 

2 

Vb < y{xi) <yt & 2/6 < y{xr) < yt 

[.^■f ,-T, ] 

{m + 1 , 76 ) 

3 

Xl < x(yb) <Xr & Xl < x{yt) < Xr 

[•T(y6),;r(2/()] 

(771 , 76 4- 1 ) 

4 

xi<x{yb)<Xr & yb<y{x,)<yi 

[x{yb).Xr] 

(711 + 1 , 71 ) 


3.3.5 Modified Newton-Raphson-Bisection Algorithm 

The root finding scheme is based on (3.52-3.53) with ro as the starting value. 
We also make use of Table 3.1 to decide the root bound and the next root square. 
Once we are in a root square, the following algorithm can be used to polish the 
root. 


Algorithm rtsafefn , X 2 , xo, f , x) 
parameteriNR = 10, BS = 50, err = l.Od - 12) 
c Solution of f{x) = 0 by combined Newton-Raphson and 

c Bisection, [ri, .r^] is the root bound, xq is the 

^It is assumed that we are currently at the (m, 7 i)th root square 
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c starting value. NR, BS are the maximum allowed number 

c of iterations for Newton-Raphson and Bisection respectively, 

c err is the maximum allowable error in |/(x)| 

T. = Xq 

do i = 1, BS 
f = f{^) 

if ( 1 /! < srr) return 

dif = f/f 

Xnew= X - dif 

if ( {Xnew > ^ 2 ) or (rCnew < xf) Or {i > NR)) then 
if (/ < 0) then 

Xi=X 

else 

X2= X 

end if 

Xnew ” (^1 d" ^ 2 )/^ 

end if 

X === Xnew 

end do 

stop 

end 


I 

This discussion has not covered every exception that occurs while solving the 
LS E/ LSM roots. However an implementation based on these ideas has been 
quite satisfactory. 
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3.4 Ez{Jz) Green’s Function as a Distribution 

3.4.1 Introduction 

Green’s function is the solution of a Partial Differential Equation (PDE) with 
the right-hand-side replaced by a Dirac delta function [42] . The solution is to 
be interpreted in the sense of distributions. The rule for getting the solution 
for an arbitrary right-hand-side, involves the integral of the product of the 
Green’s function with the right-hand-side, over the source domain. We shall 
refer to this as the transformation rule. Note that, a Green’s function is a 
generalized function, whereas the solution obtained after the transformation 
rule is a function in the classical sense.^ The dimension of the Dirac delta 
used in the derivation of Green’s function implicitly dictates the number of 
quadratures to be performed in the transformation rule. Supposing, we deal 
with current sources and fields as in EM, a 3-D delta dictates a volume integral 
in the transformation rule. 

If the source domain does not possess one of the linear dimensions, as an ap- 
proximation to the actual geometry (or as an independent assumption by itself), 
it is possible to leave out the quadrature in that direction. However, the solu- 
tion so obtained will continue to be a distribution since the legitimate number 
of quadratures have not been performed. In a boundary value problem (BVP), 
where we deal with incident and scattered fields together, the presence of a 
generalized function (scattered field) leads to awkward problems. A general- 
ized function amidst classical functions demands sophisticated interpretations 
to make sense. 

In the present section, this problem is illustrated considering the evaluation 
the E~ field from E,{Jf) Green’s function. This problem is a part of the evalu- 
ation of the RSEj; submatrix in Chapter 2, Therefore, this section serves the 


^The phrase "to make sense" will be used to mean "to make sense as a classical function". 



3.4 E.{J.) Green’s Function as a Distribution 


79 


dual purpose of an illustration as well as an elaboration of the MoM matrix 
evaluation of Chapter 2. 


3.4.2 Problem Identification 

Consider a slightly modified definition of the Ez[ ) Green's function based on 
(3.44). 


Ez ( •/- ) == Etb sinui{a - x' ) sin kyy' smui{a-x) sin kyy e I I 

" ' ^ ----- - - ^ 

Bte 


where 


"*■ sin ri (a - x' ) sin kyy' sini;i(a-x) sin kyy e 

. ... V 

Btm 

+ 5(r - r') (3.54) 


Ctb = = 

C ’ 

TM - <-■£,( J,) - 


1 fef 

+‘^Plse Vq) 

1 / Vllm V 

+2Plsm 


The Green’s function consists of two double series corresponding to LSE and 
LSM normal modes and an explicit Dirac delta term. The presence of the 
delta function definitely makes'^ (3.54) a distribution. Being a distribution, the 
Green’s function need not make sense, so that both the series in (3.54) need 
not be convergent. However, we require the E^ field to make sense after the 
transformation rule between (3.54) and an arbitrary current. Therefore, 
each of the summation in (3.54) should be convergent after performing the 
quadratures in the transformation rule. However, this does not happen if one 


*An expression need not contain a Dirac delta to qualify as a distribution! 
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of the quadratures over the source domain are dropped. 

For the printed strip problem of Chapter 2, we neglect the thickness of the 
strip and regard it as a domain, 


S'j^ = ^ ] X [ —wi , +wi ] at r' = i . 

Therefore, we leave the quadrature in the x' direction while performing the 
transformation rule. The simpHfied transformation rule is given by. 


E- 


Jz 


[ E,{J,)J,dS' 
Js' 

E Ctb I 


Bte e' 


■Ic |C- 


' I J, dS' + ^2 Bt.m [ Bxm 6 

ds' 


-7m |z-.'| 


Else 


Bls: 


M 


+ K Six - x') 


(3.55) 


where 


K = 


-1 


[ S{u — u') 5{v — v') Jz dS' 

S n 


jweo JS'^ 

= —J.iu,v) 
juea 

The following observations can be made regarding the terms of (3.55). 


(3.56) 


1. Ez\j^ field contains a one dimensional delta since the quadrature in x' is 
left undone. 


2. The series S^se is convergent. 

3. The series Slsm is divergent. Specifically, the partial sums of Sr^sM in r 7 i 
diverges. 

The term <5(r - x') makes the E^ field a distribution unlike its function 
counterparts Ez\jy etc. We cannot handle a situation where we treat 
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distributions and functions on th© sam© catcg'oiy. Tb© convsrgsnc© or tli© 
divergence of the terms Sise or Sism mre difficult to prove. We do not provide 
any rigorous convergence proofs in this regard, but outline the results to be 
used to prove such statements. 

Each series in (3.55) is double indexed with respect to m and n. The conver- 
gence properties are primarily governed by the 0 dependence^ of the coefficient 
part Ctk and the integral Btb Jz dS'^ on m and n. The function 

Bte is bounded and is composed of four trigonometric functions. Of the four 
trigonometric functions, only tr( ky y' ) participates in the integral. It is possible 
to get an approximate idea about the series by looking at the growth rates of 
the coefficient and the integral. At large values of m, there is also the result 
[33] that (/i, t'l, U 2 and V 2 tend to — m'Kl a. Therefore, by using the results in 
page 69, we have 


Cte 

1 ...... <5n 




1-2 

(3.57) 

2 bSte 

✓ 

Je 

(nf- 

*5) 



bounded part 






Ctm 

_ 1 1 2 



I'm 


(3.58) 

2 juJiQ bStm 

N 

vl 

- kl 


boxinded part 

The quadratures are boimded by the reciprocal of the coefficients so that, 


[ tT{kyy') Jz dS' 

J 


0 


ky'J ^ 


for ze S'ji- (3.59) 


Using the results in (3.57-3.59) in (3.55), we get 

E ^ 4 

m, n I 


5lf / = Oig), then, there exists suitable K such that. |/(a:)l < K g{x) for aU sufdciently large 
values of x. 
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< K' E 


sin ky y 


i'2 i- 

m,n 


= (? I) (? 


sinfcj,yl' 


The series En I sin kyy\/ky is not absolutely convergent, so that this upper 
bound for is divergent. We cannot therefore conclude that Sise is conver- 

gent. It is necessary to derive very tight upper bounds to prove the convergence 
of such series. However, for specific cases it is possible to prove the convergence 
of Slse- 



Figure3.6; Partial sums of the LS’.t/ partof 1) matrix entry under various 

conditions 


The series Slsm can be shown to be divergent without resorting to any elab- 
orate convergence proofs. Note that Ctm grows as 7 , whereas the integral in 
(3.59) is just 0(1/ (/cj, 7 )). Therefore, the series for Slsm does not have any terms 
in the denominator which increase with m. We thus have a series in m whose 
mth term does not tend to zero. We can definitely make the statement that S^sm 
is divergent in m. Fig. 3.6 shows the partial sum of the LSM part of RSE^j 
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with respect to vl. The curve in case(a) indicates an unbounded partial sum. 
We will be using this matrix entry for further illustrations, since Slsm and this 
entry are identical insofar as the growth with respect to vn is concerned. 

We are thus faced with a problem of interpreting a diverging LS M series 
and an explicit Dirac delta in a supposedly classical expression (3.55). It is 
intuitively clear that the LSM series and the Dirac delta should be handled 
together to make the expression (3.55) meaningful. This is done by representing 
the delta function in terms of a valid complete orthonormal set. 


3.4.3 Representation of Dirac Delta 

We shall use the following result for representing Dirac Delta. 

Theorem 3.1 If { Cm{x ) } is a complete orthonormal set, then 

6{x - x') = Yl ^rn{x') ■ 

m 

This result is known as the completeness expansion and is referred loosely 
as the expansion of singularity. The word expansion is misleading, since the 
right-hand-side of the above equality is not a series in the conventional sense. 
The series does not converge so that it is to be understood only in the context 
of doing a quadrature. That is, if f{x) is any function then 

m 

is a valid expression. We shall refer to this as the representation of singularity 
rather than the expansion of singularity . 

We can choose between LSE or LS M eigenfunctions for representing S(x — x'). 
Incidentally, the LSE eigenfunctions are orthogonal to begin with, whereas 
LSM eigenfunctions are not. Therefore, we opt for the LSE eigenfunctions. 
The LSE eigenfunctions are defined by { fm }, where 
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_ f smui(a-x), a>x>t 
jm — \ (u.oO) 

[ Rte sin U 2 X , t>x>0 

Note that fm satisfies all the boiindary conditions to be satisfied by the 
component in the x direction. That is, 

1. /m(0) = fm{a) = 0 

2. fmi^) is continuous at x = t 

3. f'^{x) is continuous at x = t 

The orthogonality of the set { fm} can be derived from the general orthogonality 
statement in Equation 3.28. We can therefore use, 

5{x - x') = Y. Tf ~ 7T fm{x') (3.61) 

m \Jm 1 Jm) 

in (3.55) and study the convergence of the partial sums in m. This is shown as 
case(b) in Fig. 3.6. It is to be noted that this series is divergent just as in the 
previous case where the series S^sm was taken alone. 

We explain the precise reason for this apparently strange result. If the 
expression (3.55) were to make sense, it is necessary that the series for 
cancels the delta function. It then follows that the series for Si^sai contains an 
LSM represented singularity. Therefore, what we are attempting is to cancel 
the LSM represented singularity in Sism by combining it with an appropriately 
represented 5{x ~ x'). If we do so by representing S{x - x') as an LSE series, 
then we are trying to subtract two different representations of delta function. 
This however is a meaningless operation. We will clarify this point as follows. 
Suppose {em} and { } are two (different) complete orthonormal sets, then 

5{x - x') = Y ^rn(x) em(x') (3.62) 
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^{t - x') = 22 hm{x) hm{x') (3.63) 

m 

Y1 Cm(x') - hm{x) h^{x') = 0 (3.64) 

are true statements in the sense of distributions. However, none of the state- 
ments makes sense if the summations are carried out literally. In fact, evalu- 
ating (3.64) and hoping to get zero is as meaningless as evaluating any of the 
expressions in (3.62-3.63). 

It is therefore mandatory to represent the delta function as an orthonor- 
malized LSM series so that a term wise subtraction takes place. The LSM 
eigenfunctions are however not orthogonal. Therefore, it is necessary to or- 
thonormalize the system using the Gram-Schmidt procedure. 


3.4.4 Gram-Schmidt Procedure 

We begin with the LSM eigenfunctions suitable for representing the x variation 
of the E: component. This set {fm}is defined such that. 


fm ~ 

dm 7^ 


a'mi ^) ' a>x>t 

t>x>0 

(3.65) 

cost)i(a - r) , a> X > t 

Rim COSVzX, t>x>Q 

(3.66) 


This choice satisfies all the boimdary and continuity conditions to be satisfied 
by the E, component in the x direction. That is 


1. /m(0) = fm{a) = 0 

2. fm (a:) is continuous a.tx = t 

3. amix) is continuous a.tx = t 
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The set { Om } can be shown to be orthogonal in the following inner product 

(/ , 5 ) = — / f gdx + ! f gdx. (3.67) 

67 . Jo Jt 

We shall define the inner product in { /,„ } as 


{f , 9) = ^ 3 dx + f gdx (3.68) 

so that the orthogonality properties of { } can be used to simplify several 
steps while evaluating (/^ , /n)- Gram-Schmidt procedure requires (/^ , /„) for 
all values of m and n. These can be computed by making use of (3.65—3.66), 
the orthogonality of { } and the 3 properties enumerated above. (/^^ , for 

m ^ n and m = n is as follows. 


{fm , fn) 


m^n 


A 

B 



_ V2,m A - V2,n B 

4,m - 

= sinV2,mt COSV2,nt 


= COSV2,mt SmV2,nt 


(/m , fm) = 

2 -^L 

^2,m 

Zr 

C + 

C = 


1 

sm2t^2,Tni' 


2 1 

2 ^2, mi , 

D = 

u 

f 

1 + 

sin d 


2 \ 



(3.69) 


(3.70) 


The parameters of V 2 ,m , vi^n etc. refers to V 2 , vi etc. evaluated for the corre- 
sponding index. In the next subsection we shall construct the algorithm for 
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implementing the Gram-Schmidt procedure. 


3.4.5 Algorithm Implementation 

Given a finite or countably infinite sequence { fi} of linearly independent func- 
tions in an inner product space, Gram-Schmidt procedure generates both an 
orthogonal sequence { } and an orthonormalized sequence { Cj }. This can be 

represented symbolically as 

^ M ^ M 

linearly independent orthogonal orthonormal 

It is preferable to work with {g,} so that, we shall be using the latter of the 
following representations of the delta function. 

6{x - x') = ^ ei{x) e^{x') (3.71) 

i 

5{x - x') = Yj / ^-T 9^{^) 9i{x') (3.72) 

i \9i ) 9i) 

Gram-Schmidt procedure [42) generates {g^} and { e,- } in the following se- 
quence 


91 = /i 

92 = f2 - {h . ei) ei 

93 = h - (/a . «i) “ (/a > ^ 2 ) ez 

fc'-l 

9k = fk ~ {fk^ 6p) 

p=l 


(3.73) 


Replacing with gp/\\gp\\ in (3.73), we can rewrite the sequence as, 
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91 = fi 

(3.74) 

k^l 


Qk — fk + Y2 ^>cp 9p ^ k > 1 

(3.75) 

p=i 


{fk 1 9p) 1. ^ „ 

O^kp / \ 1 k p . 

{9p , 9p) 

(3.76) 


Despite being constructive in its very definition, Gram-Schmidt procedure is 
not convenient in the form presented in (3.74—3.76). In (3.75), the function gk is 
defined recursively in terms of Qp’s which is difficult to code. We need to define 
Qk directly as a linear combination of /^s; that is, 

fc-i 

gk = /fc + E fp ■ O-^V) 

We can state the objective of the algorithm as follows: 

Objective 3.1 Given {/,} and as in (3.65), (3.69-3.70) for 

i= 1, 2 N; compute Pkp , N > k > p and {gk , gk) , k < N so that, 


9k ~ fk^ Pkp fp 

(3.78) 

p=l 


6{x ^') = . 9z{^) 

1=1 \9i ) 9i} 

(3.79) 


The algorithm is presented recursively. Assume that we have completed the 
sequence up to the (k - l)th stage. That is, Pkp and {gk , gk) are known up 
to k -1. We shall indicate the operations for getting these for kth. stage. With 
reference to (3.76), akp can be computed from {fk , pp) and {pp, pp). {fk , pp) is 
computed using, 

p-i 

{fk , gp) = {fk , /p) + E Ppm {fk , fm) (3.80) 

m=l 

whereas {pp , Pp) for p < k - 1 is already available. Having computed Okp, the 
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term {gu , gk) can be computed using 


fc-i 


{9k . gk) ~ ifk , Jk) - Y. ^Ip {g? , 9p) ■ 

P=1 


(3.81) 


The evaluation of .4,, can be conveniently carried out once we adopt the matrix 
notation. Represent (3.75) and (3.77) in matrix form so that, 


and 


9k 


fk + [o-ki Oiki ctfcs ••• akk~i 


91 

92 

93 

l9k-l 


(3.82) 


9i 


■ 1 

0 0 

... O' 


h 

92 


‘hi 

1 0 

... 0 


f2 

93 

rs 

ihi 

P 32 1 

... 0 


fs 

9k-l. 


I0k-ll 

0k~12 Pk-13 

... 1 


. fk-1 . 

9k = 

= fk + [Pki 

Pk2 PkZ ' 


r 

h 1 


/2 
/3 

fk-] 

Substituting (3.83) in (3.82) and comparing it with (3.84) we get, 


(3.83) 


(3.84) 


A;-l 


0kl = + Y 

171=2 
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fc-i 

Pk2 = OCk2 + X) Qikm (3m2 

7n=3 

k -1 

Pkp ^ (^kp "1” ^km Pmp (3,85) 

m=p+l 

’riie method of computing akp, Pkp and {gk , gk) for the kth stage from the 
k - 1 th stage has been illustrated. The following implementation is based on 
(3.76), (3.80-3.81) and (3.85). 

Algorithm gram_schmidt(A^) 
c The values (/i , fj) is available for J < N 

c 

c Algorithm start up 

{gi 1 gi) = (/i , fi) 

0^21 = —{f2,fi)/{gi,gi) 

P21 — 0:21 

{g2 ) 52) = if 2 , A) + 2 P21 (/2 , /i) + P21 {gi , gi) 

c Algorithm main body 

do k = 3, N 
do p = 1, k — 1 
sum = 0 
do m = 1, p — 1 

sum = sum + ppm {fk , fm) 
end do 

sum = {fk , fp) + sum 
akp = -sum/{gp, gp) 

end do 


sum = 0 
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do /j 1 , A- - 1 

sum - sum + al{gj,,gp) 

end do 

(//A . fik) - {fk, Ik) - sum 

do j) = 1, A: - 1 

sum = 0 

do in =/;+!, A: - 1 
sum = sum + akm 0mp 

end do 

‘k;, - (ikp + sum 

end do 
end do 
stop 
end 

In actual implementation it is necessary to use a 2-D array beta (N, N) for 
implementing ik-p- However, since k > p both the diagonal and the upper 
triangular part of beta(k,p) remains unused. We can use this space for 
storing ru,, and (.% , gk). Therefore, in the code it is possible to replace akp by 
beta (p , k) and {gk , gk) by beta (k, k) . Thus, the storage requirenaent of the 
code is exactly N^. 


3.4.6 Results and Discussion 

Cases (c) and (d) of Fig. 3.6 illustrate the result of using the LSM representation 
of delta function. It can be seen that the partial sums are comparatively well 
behaved unlike the previous two cases. Cases (c) and (d) will be used to bring 
out the subtle issues involved in subtracting the Dirac delta numerically. Note 
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that, we are trying to subtract two series which by themselves are divergent. 
Therefore, even a small mismatch between the series will add up and starts 
showing up at large values of m. Therefore, subtracting two divergent sums is 
an intrinsically unstable operation. This can be observed on close observation 
of the curves (c) and (d) of Fig. 3.6. The partial sums starts growing linearly 
at large values of m. The linear growth of partial sums is an undesired result 
since it makes the partial sums divergent. However, we will be making the 
observation that the linear growth arises from imperfect cancelation of Dirac 
delta. Of the two cases case(d) exhibits lesser slope than case(c). These cases 
will be studied carefully to validate the above observation. 

The two cases c and d arise from the different representation of Dirac delta 
(i(r — r') in (3.54). Note that, both 


5(r - r') = 5{x - x') 5{u - u') 6{v — v') (3.86) 

and 

(i(r - r') = 6{x - x') 5{y - y') 5{z - z') 

f 2 ^ \ 

= (5(x -x') sinkyXj sin ky y' j 6{z - z') (3.87) 

are valid forms of (i(r — F). If we use (3.86) in (3.55), the term corresponding to 
the delta function is obtained as just A' 5{x - x'). However, if (3.87) is used in 
(3.55), we get a series representation of A'. Thus, the last term corresponding 
to 5{t - r') in (3.55) can take two alternate forms, namely 


— ^LSE + diSM + K S(x — x') 

K = —-r—Jdu,v) 
jweo 


or 


E, 


z\Jz 


(3.88) 

(3.89) 


— Else + Elsm + Kp/5{x — x') 


(3.90) 
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A-.v 


I 2 ^ f r \ 

'.i^o K y- (3.91) 

* fltsl \ H J 


The Keric‘.s .s'; su is a double summation in both m and n. It is not possible to 
carry out either of the summation in closed form. Supposing we sum N number 
of terms <ilong //;, then the embedded singularity in will be of the form, 


^LSM = S'l^sM ~ 5{x - x') . (3.92) 

In (3.92), we explicitly show Sism containing an embedded Dirac delta. If 
this representation of *S;^.s*a/ is used in (3.88) and (3.90), we get two alternate 
forms for E.]./, . That is, 


E. 


Slse + - A W 5{x - x') + { 


Slsm 


K5{x - x') 

V ' 

case (c) 

Ke S{x — x') 

' ’V 

case (d) 


(3.93) 


In ca.se(c) we use the closed form representation of 3(r — r'), whereas in case(d) 
we use the truncated N term approximation. It can be seen that the Dirac delta 
cancels off completely in caseld); whereas a residual {K - Kn) 6{x — x') remains 
in case(c). Therefore, the difference between the curves c and d arises from 
the residual Dirac delta. This is illustrated in Fig. 3.7, where we have plotted 
the difference between curves c and d and the residual delta side by side. Note 
that, the presence of the Dirac delta causes a linear growth for the partial sxims 
in m. Ideally speaking, we should have had no residual delta in case(d), and it 
should have been flat. In actual computations perfect cancelations do not occur, 
so that (3.92) may not be valid with K^. Therefore, by the above reasoning, 
the small yet definite growth in case(d) can also be attributed to the presence of 
uncanceled Dirac delta. It is not easy to safeguard against such errors, except 
by limiting the number of modes summed to a reasonable number. 
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Figure 3.7: Explaining the linear growth of cases c and d as the presence of resid- 
ual singularity. The partial sum on the right is due to the LSM represented 
(A" - A» 5{x - x') 


3.5 Half Space Green’s Function 


The expressions for implementing the slot exterior matrix SXT is discussed. 
Fig. 3.8 shows the exterior view of the slot. A centered (;/, z) coordinate system 
is used for analyzing this problem. The slot is assumed to be embedded in 
an infinite ground plane, so that we are interested in the magnetic field //„ at 
X = 0“ in the external half space. 

The half space problem is identical to free space if the ground plane is re- 
moved and the source strength is doubled. Consequently, the electric vector 
potential for the half space is given by, (Note: bold j denotes \/^). 




( 3 . 94 ) 
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2 /j 


z 


H 

Figure 3.8: Slot exterior view 


If /, 0 is the ith basis function on the slot, then we have My — -f^ 0 and the 
j/-th matrix entry is given by, 

+ (3.95) 

jue/j. \ dy^ J ^ 

SXT{j,i) = [Hyf.dS (3.96) 

Let G denote the canonical source point integral in (3.94) so that, SXT{j,i) is 
given by, 


G = G{y,z) = j 

' fi 0 dS' 

5' r 

(3.97) 



[G" + k^G) dS 

(3.98) 

constant part 




where G" denote Note that C?(y, z) inherits the symmetry of Therefore, 
G{y, z) is even in 2 , whereas it is even or odd in y depending on whether /, is 
even or odd. Consequently SXT{j, i) = 0 if |i - = odd. We shall therefore 
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assume that j and i are both even or odd. 

The fundamental problem in evaluating (3.98) is to get rid of the partial 
differentials within the integral sign. This can be done® by evaluating the y 
integral in (3.98) via. integration by parts. After doing the integration by parts 
and using the results fj{±ls) = 0 and /" = we get, 


and 



(3.99) 


SXT{j, t) 


where 



(3.100) 



— fj{+h) — fj{ — h)G{ — ls,z). 


(3.101) 


(3.100) can be evaluated provided G{y, z) is computable at any point. The 
evaluation of G cannot be carried out directly due to the presence of the 1//' 
singularity in the integrand. Though the 1/r singularity is integrable, it is 
necessary to reorganize the quadratures in (3.100) so that a standard algorithm 
can be used. To begin with, we separate into singular and nonsingular 

parts. That is. 


e 

r 



(3.102) 


®This idea is based on Elliott [7] 
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Let and denote the corresponding part of Zzd and arising from the 

1/r term. Define, 


Z'zd = [fj f -f^4>dS'dS (3.103) 

S J S' T 

and 

^id = 2/j(+L) f f dS' dz . (3.104) 

J z J T" 

We will reorder the sequence of partial quadratures in (3.103-3.104) such 
that, both contain a term Gq of identical structure. That is, 

ZL = / /, / - A dS' dS 
Js Js' r 

= [ fid> [ fj [ - dz dy dS' (3.105) 

Js’ Jy Jz r 

S ,,, ^ 

Gq 

Z'u = 2 /;(+/,) / f h<i>dS'dz 

= 2/;(+g J^^dzds' 

= 2/j(+g [ <j> f fi [ -dzdy' dz' (3.106) 

V ^ 

Go 


Both Z' 2 d and Z'lj, contains the generic integral Gq which can be evaluated as 
follows. We carry out the z quadrature analytically using the identity 


I 


dz 


z y/ Z^ + 


= la (^z + \/ z^ + 


and combine the quadratures over [— Z*, 0] and [0, +ls]. That is. 
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f -fjdzdy = 
Js r 


/ /j In 

y=0 


■+-ls 


9{y - y') 
9{y + y') 


dy, 


for j = 2. 4. . . . 


/ /j ln[g{y-y')-g{y + y')] dy , forj 

y=o 


1,3,... 

(3.107) 


where 


g{y) = [ {ws - z') + ^J {ws - z'Y +y^ • [ (u/j + z') + 'J{ws + + y^ ^ 


(3.108) 


The quadrature in (3.107) can be evaluated using any standard routine which 
can take a logarithmic singularity at the end points. The routine dOlanf (or 
dOlajf) from Nag [36] can be used. 


3.6 Method of Moment 

3.6.1 Introduction 

Moment method is definitely the most widely used numerical technique in 
EM. Despite its widespread use, there has been few serious studies on it as a 
numerical method. MoM as a technique does not quite qualify as a numerical 
method, since it does not provide any error estimates nor any convergence rates. 
As a method, MoM does not have sound theory and it is based on the straight 
implementation of a very naive idea. In EM, MoM probably works because of 
the extremely well behaved right-hand-sides and the not so smooth nature of 
the kernels. 
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The absence of any error/convergence estimates is justified since the method 
is applied to complicated operators. However, the sequence of approximate 
solutions do not improve forever if the number of basis functions are increased. 
In fact, it can be shown that the solution of an integral equation via. MoM is 
an ill-conditioned problem, so that, the sequence of approximate solutions first 
improves and then becomes unstable. This aspect of MoM, which we think is 
the most serious of all has never been brought to the limelight. 

The motivation for the present discussion is our experiences with MoM insta- 
bilities as well as the results from S. W. Lee [30], Raj Mittra [31] and M. Leroy 
[32]. This section is primarily devoted to certain problems, which illustrate the 
fact MoM is an intrinsically unstable method. In what follows, we discuss the 
theory and practice of MoM and interpret the instability of solutions from the 
point of view of regularization. 


3.6.2 Theory of Method of Moments 

Method of Moment solves a linear operator equation of the form Tx = y. 
Algorithmically, it gives a sequence of approximate solutions which converges 
to the actual solution. MoM can also solve coupled operator equations of the 
form, 


Til + Ti2 X2 — yi 

T21 Xi -t- r22 X2 = y2 ■ 

We define three spaces related to the operator T : Domain, Range and Null 
space. The domain of T is the set of all x for which Tx is defined. The problem 
Tx — y also stipulates certain boundary condition on x which generally makes 
the domain a subspace. Therefore, any element of the domain is such that Tx 
is defined and x satisfies the boundaiy conditions of the problem. The set of all 
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elements of the form Tx constitute the range of T. The null space is defined as 
the set of all x such that Tx = 0. 

Given Tx = y, it is necessary to worry about the existence and uniqueness of 
solutions. However, such discussions can be made only for concrete realizations 
of T such as a matrix or an integral operator. This topic will be postponed until 
we encounter an integral equation. The question of existence can be ignored on 
the argument that a solution exists on physical grounds. However, uniqueness 
is more serious since it enters indirectly while attempting to solve Tx = y. 
We therefore make the readers familiar with the four equivalent statements of 
uniqueness. 

• Tx = y has unique solution. 

• Tx — 0 X = 0. 


• 0 is not an eigenvalue of T. 

• The null space of T is empty. 

We call a set { } total in a subspace V, if, given e > 0 and any x G V, there 

exist Ji, I 2 . . ,/yv such that - 2 ;|| < e. Therefore, a total set can 

approximate any element of the space by a suitable finite linear combination. 

The approximate solution of Tr = y begins by identifying a set { e, } in the 
domain of T such that {Te^} is total in the range space of T. Therefore, given 


e > 0 and any y, there exist fy, fy, • ■ • , In such that /, T 

is, 


y 


< e. That 


T (y /. Si) - 


y 


i=i 


< e. 


Therefore is an approximate solution to Tx = y. We call it as the 

e-approximate solution ofTx — y. Thus, we seek the e-approximate solution of 
Tx = y as e tends to zero. The actual solution process differs from the above 
notion. Instead of fixing e and demanding a solution, we fix N and ask for the 
best solution. That is, given A', we seek a solution xn = such that, 

II Efli IzTe,- y\\ is minimum. Therefore, we choose the coefficients fy, / 2 , . . . , //v 
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appropriately such that the error in the range space is minimum. If we define 
6 A' as, 


N 

(N = min II Y' k - y 
-Jn Ii ^ 

then 

<^N ^ ^N+l > eAf+2 ^ ^0. 


This inequality sequence is valid’^ such that tends to zero as N tends to 
infinity. Therefore, the sequence converges to the solution of Tr = y. In its 
primitive form, the solution of Tr = y thus involves a search for the right set 
of coefficients /i. / 2 , — In- This is the idealized method for solving Tx = y. 

By formulating the error Txn — y appropriately, it is possible to reduce the 
search problem to a matrix equation. MoM defines the error in terms of the 
components of the term Txn -y- If { } is a total linearly independent set in 

the range space of T, then MoM enforces the first N projections of Txat — y to 
be zero along { u, }. That is, 

^r(|;7, ei) - y , ufj = 0 forj = 1.2,...,JV (3.109) 

= {y,u,) tori = 1,2,... ,N (3.110) 

i=:l 

The sequence of approximate solutions obtained by solving (3.109—3.110) need 
not be monotonically improving. The method of point matching enforces the 
condition that Txn agrees with y on a set of N points in the range. The key idea 
of MoM is the use of matrix methods to solve the best approximation problem 
of y in the space spanned by { Tci }. 

We prefer to call this as the naive approach to MoM since no attention is 


. ^except for few equality signs this is valid in the form • - • > cn > ^ ^ n +2 > 
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payed to the performance of the algorithm in the presence of errors. We shall 
soon see that the errors in computing Tct or {Te, , Uj) lead to instabilities in 
the solution scheme. 


3.6.3 Practice of Method of Moments 

Introduction 

EM boundary value problems are generally formulated as integral equations 
of the first kind. A Fredholm integral equation of first kind (in one dimension) 
is given by, 


f K{x,y)f{x)dx = g(y), 

J a 


c <y < d 


(3.111) 


Despite the various possibilities in the kernel (such as finite rank, singular, 
higher dimensional etc.), the following result is valid for an integral equation 
of the above form. Since g{y) is defined through a quadrature on f{x), the 
right-hand-side of an integral equation is generally smoother than /(.O. Any 
numerical method which attempts to solve f{x) from (3.111), in essence, is 
doing a differentiation® of the right-hand-side. Therefore (3.111) may not have 
any solution unless the right-hand-side possess a certain degree of smoothness. 
This is especially true if the right-hand-side is known with some error as in 
a measured data. Therefore, the existence of solutions is closely linked to the 
smoothness of the right-hand-side. In the context of EM, where we deal with 
incident fields on the right-hand-side, this problem does not arise. The incident 

®See the following example from Tikhonov [43]: The nth derivative z(t) of the function u{t) 
is the solution of the integral equation, 


in - 1 )! 


2{t) dr = u{t) 


0 
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field is infinitely diflerentiable and often anal 3 d:ic, so that it can be regarded as 
sufficiently smooth. 

Even in cases where the right-hand-side is known to be smooth, and a so- 
lution exists, the direct solution of first kind integral equation is regarded as 
ill-conditioned. Note that ill-conditioning is intrinsic to the smoothing property 
of the quadrature process, and it does not depend on what appears on the right- 
hand-side. A right-hand-side with errors merely brings out the badness of the 
left side. The smoother is the kernel, the worse is the numerical evaluation of 
/(.r). The numerical solution of inverse Laplace transform and the problem of 
numerical differentiation etc. also fall into this category. 

The purpose of the present discussion is to illustrate this problem via. a set of 
examples. We use our results as well as the results from the following sources 
[30, 31, 32, 44]. Despite its widespread use, the stability issues of MoM has 
not received any major attention in EM. The three works mentioned herein 
[30, 31, 32] and the pioneering work of R. Mittra on "Relative Convergence" 
remain the major EM investigations on stability. 

This is in direct contrast with Mathematics where the stability issues of 
first type equations have received serious attention. In fact, the theory of 
regularization and inverse methods evolved based on the difficulties in solving 
first kind integral equations. 

Literature on Stability Issues 

S. W. Lee [30], R^* Mittra [31] and M. Leroy [32] studied the convergence of 
certain discontinuity problems in rectangular and parallel plate waveguides. 
Their investigations were on transverse discontinuities. Note that, this includes 
the iris, diaphragm and bifurcation problems. The ntunerical performance of 
both mode matching and the solution of integral equations via. MoM were 
studied. The first observation was that both the methods lead to the same set 
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of equations and are therefore identical. 

We shall therefore assume that the transverse discontinuity problems are 
solved using integral equation MoM. The Green’s function for transverse dis- 
continuity problems does not have the factor. Therefore, the authors 

have used a waveguide Green’s function of the form, (Note: tr stands for sin or 
cos). 


OO CO 

K{x,y,x',y') = ^ tr tr tr tr (3.112) 

m=l n=l 

The computations were done using a truncated version of the kernel. Let 
P denote® the number of modes summed for the kernel, while N denote the 
number of basis functions used for the approximation of the unknown. 

The authors used an iris in a parallel plate waveguide as the discontinuity. 
This problem was used as a benchmark for the study since it had an exact solu- 
tion via. Wiener Hopf technique. Let s denote the ratio of the plate spacing^°to 
the distance below the iris. 

The authors make the following fundamental observation. For a given P, the 
computed results improves with increasing N and then starts deteriorating. By 
the word deteriorating, we understand that either the results becomes unstable 
and drifts off, or the result converges to the wrong value. Conversely, for a 
given N, increasing the number of modes P takes the solution from instability 
to inaccuracy. Therefore, for a given P, there exists an optimum value of N such 
that the results are not too bad either way. We can represent this symbolically 
as, 

. N 

Given P : inaccurate )■ unstable . 

Optimum N 

It was proved via. rigorous arguments that the optimum value of N for a given 

double summation can be regarded as a single summation via. a suitable ordering. 

^°If the iris does not exist then s = 1. 
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P is such that P/ K -- s. Therefore, the sequence of solutions do not indefinitely 
improve. 

These obsei'vations contain sufficient pointers to all the aspects of stability. 
However, the conclusions are based on an integral equation with a finite rank 
kernel. A kernel is finite rank (also referred to as separable, degenerate) if it is 
of the following form, 

K{x, y, .r', y') = ipkix, y) bfc(x', y'). 

k^l 

Consequently an integral equation with a finite rank kernel takes the following 
form, 


n%K{x\y')f[x',y')dx'dy' = g{x,y) 

k:.l 


(3.113) 


A standard result on integral equations with finite rank kernels is that, they 
have infinitely many zero eigenvalues. That is, there exists a nonzero / such, 
that the homogeneous equation has a solution. To prove this, all we need to do 
is to choose an / such that the P quadratures in (3.113) are zero. For example, 
we can orthonormalize the set { d>k{x\ y ') } in the rectangle [a, b] x [c, d], and then 
choose / from the orthonormalized family for k > P. 

Since the (truncated) kernel for the transverse discontinuity problem is of 
finite rank, we have the result that the numerical solution of transverse dis- 
continuity problem is not unique. Consequently, the ill-conditioning which the 
authors observe does not have a rigorous basis because the solution itself is 
not unique. In fact, there exists a suitable choice of basis such that the MoM 
matrix can always be made singular. 

Except for this unfortimate choice of the problem, the results in these refer- 
ences have significant qualitative validity. 

Our discussion will be based on the computation of scattering from a prmted 
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vertical strip within, a waveguide. This does not involve a separable kernel. 

The reader may refer to Delves and Mohamed [44] or C. T H. Baker [45] for 
further examples. 

Illustration of Instabilities 

Consider the scattering from a printed vertical strip within a waveguide. See 
Fig. 2.1 for the geometrical description of the problem. The integral equation 
is of the form, 

rvi+h 

/ / Ey{Jy) Jy{y',z')dy'dz' = (3.114) 

where 

^y{Jy) = Yu sinui(a-t) coskyy' sinui{a - x) coskyy 

+ Y smvi{a-t) coskyy' sinvi(a - x) coskytj 

(3.114) is the statement that the scattered Ey component cancels the incident 
Ey component on the strip. This integral equation is valid on the strip domain 
Sr, namely 

Vi-k < y < yi+ h 
Sr : -wi < z < +wi 

X = t 

(3.114) is solved using the same basis and testing functions as in Chapter 2. 
Let Mmax (Nmax) denote number of modes summed in the x (y) direction, whereas 
Ni denote the number of basis functions used for the strip current. We have 
measured two cases having different strip widths via. TRL standards. The 
fabrication and the measurements were done carefully so that the measured 
results can be regarded as reliable to check the computed results. 

Fig. 3.9-3.10 show the measured S parameters along with computed results 
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Frequency (GHz) Frequency (GHz) 


Figure 3.9; S parameters of strip scattering, a = 22.86 mm, b = 10.16 mm, 
t = 1.57 mm, fr = 2.5, 2li = 9.0 mm, 2vji = 0 8 mm, yi = 5.08 mm, 9 = 

90°, Mmax = 100, N„iax = 100. Measurements were done using TRL 

standards 

for a specific choice of Mmax, ^max and Ni. Specifically, we have used 5 basis on 
the strip, whereas 100 modes were used in each direction. The results seem to 
suggest that the computations give accurate results and are stable. However, 
this is not quite the case. 

The scattering behavior of the strip is as follows. The strip is a reactive 
element and it appears as a series L-C circuit in shunt. Consequently, there is 
complete reflection (Sn = -1) and no transmission (^21 = 0) at resonance. The 
measured results do not really reach these two values at resonance due to the 
losses in the structure. Fig. 3.10, which represents the wider strip, has a lower 
resonant frequency and has broader bandwidth. 

We present two distinct cases of instability depending on whether the strip 
is offset on the narrow wall or not. An offset strip is one for which the center of 
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Figure 3 . 10 ; S parameters of strip scattering. 2u!i = 1.5 mm. See Fig. 3.9 for other 
details. 

the strip yi is different from 6/2. 

Fig. 3.11 shows the S parameter of an offset strip with respect the number of 
basis functions. The results are computed sufficiently away from the resonance. 
The results are outright divergent for 50 modes^^ and it seems to improve on 
switching to 100 or 200 modes. Even for 100 or 200 modes the results exhibit 
a small yet definite linear growth which indicates instability. We do not have 
measxired data for this case. However, the result seems acceptable at low values 
of Ni before the instability creeps in. 

The computed results for the centered strip of Fig. 3.10 is more interesting. 
Fig. 3.12 and Fig. 3.13 shows the magnitude and phase of Sn with respect to 
the number of basis function. Since the strip is centered on the narrow wall, 
the even indexed basis functions do not contribute. This explains the stair like 
appearance of these curves. We have made these computations right at the 

^^"for N modes" mean "for N modes in each a: and y direction" 
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Figure 3.11: Illustration of instability with respect to the bases, a = 22.86 mm, 
h = 10.16 mm, t = 1.57 mm, = 2.5, 2/; = 8.0 mm, 2wi = 1.5 mm, 
!/, = 6.08 mm, 0 = 90°. 


resonant point so that 5ii = -1. 

It can be seen that the computed result improves, reaches the true value 
and then starts deteriorating. For a given choice of Mmax and Nmax there exist 
and optimum value of Nt such that the results are accurate. An increase in the 
number of basis functions has to be compensated by a commensurate increase 
in the number of modes for maintaining the accuracy. But the relationship 
between modes and basis is not linear unlike the iris problem. 


3.6.4 Conclusions 

The numerical solution of integral equations of the first kind is an ill-conditioned 
problem. Algorithms based on the naive methods presented in Section 3.6 be- 
comes unstable beyond a certain accuracy. 
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Figure 3 . 12 : 5ii magnitude versus number of basis functions at variou.s values of 
number of modes. The number of modes used in each direction (A = A',,,,,,.) is 
indicated on the curve. See Fig. 3.10 for details. 

In EM, we encounter predominantly first kind integral equations with ex- 
tremely well behaved right-hand sides. Consequently, it is necessary to pay 
attention to the quadrature/truncation errors of the lea-hand-side, since there 
are no errors on the right-hand-side. Waveguide Green’s function is an example 
of a kernel, where the quadratures can be done analytically, whereas the series 
can not be summed in closed form. In EM, we also have instances where there 
are quadrature errors and no truncation error. 

The sequence of approximate solutions does not improve indefinitely as in 
theory. Instead, it first improves and then starts deteriorating. This behavior 
seems to depend on a parameter k, such that, the results are best if the right 
value of K- is chosen. For example, the ratio P/N suggested by [30, 31, 321 serves 

this purpose. The optimum choice for P/N was found to be s for the iris problem 
in a parallel plate guide. 
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Figure 3 . 13 ‘ Phase of 5ii versus number of basis functions at various values of 
number of modes. See Fig. 3.12 for details. 

These observations can be viewed in the light of the theory of regularization 
proposed by Tikhonov [43]. Regularization belongs to one school of solving 
ill-conditioned problems. An elementary description of regularization is as 
follows. 

Before we attempt to solve Tx = y, we expect the solution to belong to a 
specific smoothness class. For example, x can be assumed to be square inte- 
grable. Consequently, while generating the approximate solution x^s,, we also 
look whether the sequence ceases to belong to this class. 

If the solutions x^ = li has to be convergent, we expect the sequence 
{ /, } to decrease fast enough. Consequently, instability is caused by an insuf- 
ficient decay rate for the sequence { J, }. For example, if x is square integrable 
and { e, } orthonormal, then 4 = o(l/ Vk). 

If the sequence becomes unstable, it is to be understood that the higher order 
terms are not deceasing sufficiently fast. Therefore, smoothness and stability 
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axe pynrinymrms. Let L be a lineax operator such that \\f--i .v]| is a mcasuro of 
the smoothness of XjV- A. regularized sequence of solutions is goneiutod based 
on both the accuracy l|Tx,v - yll and the smoothness \\Lry\\. A parameter k ties 
up the two such that, we minimize 

llTx.v-yll^ + 

The parameter k is called the regularizer of the problem T.r - y. 

Therefore, a regularized sequence of solutions is a tradeoffbetween accuracy 
and stability. For a given problem, there exists an optimum choice of s. Too 
low a K will lead to very unstable solutions, whereas too high a value will lead 
to inaccurate solutions. 

The parameter P/N is a regularizer for the iris discontinuity. It is possible 
to identify regularizers for other type of problems. For example, If the MoM 
matrix has quadrature errors such that c is the maximum absolute error for 
any MoM matrix entry, then s.V may serve as a regularizer. 

The numerical behavior of MoM solutions leads to certain significant af- 
terthoughts. For an unstable problem, the practitioner of the method arrives 
at some intuitive criterion to stop at a good solution. This is very close to the 
choice of a regularizer and therefore justified. On the other hand, it also makes 
it difficult to distinguish between a formulation which is erroneously modeled 
and one that is unstable. 

Because of these reasons, the computed results presented in Chapter 5 are 
generated for a suitable choice of basis and modes, which were found to give 
acceptable numbers in most of the computational experiments. The number of 
basis functions (A'; and AJ and the number of modes (Mmax and Nmax) ^ased for 
every computed result are indicated in the figure captions. 


derivative can serve as a measure of the smoothness. 



3.7 Summary and Conclusions 


113 


3.7 Summary and Conclusions 

This chapter has discussed various issues of independent interest. Green’s 
functions in inhomogeneously filled waveguides were derived, primarily to con- 
solidate the notations and to understand the structure and properties in depth. 
The derivation made use of the method of normal mode expansion from "Field 
Theory of Guided Waves" by R. E. Collin. The normal modes being LSE and 
LSM, it became necessary to solve the LSE/LSM eigenvalues. The solution of 
eigenvalues were based on certain simple results, which made the salient prop- 
erties of the curves analytically known. This led to the choice of good starting 
values and tight bounds for the roots. 

A Green’s function is a distribution so that it leads to various inconsistent 
results if we attempt to use it as a classical function. This theme was illustrated 
using the standard EaJ^) Green’s function, which is a distribution due to the 
presence of a Dirac delta. It was found necessary to represent this Dirac delta as 
an LSM series so that the E^ field was computable. This section also developed 
an algorithm for implementing the Gram-Schmidt procedure. 

Section 3.6 dealt with an interesting topic, namely, the theory and practice 
of method of moment. MoM is a technique which attempts to solve a first kind 
integral equation directly. Since this problem is fundamentally illconditioned, 
instability creeps in approximately at the same time we tend to get acceptable 
solutions. The practitioners of MoM use various thumb rules to decide when to 
stop while solving a problem. This was interpreted as the choice of a regularizer 
based on regularization theory. 



Chapter 4 


Measurement Procedure 


4.1 Introduction 

In order to verify the computations of the previous two chapters, an experiment 
was designed to measure the equivalent network of the strip excited slot in a 
waveguide. Since impedances cannot be m(?asur(‘d directly, it i.s obtained by 
measuring the S-parameters and then converting this to an (‘(puvalent two- 
port network. The strip excited slot in an inhoraogeneoualy filled wav<‘guicle, is 
treated as a two-port network, with the ports defined at a large distance (few 
wavelengths) from the slot, so that only the propagating mode exists at the 
measurement planes. A vector network analyzer (VNA) is used for measuring 
the S-parameters. For accurate measurements it is necessary to calibrate the 
VNA. The standard waveguide calibration kit consists of short, shifted short, 
load and a standard length waveguide- These are precision fabricated, and 
can be used to establish measurement planes between which the waveguide 
two-port network can be inserted. 
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Since the computations give equivalent network impedances or the 
S-parameters referred to the center of the slot, it is required to shift the mea- 
sured results to the center of the slot for comparison. This can be done theo- 
retically, but in practice, this procedure was found to be unsatisfactory, for the 
reason that the shift is accurate only to the extent of accuracy with which the 
parameters of the waveguide section (inhomogeneously filled) are known. For 
improved accuracy of measurement, Thru-Reflect-Line (TRL) calibration was 
adopted in later measurements. In this procedure, the measurement planes are 
defined at the slot center itself, and the entire inhomogeneously filled waveguide 
section, along with the transition to normal waveguide is de-embedded auto- 
matically, to give the S-parameters of the strip excited slot at the slot center. Of 
course, for this the TRL standards have to be specially designed and fabricated. 
This chapter discusses the various aspects of experiment design, fabrication of 
standards, TRL calibration and some representative measurements. 


4.2 Fabrication of Samples 

The analysis of the strip excited waveguide slot assumes an infinite inhomo- 
geneously filled waveguide. The strip and the slot are assumed to have no 
thickness whatsoever, and the conductor and dielectric losses are regarded as 
zero. The slot is assumed to be embedded in an infinite ground plane. There- 
fore, the test pieces should be fabricated in such a way that these idealizations 
are met at least approximately. 

If the exciting strip and the slot are printed on either side of a low loss 
dielectric substrate, and the waveguide narrow wall is substituted with it, then 
the requirements of theoretical model are nearly met. A ground plane of certain 
size can be attached in flush with the narrow wall of the waveguide to simulate 
the effect of an infinite ground plane. The length of the dielectric substrate can 
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be chosen to be at least 2Xg, so that in the vicinity of the slot, the structure 
resembles a uniform inhomogeneously filled waveguidi'. For low filling factors 
(low er and thin substrate), the length of the substrate us<Ki to simulate an 
inhomogeneously filled waveguide is not very critical. 


WR 90 Wave Guide 



160 

All Dimonstons tn mm 
Not to seals 

Figure 4.1: Details of waveguide test piece 

In order to implement the above scheme of ideas, WR-90 Brass wav<‘guid(‘.s 
of length 160 mm are fabricated. [See Fig. 4.11. The wavogui(l(‘.s are fitted 
with X-Band brass fianges (UCx-39/U). A portion of the narrow wall of size 
66 X 8 mm is machined off and the narrow wall is thinned down to 0.3 mm. The 
waveguide pieces are silver plated after machining so as to keep the waveguide 
wall losses to a minimum. The test samples for measurements are made using 
this machined waveguide pieces. In order to fabricate the strip excited slot 
antenna, we proceed as follows. [See Fig. 4.2]. 



4.2 Fabrication of Samples 


117 



All dimensions m mm. 

Not to scale 

Figure 4.2: Substrate geometry for measurement using standard waveguide calibra- 
tion 


A dielectric substrate is etched with the slot and the exciting line on either 
side. It is cut accurately to the dimensions indicated in the figure, so that the 
substrate fits snugly between the broad walls of the waveguide. The substrate 
is kept pressed from the waveguide interior to the machined narrow wall side, 
and the edge joining the substrate with the narrow wall is soldered around. 
[See Fig. 4.3] . The ends of the substrate for about 20 mm is kept free of copper, 
since this region is not involved while soldering to the narrow wall. 
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Figure 4.3: Details of the test antenna element 
This procedure is carried out carefully, so that, the substrate is mounted 
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sufficiently close to the narrow wall of the waveguide to avoid any air gaps. It 
is necessary to align the slot with the center of the waveguide, so that mea- 
surement reference plane coincides with the slot center. 1 his i.s done with aid 
of a traveling microscope. Since the width of machined off narrow wall is only 
8 mm, it is necessary to leave more space around the slot ends, especially if 
the slot is more than 8 mm long. This is done by removing the remaining 
part of the narrow wall in the neighborhood of the slot, such that there is just 
sufficient narrow wall left for soldering. It is also necessary to remove any 
extra solder remaining on the substrate or the waveguide walls prior to making 
measurements. 

A ground plane fixture is fabricated by attaching two L shaped Aluminum 
plates to a 1 inch ebonite spacer as shown in Fig. 4.4. The waveguide test 
piece can be lowered into the channel between the Aluminum plates and the 
spacer, so that the machined narrow wall plane becomes flush with the ground 
planes. Copper foil tape with a conducting acrylic adhesive is used to join the 
waveguide wall with the ground plane. 

This method of fabrication, though unsuitable for a practical array, is found 
to be good enough for laboratory evaluation. If the calibration used is the 
standard waveguide cal-set, setting the reference plane at the flange planes, 
it is necessary to taper the substrate ends so as to minimize the reflections 
from the dielectric step discontinuity. The reflections from the dielectric step 
can be of the same order of magnitude as the reflections from the strip excited 
slot, especially if the latter is operated at lower conductance ranges. This is 
especially the case in an array of 30 or more elements, where each slot appears 
as a small discontinuity. If TRL calibration is used for measurements, it is not 
necessary to taper the substrate ends since the effect of the step discontinuity 
is automatically de-embedded in the calibration procedure. 

The measurements are done on a HP8510C Vector Network Analyzer. The 
measurements are done on two types of substrates, namely, Taconic TLY-5- 




Figure 4.4: Ground plane fixture 


0310-cl/cl (cr = 2.18, substrate thickness= 1/32 inch) and ET/duroid 5550 
(fr = 2.5, substrate thickness= 0.062 inch). 

Fig. 4.5 illustrates the S-parameters measured using standard waveguide 
calibration. It can be seen that S 21 phase does not pass through zero as required 
at resonance. This is due to the errors in shifting the reference plane from the 
flange to the slot center. The phase error is caused mainly by the errors in 
the distance involved as well as in the propagation constant. A distance error 
of 0.1 mm corresponds to approximately 1.8° in ^21 at 10 GHz. The phase 
of Sii also contains shifting errors of the same order as in 52i. However, 5ii 
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3 degree phase error in 52i due to the errors in shifting the reference plane from the 
flange plane to the slot center. No ground plane is used. 0 — 50°, 2// = 8.0 mm, 
2xui = 0.8 mm, yi = 5.08 mm, z; = 0 mm, 2^* = 8.5 mm, 2a», = 0.8 mm, 
er = 2.20, t = 0.8 mm 


phase error is not perceptible since it is in a scale of 180 degrees. Note that, 
while shifting the plane of reference, we encounter three regions starting from 
the flange plane to the slot center. There is a distance of 27 mm where the 
waveguide is unfilled, followed by a tapered region of 20 mm, and then a length 
of 33 mm where the substrate has full thickness. Since the propagation constant 
of the tapered region is not known exactly, either the value in the unfilled region 
or the value in the filled region has to be used in this section. This causes some 
error in the shifted S-parameters. Since, the reflections from the dielectric 
step are part of the measurement, the measured S-parameters will be slightly 
erroneous. The tapering of the substrate makes the substrate ends flexible and 
the chances of air gaps are increased. 
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4.3 Calibration Schemes for HP8510C VNA 

A nstwork mcasuremont using VNA is reliable only if the errors in the mea- 
surement circuitry are identified and corrected for. The errors encoxmtered 
in a typical measurement circuit are of varied t 5 pes and affect the measure- 
ment to varying degrees. However, some class of errors are repeatable in a 
measurement session. The errors are modeled by introducing an error model 
which involves a certain number of error terms. In a calibration procedure, a 
set of devices (standards) with known characteristics are connected to the mea- 
surement ports of the VNA. The errors parameters are computed based on the 
measurements using the standards. The estimation of the error terms involve 
the solution of a set of complex, non-linear equations. A formal treatment of 
the theory of calibration schemes can be found in [46]. 

In the conventional waveguide calibration scheme, a set of three discrete 
impedance standards (namely, short, shifted short and load), and a transmis- 
sion standard which is a standard length waveguide are used. These standards 
are precision fabricated and should have properties known to high degree of 
accuracy. This technique can be used to establish a reference plane at the 
center of an inhomogeneously filled waveguide. However, this requires high 
accuracy on the test pieces, in addition to having to fabricate and characterize 
a number of standards. The next two subsections illustrate the technique of 
establishing the reference plane at the center of slot using conventional as well 
as the Thru-Reflect-Line standards. 


4.3.1 Standard Waveguide Calibration 

Fig. 4.6 illustrates the necessary standards required to establish the reference 
plane at the center of the test waveguide. SHORT is a half length waveguide 
fitted with a dielectric substrate such that it is identical with one s ymm etric half 
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Figure 4.6: Standard waveguide calibration modified to establish the measun.'metJt 
plane at the center of an inhomogeneously filled waveguide 


of the test waveguide piece. SHIFTED SHORT is identical with the SHORT, except 
that the former is slightly longer than the SHORT waveguide. MATCHED LOAD 
is implemented by connecting a matched termination at the end of the SHORT 
waveguide. THRU is a definite length of inhomogeneously filled waveguide 
which is introduced between the two half length waveguides. Note that, SHORT, 
SHIFTED SHORT and LOAD constitute three single port measurements (totaling 
to six), while THRU constitute a full two port measurement (totaling to four). It 
IS necessary to have two half length waveguides since the THRU requires a half 
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length waveguide on either side. The shorting plate can be from a standard 
waveguide calibration set. However, the matched load is to be specifically 
fabricated such that it is a matched termination for an inhomogeneously filled 
waveguide. It is necessary to fabricate four standards namely SHORT (2 Nos. 
), SHIFTED SHORT (1 No. ), MATCHED LOAD (1 No. ) and THRU (1 No. ). It is 
also necessary to characterize the length of the SHIFTED SHORT, THRU and the 
reflection coefficient of the LOAD very accurately to implement this calibration 
scheme. 


4.3.2 The TRL Calibration for Waveguide 

TRL calibration makes use of the properties of a short piece of transmission 
line and does not require a set of three precision impedance standards [47]. 
It assumes an 8 term error model for the measurement circuit. To do TRL 
calibration, it is necessary to carry out the following three steps, which amounts 
a total of ten measurements. 

• THRU: connection of port 1 and port 2, directly or with a short length of 
transmission line. If directly connected, it is called zero length THRU 

• REFLECT; connect identical one-port high reflection coefficient devices to each 
port 

• LINE: insert a short length of transmission line between port 1 and port 2. 

The lengths of THRU and LINE should be tmequal. The difference in line length 
between LINE and THRU should preferably Xg/A at the center of the frequency 
range of calibration. For X-band waveguide Xg/A corresponds to about 10 mm 
at 10 GHz. 

Since the TRL error model has only 8 terms and the number of measurements 
total to ten, it is possible to compute the two constants defining the REFLECT and 
the LINE standards. Therefore, the complex reflection coefficient of the REFLECT 
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and the propagation constant of the LUSnE need not be specified whil(‘ specifying 
the cahbration standards. It is only necessary to specify the characteristic 
impedance of the LINE, and this becomes the measurement reference. In the 
case of waveguide TEL, all the measurements are carried out with respect to a 
- normalized characteristic impedance of unity, and the characteristic impedance 
entry is 1 Ohm. 

4.3.3 TRL Calibration Standards 



Figure 4.7: TRL calibration standards 


Fig. 4.7 indicates the necessary calibration standards required for setting the 
measiarement plane at the center of the slot [See also Fig. 4.3|. Two machined 
160 mm waveguide pieces are used to make the THRU and the LINE. In order 
to fabricate the THRU, a dielectric substrate of dimensions 106 x 10.16 mm is 
fitted symmetrically inside the waveguide where the narrow wail is machined 
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off. The LINE is identical to the THRU except for the increase in substrate length 
by 10 mm. The substrate can be attached to the waveguide narrow wall by 
soldering from outside as in the case of the test pieces. Note that the facing 
side of the substrate should be copper clad. 

WR 90 Wave Guide 


UG-39/U Flange 

/ 



All Dimensions in mm 
Not to scale 

Figure 4.8: Half length waveguide with flanges 

A half length waveguide of 80 mm is fabricated for use as a REFLECT. [See 
Fig. 4.81 It is fitted with Brass X-Band flanges UG-39/U and is silver plated. A 
dielectric substrate of 53 mm is fitted inside this guide, such that the su 
end is in flush with the flange plane. This waveguide along with a shorting 
plate act as the REFLECT. Since the REFLECT simulates one symmetric half of 
the THRU, it can be used to establish a short at the measurement reference 

Two 5 mm waveguide extenaiona are also fabricated as indicated in Fig. 4.9. 
Tire choice of 5 mm correaponda approximately to A,* at 10 GHx. Therefore, 
attaching the 5 mm extenaiona to either sides of the line wavegmde correspon a 
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Ail Dimensions in mm 
Not to scale 

Figure 4.9; Details of the 5 mm waveguide extonsicm 


to a Ap/4 extension at the measurement plane. 


Fig. 4.10 indicates the connection schematics for TRL standards. The short- 
ing plate for the REFLECT can be from the standard wav(‘guidt! calibration set. 
Note that, this scheme corresponds to a zero length THRU and a \/4 LINE im- 
plicitly. Since the waveguides used for THRU and test pieces are identical, the 
measurement reference at the THRU center corresponds to the slot center as 
well. The waveguides fabricated for making the test pieces are found to have 
an error of 160 ± 0.3 mm. Therefore, it is necessary to introduce small shifts in 
reference plane in TRL measurements as well. However, this is of the order of 
1—2 degrees unHke the measurements using standard waveguide calibration. 
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Short Short Coax to waveguide adapters 



Figure 4.10: Connection schematic for TRL standards 


4.4 Few Representative Measurements 

This section contains few representative measurements which are used to illus- 
trate certain aspects of the measurements. The examples also serve to illustrate 

the behavior of the strip excited slot. 

Fig. 4.11 presents the S-parametersofa typical dot measured using the TEL 

calibration. The solid lines indicate measurement with the ground plane, and 
the dotted lines correspond to the measurement without the ground plane. It 
can be observed that S 21 phase passes through zero at resonance, 
measurements using the standard waveguide caUbration. The effect of the 
ground plane is a small shift in resonant frequency, approximately of the order 

of 0.5%. . 

The S-paxameters can be converted to normalized Z-parameters using 
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Frequency (GHz) Frequency (GHz) 


Fignre 4.11; S-parameters measured by TRL standard with and without {'round 
plane. 6 = 70°, 2/; = 9.4 mm, 2wi — 0.8 mm, //; = 5.08 mm, r.t () mm, 

2ls = 8.0 mm, 2ws = 0.8 mm, = 2.2, t = 0.8 mm 

transformation rule Z = {I + 5)(/ - SyK The T-equivalent network pves a 
better insight into the behavior of the discontinuity than the S- Parameters. 
With respect to Fig. 4.12, the impedances and :2 form the Heriti.s arms, while 
the impedance zt forms the shunt arm. Note that, all the impedances are 
normalized with respect to the characteristic impedance (wave impedance) of 
the waveguide for the propagating mode. 

Fig. 4.13 gives a plot of the network parameters of the strip excited slot. 
The series arm impedances are seen to be small compared with the shunt 
arm impedance. Therefore, the strip excited slot is predominantly a shunt 
discontinuity. Note that, the real part of the shunt impedance is a constant 
while the imaginary part goes from negative to positive values. This is an 
example of a series resonance. It can therefore be stated that the strip excited 
slot is a shunt admittance, which goes through a series resonance. 






Impedance (normalized) Impedance (normalized) 
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Figure 4.12: T-equivalent network 



Frequency ( GHz) Frequency ( GHz) 

Figure 4.13: T-network parameters of the strip excited slot with ground plane. See 
Pig. 4.11 for fiirther details. 
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Frequency (GHz) 

Figure 4.14: Multiple resonance exhibited by the atrip excited slot. 0 72‘’, 2// ~ 

6.95 mm, 2wi = 0.6 mm, yi = 5.08 mm, C/ = 0 mm, 2/, •- 9.0 mm, 2ti\ = 
1.5 mm, ys = 5.08 mm, €r = 2.5, t = 1.6 mm 
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Frequency (GHz) Frequency (GHz) 


Fifrurc 4,15; Normalized shunt admittance of the strip excited slot, illustratmg the 
effect of clo.sing the slot. 9 = 40°, 2li = 7.37 mm, 2wi = 0.8 mm, yi = 5.08 mm, 

Zt = 0 mm, 2/., = 8.0 mm, 2ws = 0.8 mm, = 2.5, t = 1.6 mm 

The strip excited slot can exhibit multiple resonances for certain choice of 
parameters. This is illustrated in the following Fig. 4.14. There exists 3 res- 
onances starting with the series resonance at 11.5 GHz, followed by a shunt 
resonance at 12.25 GHz, and a series resonance at 13 GHz. All resonances 
except the first are found to radiate very little. The bandwidth of the struc- 
ture is considerably improved [compare Fig. 4.14 with Fig. 4.11], if the slot is 
made longer and made to resonate at the higher end of the frequency band. 
This amounts to the statement that the slot should be operated closer to its 
unloaded resonance for broader bandwidth. Both the conductance and the sus- 
ceptance exhibit asymmetry around the first resonance. It can also be observed 
that the conductance peak does not coincide with the zero of the susceptance. 

Measurements are carried out to study the behavior of the scattering from 
the strip in the absence of the slot. [See Fig. 4.15]. The slot can be covered 
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by a conducting tape while making these measurements. It is observed that 
the strip resonance is narrow and it occurs at another frequ(;ncy. Note that, 
the strip admittance is not imaginary and has in fact a significant value for 
the conductance, especially around the resonance. This conductance accounts 
for the losses in the system and it is significant around the strip resonance. 
The magnitude of the induced current on the strip is quite large if the strip is 
resonant, and this leads to conductor losses. Since losses are directly related to 
the magnitude of the strip current, it is not significant if the slot is radiating. 


4.5 Summary and Conclusions 

This chapter dealt with the topic of measurements. To validatti the computa- 
tions, it is found that the experimental model should meet the assumptions of 
the analysis at least partially. A method of fabrication of the strip excited slot 
was evolved which met this criterion. The technique was to attach a printed 
dielectric substrate to the machined off narrow wall of a wuvi'j’uido. A sim- 
ple ground plane fixture was fabricated which could simulate the effect of an 
infinite ground plane. 

Measurements using the standard waveguide calibration were found to be 
inaccurate due to the errors in shifting the mea.sunmu-nt plane. The substrate 
used for making the test pieces were tapered to reduce the reflections from 
the dielectric step discontinuity. A calibration scheme such as TRL is found 
to be excellent in setting up the measurement plane at the center of the slot. 
The design and fabrication of the necessary TRL calibration accessories were 
discussed. 

A few representative measurements were presented which throw light on 
the various facets of the strip excited slot. The set of measurements indicate 
that TRL calibration gives accurate results. The ground plane is found to shift 
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the resonant frequency by a small amount. The strip excited slot is found to be 
a shunt discontinuity which exhibits a radiating series resonance. It can also 
exhibit multiple resonances in which case the admittance curves are skewed 
and asymmetric. The strip excited slot has a larger bandwidth if the strip 
loading on the slot is kept to a minimum. The conductor losses were found to 
be significant if the strip was made resonant in the absence of the slot. 



Chapter 5 


Characteristics of Strip Excited 
Slot 


5.1 Introduction 

The strip excited narrow wall waveguide slot was conceivcti mainly for use as 
an array element. Details of the analysis procedure and experimental tech- 
niques have been discussed in the previous three chapters. The current/field 
distributions on the strip and the slot were obtained via. solution of a set of 
coupled integral equations, using MoM procedure. Experiments were designed 
and carried out to characterize the strip excited slot experimentally. How- 
ever, experiments are expensive and time consuming to perform, for complete 
characterization. Therefore, few selected measurements were carried out to 
validate the theoretical computations and also to ascertain the errors involved 
in the computations. Knowledge of the accuracy and reliability of the theoreti- 
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cally computc-d characteristics enable one to evolve a mechanism of correcting 
these errors, so that, the computed results are useful in the design of high 
performance antenna arrays. 


Precise control of the excitation coefficients is the key to a high performance 
antenna array. Practical realizations of this antenna is directly dependent on 
the ease with which these excitation coefficients can be realized and maintained 
over the bandwidth. Therefore, it is necessary to choose one or two of the 
parameters of the strip excited slot as variables of control and keep the rest of 
the parameters constant. Of course, appropriate choice of constant parameters 
is also important. Strip excited slot has several parameters, such as, slot length 
and width, atrip length and width, strip inclination, strip offset, dielectric 
constant and thickness of the substrate, waveguide dimensions etc. Certain 
parameters such as the vertical offset do not significantly affect the behavior 
of the slot. On the other hand, there are others with respect to which the 
impedance characteristics is exceptionally sensitive. Therefore, it is necessary 
to identify the behavior of the strip excited slot with respect to each of the 
parameter, and single out those which are relevant in the context of array 
design. This calls for an enormous amount of data, representing the behavior 
of the slot with respect to various parameters. In order to generate this data, 
the analysis program is preferred to an extensive series of measurements. It , 
then becomes essential to validate the program with select measurements. 


A 

The chapter therefore begins with a comparison of theoretical results with 
the measured results. Once the validity and limitations of the theoretical model | 
are understood, it is possible to generate the parametric data through the 
program. The data sets generated are governed primarily by the requirements 
of an array design algorithm. 
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5.2 Comparing Theory with Experiment 



Frequency (GHz) Frequency (GHz) 


Figure 5.1: S-parameters of the strip excited slot. Mcfa-surements W(‘rc don<* using 
TRL standards with ground plane, waveguide: (i •- 22.86 mm, h 10. 16 mm; 
strip: 9.4X0.8 mm, 80°; slot: 8.0 x 0.8 mm, offsets: /// -- 5.08 mm, :/ ■ 0 mm; 
substrate: = 2.2, t = 0.8 mm; modes, basis: M,ni,x ■" ‘^'nuu - 100, Xi .V., — 

9 

Various sources of errors exist in the process of nu>:is\iri«in(‘nt us well as com- 
putation. It is therefore necessary to have an estimate of the agreement to 
be expected between theoretical and measured results before arriving at any 
definite conclusion. In practice, various parameters of the test antenna are 
known within a tolerance. For example, the substrate thickness as well as the 
permittivity are known within 2% of their actual values. The accuracy of the 
printed circuit fabrication in the laboratory is of the order of ±0.05 mm. This 
affects the linear dimensions and the inclination of the strip. The computed 
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Fifiure 5.2: S-parameters of the strip excited slot. Measurements were done us- 
ing TRL standards with ground plane, a = 22.86 mm, b = 10 16 mm; 
strip: 9.4 y 0,8 mm, 70°; slot: 8.0X0.8 mm; offsets: yi = 5.08 mm, Zi = 0 mm; sub- 
strate: <r ~ 2.2, / = 0.8 mm; modes, basis: M„iax = Mnai = 100, Ni = = 9 


results depend on the assumptions made in the formulation. We have, for 
example, neglected the transverse current flow on the strip and assumed the 
longitudinal current to have a symmetric transverse variation. The extent to 
which these approximations affect the computed results will come to light in 
the comparison. 

Fig. 5.1 to 5.6 will be used for comparing theory with measmement. First 
three figures give measured and computed Fn and F 21 , and the other three 
present normalized admittance. All the parameters are shown as a function 
of frequency. The first three cases were measured using TRL standards with 
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Frequency (GHz) Frequency (GHz) 

Figure 5.3: S-parameters of the strip excited slot. Measurements were dont* using 
TRL standards with ground plane, strip: 10.0x0.8 mm, 70'^; slot; B.O "OH mm; 
offsets: yi = 5.08 mm, Z; = 0 mm; substrate: = 2.2, t - 0.8 nuts; niuties. bst- 

Sis: Mmax — = 100) M = 9 


ground plane, whereas the latter three were measured using standtird waveg- 
uide calibration without ground plane. Consequently, for the last three cases 
we choose to present the normalized admittance {The admittance is normal' 
ized with respect to the LSEiq wave admittance), rather than the S-parameters, 
since the latter contains appreciable errors due to shifting of the measurement 
plane. The susceptance curves of Fig. 5.4 to 5.6 have been shifted vertically to 
account for the errors in measurement plane shifting. 

Referring to Fig. 5.1 it can be seen that the agreement between theory and 
experiment is excellent. This corresponds to the case when the strip angle 9 is 
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Fi^run* r).4; Normalized shunt admittance of the strip excited slot Measure- 
nufiits W(;rt< doiuf u.sing standard wave guide calibration without ground plane, 
strip: H,0 V 0.8 mm. oO'’; slot: 8.5X0.8 mm; offsets: yi = 5.08 mm, c/ = 0 mm; sub- 
strate: ( r ^’2, i ~ 0.8 mm; modes, basis; Mmax = -Tmai = — 9 


80® which is nearly vertical. Next two cases, namely Fig. 5.2 to 5.3 corresp 
to a strip angle of 70», but different strip lengths. It can be seen that the 
computed resonant frequency is lower than the measured by about 1% in these 
2 oases. Referring to Fig. 6.4 to 5.6, with corresponding stnp mclmations o , 

30°and20°,theshiftinresonantfrequencyisseentobemcreasing. ote a, 

for these three cases, the computed values of resonant conduct^ce as weU as 
the bandwidth are considerably off compared with the measured resu s. 

The following general observations can be made based on the 6 figures 
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-o 



Figure 5.5: Normalized shunt admittance of the .strip excited .slot. Mea.surt!- 
ments were done using standard wave guide calibration without ground {)lane. 
strip: 10.0X0.8 mm, 30°; slot: 9.0X0.8 mm; offs(!t.s: /// ■■ 5.08 mm. :/ () mm; .sub- 
strate: e, = 2.2, t = 0.8 mm; modes, basis: M„„„ 100, .\) . .V, - 9 


computed resonant frequency is always lower than the measured; the error in 
resonant frequency appears directly related to the strip inclintition and can be 
seen increasing progressively through the 6 cases; Fig, 5.7 shows the error io 
computed resonant frequency as a function of strip angle 0. The error is seen 
to be about 4% between 0<9< 60°, and falls off to a small value at = 80°. 

In order to identify the source of this error, a series of experiments were 
conducted to eliminate various possible factors. The substrate RT/duroid 5550 
has the following material specifications provided by the manufacturer; f. = 
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Figure 5,6; Normalized shunt admittance of the strip excited slot. Measure- 
jnents were done using standard wave guide calibration without ground plane. 

(I - 22,86 mm, b ••= 10.16 mm; strip: 10.0X0.8 mm, 20°; slot: 9 0x0.8 mm; 
offsets: iji ™ ,5.08 mm, Zt = 0 mm; substrate: €r = 2.2, t = 0 8 mm; modes, ba- 
sis: ~ A'rdflx “ 100, A/ ~ Ng ~ 9 

2.5 ± 0.04, substrate thickness= 0.062 ± 0.002 inches. The accuracy to which 
the linear dimensions could be controlled are about ±0.05 mm in the laboratory. 
This, incidentally, affects the dimensions and the inclination of both the strip 
and the slot. There also exist small alignment errors between the strip and the 
slot while making a double sided exposure of the substrate. This appears as a 
small strip offset. The program was run with slightly different parameters to 
check for the effect of parametric uncertainty. The computed results showed a 
maximum error of about 1 % in the resonant frequency due the uncertainty of 






strip angle 0 (degrees) 


Figure 5.7: Percentage difference between measured and c(jmputed resonant fre- 
quency {fmeas " f theory) ! fmeas. VerSUS Strip anglo. 

parameters. The error between theory and mciasureinent below d 50'^ is more 
than this, and hence it can be concluded that the uncertainty in the material 
parameters and fabrication tolerances are not the only cause of error. 

There exists a 0.3 mm step in the region where the substrate is joined with 
the waveguide narrow wall. If this affects the slot as an end loading, then 
the measured resonant frequency could be different from the cniuputcMl values. 
Measurements were repeated after thinning down the machined off narrow wall 
completely along with a ground plane. Mea.suroment.s were also carried out on 
very narrow slots in order to reduce the effect the longitudinal component of 
the slot E-field. However, these additional set of oxperimentH could not account 
for the order of magnitude errors appearing in Fig. 5.4 to 5.6. Therefore, it 
was concluded that the errors arise from the formulation rather than from the 
measurement. 

It is possible to solve for the scattering from the inclined strip in the absence 
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of the slot. If we compare this result with the measurements on the strip 
excited slot in the absence of the slot (slot is closed with a conducting tape), 
then it is po.ssible to chock whether the problem arises from the modeling of 
the interaction of the strip with itself This experiment, interestingly lead to 
identical obseiwations as in the case when the slot was present. Thus there 
exist errors in the inclined strip scattering formulation itself 

We now accept the premise that the errors arise from the formulation and 
try to identify the reasons. An inclined strip inside a waveguide carries both 
longitudinal and transverse currents. Both the current components has vari- 
ations in longitudinal and transverse direction. If the length to width ratio of 
the strip is largo, the following assumptions are generally made. 

• The transverse variation of the longitudinal current is assumed to be 
either a constant or an even edge current distribution. In the present 
case*, an even edge current distribution is assumed. 

• The trausvcr.M' current is assumed to be zero. This is justified so long as 
its magnitude is much i^maller than the longitudinal current. 

Since these two a{)proximations were made in the formulation, it is believed 
that the errors arise from these approximations. The incident E-field in an 
inhomogeneously filled waveguide has only the y component. This can be re- 
solved into two components, namely along the strip and transverse to the strip. 
It is intuitively clear that the longitudinal E-field induces the longitudinal strip 
current, while the transverse E-field induces the transverse current. As the 
strip angle (I is reduced, the longitudinal E-field decreases as sin^, whereas 
the transverse E-field increases as cos^. Therefore, the longitudinal current 
decreases and transverse current increases with strip inclination. In fact, at 
small strip angles, the values of the transverse and longitudinal currents can 
become comparable. 
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The above reasoning is partly influenced by the following observation made 
by McNamara [20], and Espino [19], regarding the longitudinal ofiset broad wall 
slot— A longitudinal slot is excited primarily by the waveguide wail currents 
which flow transverse to it. This transverse current flow is odd about the 
broad wall center line. Therefore, increasing the slot offset from the center line 
increases the slot excitation. However, in practice (in an array environment) 
the excitation required for a longitudinal slot is quite small which results in 
small offset slots. Such slots often have offsets of the order of the slot width and 
are called weakly excited . The weakly excited slot is found to bc' different from 
its strongly excited counterpart in various aspects. Irrespective of the offset, 
every longitudinal slot is also excited by a longitudinal wall current causing 
a longitudinal E-field in the slot. This current is even about the centerline 
and therefore its effect on the slot scattering is negligible at larg(? offsets. In 
a weakly excited slot, the magnitude of the longitudinal E-fiold is found to be 
comparable to that of the transverse E-field. This is true even in the case when 
the slot length to width ratio is large. It is also obseived that tlu.* transverse 
variation of the transverse E-field is asymmetric about the slot ccmter line. 
Thus the classical assumptions for the large offset slots are no longer valid in 
the case of a weakly excited slot — An analysis of the weakly t‘.xcited .slot (with 
the assumptions made in the case of a large offset slot) tends to give results 
which are off from the measured results. 


The modeling of the strip excited slot, taking into account both the current 
distributions and allowing for variation in both the directions will be exception- 
ally difficult. The present analysis however predicts the trend in the behavior 
correctly. Consequently, it can be used to conduct a systematic preliminary set 
of experiments. It is perhaps worthwhile to characterize this error and correct 
for it, rather than attempting to solve the problem with complete rigor. 
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5.3 Behavior of the Strip Excited Slot 

The next 8 .sections describe various aspects of the strip excited slot. Since, 
the discus.sion is based on graphs, certain remarks are made to begin with. 
Standard waveguide calibration was used in the initial set of measurements, 
whereas TRL calibration was used for subsequent measurements. The mea- 
sured susceptance curves using standard waveguide calibration are shifted 
appropriately to account for the errors in shifting the measurement plane. It 
is also indicated, whether a ground plane was used while making the measure- 
ments. 

b’or tht‘ computed results, M,nar (Awar) is the number of modes summed in the 
.r (//) direction; (.VJ denote the number of basis functions used on the strip 
(slot). The measured and computed results are presented together (if possible) 
to indicate the trends. This however should not be taken as an indication of the 
agreement between theory and experiment. This aspect has been discussed in 
the previous section. 


5.3.1 S-parameters, Equivalent Circuit 

The strip excited slot is a resonant structure and is relatively narrow band. 
Off resonance, the structure scatters very little, with the result that the S- 
paramefor.s resemble that of a waveguide section; Sn « 0 and Sai 1, indicat- 
ing very little reflection and almost complete transmission. However, around 
the resonance, there is appreciable scattering, and both and Sai attain non- 
trivial values. The resonant point is defined as the frequency at which the 
S-parameter 3 are real at center of the slot. The phase of Sn and S„ 180 
and 0“ respectively at resonance. The peak values of |Sni and |Sn|. however 
occurs at slightly lower frequency. In fact, |Snl and ISnl are not exact y sym- 
metrical about their peak values and are slightly skewed. The strip exci e so 
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Fi^re 5 8; Normalized radiated power ( 1 - |.Siip - P ) vh. frequency at variouK 
strip angles, waveguide: a — 22.86 mm, b = 10.16 mm; strip: 10.0 0.8 mm; 
slot: 9. 0x0. 8 mm, offsets: yi = 5.08 mm, c; = 0 mm; substratt': t, 2.2, 
t 0.8 mm; modes, basis: 50, j\i * <5, 

is also a syrtimetric passive structure so that .sVi -- ontl .sV,! -Si • 

The value of 1 - |S'ii|^ - |52if'^ gives the ratio of the powttr Itmt in.Hide the two- 
port to the incident power. If the losses in the structure are .small compared 
with the radiated power, then this represents the norinali/cd radiated power. 
[See Fig. 5.8]. If the structure is loss less (as in the case of scattering from the 
strip in the absence of the slot), then the identity |5ul”'^ -l- l-S'aii'** - I can be used 
to check whether the formulation is erroneous. However, this is a necessary 
condition and is not sufficient to comment on anything else. Therefore, if a 
scattering formulation does not satisfy this identity, we can conclude that the 
results are incorrect. 

The T-equivalent network gives better physical insight into the behavior 
than S-parameters. The T-network is drawn with respect to the propagating 
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LSI-ho mode in an inhomogeneously filled waveguide. Referring to Fig. 4.12 
and iMg. AAA iSrr (Juipter 4), it can be seen that the series arm impedances 
are about I/IO that of the shunt arm. The variation of the shunt impedance 
with rc>spc>ct to frc‘qucncy, suggests that it is a series resonance. Thus the strip 
excited .slot i.s predominantly a shunt admittance which goes through a series 
resonance. Off i(*sonance, the admittance is very small indicating the absence 
of any .scattering. 



Figure fi.t): Slot voltage magnitude vs. frequency at various strip angles. See Fig. 5.8 
for other details. 


If a discontinuity is a shunt admittance, it can be shown that the maximum 
power transfer takes place when the normalized admittance is 2 + jO. Thus 
a normalized shunt admittance of 2 corresponds to maximum radiation firom 
the slot. It is also obvious that this corresponds to the case when half the 
incident pow'er gels radiated, while the other half gets distributed as reflected 
and transmitted power. It can be concluded that, for a given incident field, 
the rnaximum slot voltage occurs at this value of admittance. With respect to 




Frequency ( GH2 ) 


Figure 5.10: Slot voltage phase vs. frequency at variuu.s ,s!np anjjh-s So* Kig fj.H for 
other details. 

Fig. 5.9 to 5.10, the magnitude of the slot voltagi* is seen to he inaxintum, while 
the phase of the slot voltage is 90'* at re.stjnance. Tin* slot voitagi? is specifujd 
with the assumption that the incident E,, component is ! » J ‘ at the centtjr of 
the slot atx = t. 

The strip excited slot exhibits multiple rf.'.nn.tnci*. duo to tho various inter- 
actions between the strip and the slot. Fig. 4.14 indicates a .seciuencc* of three 
resonances series-shunt-series starting from the lower end of t he band. Exper- 
iments indicate that only the first series resonance t ttrro -ptuid.- to appreciable 
radiation from the slot. This suggests the following approximafo equivalent 
circuit for the shunt admittance. With respect to Fig. 5. 1 1, one can identify two 
series R-L-C circuits connected in parallel. The first R-L-C circuit corresponds 
to the slot, whereas the second corresponds to the strip. The resistance /?(„„ 
is introduced to account for the losses in the strip. Each R-L-C circuit has its 
own series resonance and the combined admittance behavior is identical to that 
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f), 1 1 : Approximate equivalent circuit of the strip excited slot 

shown in Fig. 4.14. This equivalent circuit is approximate since the two circuits 
interact to some extent. 

In th(‘ absence of the slot, the strip scatters on its own without any radia- 
tion. The strip appears as a shunt susceptance which goes through a series 
resonance. Th(‘ measurements, however indicate that the strip admittance has 
an appreciable real part. This is due to the losses in the structure. In the 
absencx* of the slot, the resonant strip current reaches much larger values than 
what is induced in the presence of the slot at the combined resonant frequency, 
which is away from the strip resonance. This leads to increased conductor 
losses in the strip. 


5.3.2 Strip Angle 

The incident wave does not couple directly to the slot and cause radiation. The 
inclined strip acts as the intermediary between the vertical slot and the incident 
propagating mode. Specifically, this involves two coupling stages, namely the 
incident wave to the strip, followed by the strip to slot coupling. The strip 
scattering is maximum when it is kept vertical and it is zero when the strip 



150 


Characteristics of Strip Excited Slot 



Figure 5.12: Normalized shunt conductance vs. frequcrncy at varioti.s .strip angles. 

Measurements were done using standard waveguide calibratiutt and nn gnaind plane. 

See Fig. 5.8 for other details. 

is made horizontal. However, the strip to slot coupling is nia.xiinuni when the 
two are orthogonal (i. e. when the strip is hotizonhil) and whan t.h(‘ two 
are made parallel. Therefore, the strip excited slot cease.s to radiate at both 
9 = 0 and 9 = 90. The reduction in coupling around (I 90 is extremely 
fast, and it happens in a very small angular sector. The analysis pnigram has 
a restriction and cannot be used to generate data for nearly vertical strips. 
Therefore, no theoretical data is provided to substantiate this claim. However 
there are measured data in this region and it confirms this observation. Since 
this data does not contain any other information it is not present here. The 
S-parameters of the strip excited slot is identical to that of the strip inside the 
waveguide when the strip is made vertical 

Fig. 5.12 to 5.14 illustrate the effect of strip angle at lower angles, whereas 
Fig. 5.15 to 5.17 illustrate the effect of the strip inclination at higher angles. The 
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Figure 5.15- Normalized shunt conductance vs. frequency at various strip angles. 
Measurements were done using TRL calibration with ground plane, waveguide, tt 
22.86 mm, b = 10.16 mm; strip: 9.4x0. 8 mm; slot; H.OsO.H mm, offsets: j// '■ 

5.08 mm, zt — 0 mm; substrate: f,. = 2.2. t •= 0.8 mm; mode.s, basis. '■ 

Nmax = 100, Ni = 5, = 3 

following observations can be made regarding tlu* effect of ilit* .strip angle. A 
very wide range of conductances can be realized within O'* to HO". Tht? resonant 
frequency increases with increase in the strip inclination. It can bi* .seen that 
the change in resonant frequency with respect to strip indintition i.s faster in 
the higher angle ranges. 

Altering the strip angle from 9 to -0 changes the sign of the slot voltage. 
This is analogous to flipping the offset in a longitudinal slot or altering the tilt 
in an inclined narrow wall slot. It is required the flip the strip angle to cancel 
the propagation phase shift of 180° between two ac^acent slots, so that the slot 
fields are in the same direction. 
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5.3.3 Strip Length 



Figure 5.18: Normalized shunt conductance v.s. fn;qiu-ncy :.t various .strip lengths. 
Measurements were done using TRL calibration with gmuntl plam*. wavi-guule. n . 
22.86 mm, b = 10.16 mm; strip: 2ivt - 0,8 mm. 70", .slut; HO • OH nmi; ofO 
sets. yi 5.08 mm, 0 mm; substrate: t f, 2.2, t ■ OH nun; nujtles, ha* 
sis. Mfnax ~ ^max ~ 100, Ni — 5, Ng “ 3 


As shown in Fig. 5.18 to 5. 19 the strip length affects both the roHonamt frequency 
and the conductance of the element. A longer strip interacts more with the 
incident field and leads to increased radiation from the slot. The resonant 
fi-equency, on the other hand, decreases linearly with the strip length as shown 
m Fig. 5.20. By adjusting the strip length along with strip inclination, it is 

possible to operate the strip excited slot at any given frequency at a specific 
admittance level. 
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5.3.4 Slot Length and its Effect on Bandwidth 



Slot length 2 (mm) 


Figure 5.21: Variation of resonant frequency and roKonant conductance witii roupt-ct 
to slot length. Q = 80°, 2/, = 9.4 mm, 2wi = 0.8 mm, /// 5.0H mm. . 0 imn, 

2w, = 0.8 mm, e, = 2.2, t = 0.8 mm, 100. .V, • ,V, ■ 9 

The slot length does not affect the admittance behavior a.H mticit a.H the strip 
length. With respect to Fig. 5.21 it can be seen thiit the r(‘.sonant conductance 
decreases with slot length. The resonant frequency also denv-a.sr.s linmriy with 
slot length. The effect of slot length on resonant frequency is much less than 

that of the strip length. The slot length however affects the handwidtli of the 
strip excited slot. 

Fig. 5.22 to 5.23 shows the effect of slot length on the bandwidth of the strip 
excited slot. The bandwidth of the strip excited slot depends on how close the 
element is being operated to the slot’s unloaded resonant frequency. It is to 
be noted that a narrow wall slot in a TO 90 waveguide is much smaller than 
its resonant length. Fig. 5.22 illustrates the conductance behavior of 4 slots 
A, B, C and D made to operate at 8.0, 9.0, 10.0 and 11.5 GHz respectively. 
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Figure; 5.22; Illustrating the bandwidth improvement with respect to the slot 
length, slot: — 0.8 mm; offsets. yi = 5.08 mm, zi = 0 mm; sub- 

strate; f , = 2.2, t = 0.8 mm; basis, modes: Mmax — ^^max — 100, = A^s = 9 
(A) strip: 10.38x0.8 mm, 9°; (B) strip: 8.30x0.8 mm, 20°; (C) strip: 6 91x0.8 mm, 

47°; (D) strip: 7.24x0.6 mm, 78°. 

The strip length and strip angle are adjusted appropriately for this. Note that 
even a 9.5 mm slot at 11.5 GHz is only about 0.36Ao which is shorter than the 
unloaded resonant length (approximately 0.5Ao). The bandwidth achievable 
for any narrow wall vertical slot is limited by the narrow wall dimension of the 
waveguide. This can be improved by opting for a nonstandard waveguide such 
that the nairow wall is wider than the standard waveguide dimension. 

5.3.5 Strip and Slot Width 

The eflfect of strip width is similar to that of strip length. Fig. 5.24 to 5.25 
illustrate 2 cases where the effect of the strip width is presented. It can be 
observed that the resonant frequency decreases and the resonant conductance 
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Figure 5.24: Variation of resonant frequency and resonant conductance with respect 
to strip width. See Fig. 5.4 for other details. 


ictanci 







Strip width 2 (mm) 

Figurtj 5.25: Variation of resonant frequency and resonant conductance with respect 
. to strip width. See Fig. 5.1 for other details. 

increases with strip width. The change in resonant frequency is not linear with 
respect to the strip width but the conductance variation is noticeably linear, 
i Though the effect of strip width can be achieved by adjusting the strip length, 

i it is possible to use strip width as a controlling parameter in a strip excited 

i slot element. Strip widths of the order of 0.4 mm can be fabricated easily in 

the laboratory. We have for example reduced the strip width from 0.8 mm to 
; 0.6 mm while generating case D of Fig. 5.22. 

] Fig, 5.26 shows the effect of slot width on the admittance behavior of the strip 

i excited slot. There is a reduction in resonant frequency with increasing slot 

I * 
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Figure 5 26: Variation of resonant frequency and resonant comiuctanct; with resfifct 
to slot width. See Fig. 5.1 for other details. 


width. This is nonlinear and less prominent than due to the strip width. The 
normalized conductance of the slot is above 2.0 in the graph and it decreases to 
2.0 with slot width. If the slot width were to be increased beyond 0.8 mm, we can 
expect the conductance to fall sufficiently below 2.0, which implies a reduction in 
the radiated power. However, for an array design, it is not necessary to change 
the slot dimensions at all since the required resonant point and admittance can 
be S5mthesized with the strip parameters. 
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5.3.6 Strip Offsets 



Figure 5.27- Normalized shiint conductance vs. frequency at various strip offsets Zi- 
Measurements were done using standard waveguide calibration and no ground plane, 
strip: 8.0 X 0.8 mm, 50°; slot: 8 5x0.8 mm; offsets: yi = 5.08 mm; substrate: = 

2.2, t = 0.8 mm; basis, modes: -1/max = = 50, Ni = 5, = 3 

The cross-over point of the strip with respect to the slot can be altered by offset- 
ting the strip. The strip can be moved vertically, horizontally or both ways. We 
shall call zi as the horizontal offset and yi as the vertical offset. Fig 5.27 to 5.29 
show the effect of horizontal offset on the radiation characteristics. The reso- 
nant frequency increases and resonant conductance decreases with horizontal 
offset. The change in resonant frequency is faster at larger offset values. 

Fig. 5.30 presents the variation of resonant frequency and resonant conduc- 
tance with respect to vertical strip offset. Note that the effect of vertical offset 
is very small. Therefore, it can be assxnned that the effect of a compound offset 
will be governed by the horizontal offset rather than the vertical offset. For 
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Figure 5.28: Normalized shunt susceptance vs. frequency at various strip offsets t;. 
See Fig. 5.27 for details. 



Strip offset Z, (mm) 

Figure 5.29: Variation of resonant frequency and resonant conductance with respect 
to the strip offset Zi. See Fig. 5.27 for details. 


Resonant conductance (norm.) 
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Strip offset (mm) 


Figure 5.30- Variation of resonant frequency and resonant conductance with respect 
to the strip offset yi. offset: Zi = Q mm, See Fig. 5.27 for other details 


a test antenna with a centered strip, we can expect the measured resonant 
frequency to be slightly higher due to alignment errors. 


5.3.7 Substrate Thickness 


Changing the substrate thickness alters the propagation constant and the mode 
patterns of an inhomogeneously filled waveguide. However, for the strip excited 
slot, a change in substrate thickness also amounts to changing the distance 
between the strip and the slot. We can expect an increase in coupling between 
the incident field and the strip with increasing substrate thickness, since the 
strip is now moving into stronger incident field region. However, the strip to 
slot coupling might be decreasing with substrate thickness leading to an overall 
• reduction in slot radiation. It can therefore be expected that the radiated power 

: first increases and then decreases with substrate thickness. Fig. 5.31 confirms 

that this is indeed true. It is to be noted that the normalized conductance first 
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Substrate thickness (mm) 

Figure 5.31. Variation of resonant frequency and resonant conductance with respect 
to substrate thickness, strip: 8. 0x0.8 mm, 50°; slot: 8.5 x 0.8 mm; offsets: iji = 

5.08 mm, Zi = 0mm, substrate: Cr = 2.2; basis, modes: Mnxax = ^\nax = 100, 

.\^ = TV, = 9 

increases up to 2.0 and then increases beyond it. Note that, an increase in 
conductance beyond 2.0 is actually a reduction in radiated power. The resonant 
frequency also increases with substrate thickness and then starts decreasing 
approximately around the same point where the radiation starts dropping. 

5.3.8 Substrate Permittivity 

Fig. 5.32 shows the effect of changing the substrate permittivity. The permittiv- 
ity affects both the resonant frequency as well as the conductance considerably. 
Note that, a change of €r from 1.5 to 3.5 takes the resonant frequency all the 
way from 11.75 GHz to 8.5 GHz. The variation of the conductance with respect 
to Cr is similar to that of substrate thickness. The radiation first increases and 
then decreases with Cr. The relative permittivity also changes the equivalent 




5.4 Array Design Aspects 


165 



Relative permittivity 

Figure 5.32: Variation of resonant frequency and resonant conductance with respect 
to the relative permittivity of the dielectric substrate, strip 9 4x0.8 mm, 80^; 
slot: 8.0 X 0.8 mm, offsets: yi =, 5.08mm, Zi = 0 mm, substrate: t = 0.8 mm; 
modes, basis: Mmax = Nmax — 100, Ni = Ns = 9 

slot length. 


5.4 Array Design Aspects 

The principles of waveguide fed slot array design is exhaustively discussed in 
[48, 17, 16, 18]. A qualitative introduction of array design using the strip excited 
slot will be taken up based on these references. A linear array of equispaced 
slots is considered to make the discussion simple. 

A waveguide fed slot array consists of a wavegmde with a series of equally 
spaced slot elements excited from one end. Each slot radiator is a strip excited 
slot. The array is termed stctTidiTig wavs fed or tTctveliTig wave fed depending on 
the load end termination. A travehng wave array is terminated in a matched 
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load, while a standing wave array is terminated in a short circuit A5/4 away 
from the last slot. In addition, the slots are spaced \g/2 apart in a standing 
wave array, whereas the slot spacing is different in a traveling wave array. 
Irrespective of these differences, the following results are true for any type of 
serially fed array of radiators. 

When an array of slots is fed from one end, each slot radiates a fraction of 
the incident power with the result that a set of slot voltages gets established on 
the slots. In an array, the slot radiators are designed such that the required set 
of slot voltages are established (and maintained) over a designated frequency 
band. When a propagating mode interacts with a slot, the slot radiates as well 
as scatters the incident power to either side inside the waveguide. The scattered 
wave in turn propagates to either end of the array, thereby influencing every 
other slot. This can be regarded as an internal excitation regardless of whether 
the mode is due to the scattering from a slot or due to the array driving source. 
The first few evanescent modes also affect the slot to some extent, and it is 
necessary to take their effect into account if the array is to meet specific design 
goals. Each slot also interacts with every other slot through the exterior of the 
array and this is termed external mutual coupling . 

The situation can be summarized as follows. A slot is excited by three dis- 
tinct mechanisms — the two waves incident from either side internally, and the 
wave incident externally due to the mutual coupling via. other slots. Internal 
interaction is predominantly through the propagating mode, but it is necessary 
to take into account the effect of at least the nearest evanescent mode namely 
LSE20, in addition to the propagating mode. 

If the strip excited slot is represented as a shunt admittance, then the array 
can be represented to a first approximation as an array of shunt admittances in 
an equivalent transmission line. This presumes that we are neglecting both the 
external mutual coupling and higher order internal mutual coupling. Elliott's 
procedure [18] introduces the concept of active admittances to account for mu- 
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tual coupling. Instead of regarding the array as an array of shunt admittances, 
it regards it as an array of active admittances. The active admittance is the 
sum of self admittance along with a term which linearly depends on all the slot 
voltages. This is termed mutual coupling term and it includes the effect of both 
internal and external mutual coupling. What is required of an array, element 
is the ability to achieve any specified shunt admittance by a suitable choice 
of parameters. The shunt admittance is to be adjusted such that the active 
admittance of each slot is made real. Note that, if we were to ignore mutual 
coupling, the resonant conductance variation as a function of the parameters 
of the strip excited slot would have been sufficient for the design. 

In an array design, one starts with a pattern and chooses an appropriate 
number of elements and excitation coefficients to meet this. It is also necessary 
to have an input match and a specified efficiency for the array. The array 
coefficients are identical to the set of slot voltages within a scale factor. The 
array design scheme comes out with the specification of each slot radiator 
iteratively, such that the set of slot voltages are approximately realized. In 
the case of strip excited slot, it is possible to achieve any shunt admittance by 
adjusting the strip length and strip inclination alone. All other parameters can 
be fixed once for all. Therefore, the external mutual coupling effects are not 
changed while iterating to get the appropriate strip length-angle combination. 
The following results are required for designing an array. 

• The shunt admittance variation as a function of strip angle and strip 
length. 

• The voltage induced across the slot, if the slot were excited by LSExq or 
LSE20 mode fi-om either side. 

• The scattered wave amplitudes of both LSEiq and LSE20 modes, had the 
slot been excited by LSEio or LSE20 mode firom either side. 
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• The scattered wave amplitudes of both LSEw and LSE 2 Q modes, if the slot 
were to be excited by another slot having a specified voltage across it at 
an arbitrary location in the array plane. 

The above data set can be generated based on the methods developed in 
Chapter 2. 


5.5 Design Data Generation 

The behavior of the strip excited slot with respect to each of the parameter 
suggests numerous ways of tailoring the excitation level and resonant point. It 
appears sufficient to use only the strip length and inclination in adjusting the 
admittance level of the slot. In order to get better bandwidth, it is necessary 
to choose the slot length close to its unloaded resonant length. In general, this 
cannot be achieved on the narrow wall of a standard waveguide. It appears 
necessary to use waveguides with increased narrow wall dimensions such that 
the LSE 20 cut off is still sufficiently far. 

It is necessary to realize any shunt admittance if the strip excited slot were 
to be used as an array element. However, if the mutual coupling is moderate, 
it is not necessary to realize admittances which are far from resonance. It 
is generally sufficient to generate the admittance values around the resonant 
point. This implies that a slot radiator in an array environment is operated off 
from the resonant point, but not too far from it. 

Fig. 5.33 shows the normalized admittance variation with respect to strip 
length and strip inclination. This data is generated at 11.0 GHz, and the 
normalized conductance range from 0.015 to 0.10 was selected, for the design 
of a 21 element traveling wave array. The strip inclination is varied from 52° 
to 76° in steps of 4° while the strip length is varied from 4.0 mm to 10.0 mm in 
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conductance (normalized) 


Figure 5.33: Admittance design curve: Shunt admittance variation with respect to 
strip angle 6 and strip length 2// at 11.0 GHz. 9 is varied from 52*^ to 76^ in steps 
of 4®, while 2(/ is varied from 4 mm to 10 mm in steps of 0.30 mm. strip: 2iWi = 

0.6 mm; slot: 9 0X1.5 mm; offsets: yi = 5.08 mm, Z; = 0 mm, substrate: Cr = 2.5, 
i = 1.6 mm; modes, basis: Mmax = ^max — 100, Ni = Ns = 9 

steps of 0.3 mm. The following points are worth amplifying: for a given strip 
length, it is possible to adjust the inclination to make it resonant; conversely, 
for a given strip inclination, resonance can be achieved by adjusting the strip 
length. 

Five points are identified from Fig. 5.34, such that they correspond to strip 
inclinations 56°, 60°, 64°, 68.8° and 72°. The strip lengths are chosen such that 
all of them are resonant at 11.0 GHz. This corresponds to the case when all 
the points lie on the zero susceptance line of Fig. 5.33. The corresponding strip 
lengths are foimd to be 5.49 mm, 5.80 mm, 6.14 mm, 6.62 mm and 6.95 mm 
ihom Fig. 5.34. Note that it is necessary to increase the strip length along with 
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Figure 5.34: Variation of resonant conductance and resonant strip length with respect 
to strip angle, from the Admittance Design curve of Fig. 5.33 



Figure 5.35: Comparison of computed and measured conductance values based on 
Fig. 5.33. Measurements were done using TRL standards and no ground plane. 


Resonant strip length (mm) 
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strip inclination to maintain the resonant frequency. 

Fig. 5.35 shows the computed and measured conductance for these 5 cases 
of length-angle combinations. It can be seen that the measured resonant fre- 
quency is about 4.5% higher than the computed value. Interestingly, the reso- 
nant frequencies for the last 3 cases is progressively decreasing. This seems to 
support the conclusion that the error in resonant frequency is directly related to 
the strip inclination. The conductance values are in moderate agreement in all 
the 5 cases. The measured conductances are slightly larger than the computed 
ones. This can be attributed to the losses in the antenna. Note that, the peak 
conductance values in Fig. 5.35 are larger than the resonant conductances from 
Fig. 5.33. This is because the resonant point does not coincide with the peak of 
the conductance value. 

This experiment suggests the relevance of the present analysis in conducting 
a systematic preliminary set of experiments. 


5.6 Summary and conclusions 

A comprehensive discussion of the properties of the strip excited slot was made 
in this chapter. The strip excited slot has numerous parameters which affect 
its behavior. Therefore a proper understanding of the nature of this antenna 
required extensive parameter-wise characterization. The analysis program was 
preferred for generating this data to an extensive series of measurements. This 
required validating the theoretical analysis by comparing it with a select set of 
measurements. The computed results were found to be inaccurate as the strip 
inclination was reduced. A series of experiments were conducted to identify 
the reasons for this behavior. It was found that the approximations used in the 
formulation was not quite true at lower strip inclinations, and this lead to the 


errors. 
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Computed results with respect to 8 parameters of the strip excited slot were 
made. Certain observations are as follows. The strip angle could change the 
conductance of the element over a veiy wide range. By changing the strip 
inclination from 9 to -9 the slot voltage can be flipped. It is found that the 
slot ceases to radiate at both 0° and 90° strip angle. The strip length could 
be adjusted to change the resonant frequency of the element. The variation of 
resonant frequency with respect to strip length was found to be nearly linear. 
The substrate permittivity changes the resonant frequency and conductance 
considerably. The radiated power first increases and then decreases with the 
substrate thickness. Similar effect was found for the slot width. Any relative 
offset between strip and the slot was found to increase the resonant frequency. 
The bandwidth of the element increases as the slot was made to approach its 
unloaded resonant length. 

Array design aspects using the strip excited slot were discussed based on 
Elliott’s work on longitudinal slot arrays. An array element should possess 
the ability to achieve any shunt admittance, if an array is to designed in the 
presence of mutual coupling. For the strip excited slot, a suitable combination 
of strip length and strip angle exist such that any admittance could be real- 
ized. A design curve was generated from which the required strip parameters 
could be read for a given admittance. An experiment was conducted to verify 
the validity of the design curve at a set of 5 points. The theoretical results 
shows moderate agreement in conductance values, however the resonant fre- 
quency was considerably off. This experiment confirms the observation that 
the computed results are inaccurate at lower strip angles. This experiment 
also suggests the significance of theoretically generated results in conducting 
systematic initial trials. 



Chapter 6 


Summary and Conclusions 


A method of exciting the narrow wall vertical non-wrap around slot was pre- 
sented in this thesis. The method makes use of an inclined or an L shaped strip 
for exciting the slot. The strip and the slot can be etched together on either 
sides of a dielectric substrate and can be substituted for the narrow wall. It is 
possible to get good accuracy and square ended slots, since this method makes 
use of printed circuit technology. 


This thesis has presented a complete investigation of the strip excited slot 
as a prospective array element. Both computations and measurements were 
used for an in-depth parameter wise study. The measured results were used 
primarily for illustrating certain effects — ^that of ground plane, losses etc. — ^for 
which computations could not be carried out. 
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6.1 Strip Excited Slot — ^Major Characteristics 

• The strip excited slot is predominantly a shunt discontinuity which goes 
through a radiating series resonance. It also exhibits multiple resonances for 
certain choice of parameters, in which case it has broader bandwidth. 

• An L shaped strip can also be used to couple the slot to an incident LSEw 
mode. The slot was found to get excited (for both inclined strip and L-strip) 
even in the extreme, when the strip does not overlap with the slot. 

• In the absence of the slot, the strip scatters on its own and appears as a series 
L-C circuit in shunt. The measurement, however indicated a significant real 
part which is indicative of the losses in the antenna. This happens, since the 
resonant strip current is much larger than its value at the combined resonant 
frequency when the strip and the slot are present simultaneously. 

• The strip length and the strip angle affect the resonant frequency and admit- 
tance of the element considerably. The resonant conductance increases with 
strip angle. The resonant frequency also goes up with strip angle. Therefore, 
the coupling and resonant point are tied up. However, the resonant frequency 
is found to scale linearly with strip length, so that, a given admittance can be 
realized by a suitable choice of strip length and angle. 

• The length of the slot decides the bandwidth of the element. Bandwidth 
increases (approximately) linearly with the normalized slot length, and for 
practical applications it is necessary to choose the slot length as close to its un- 
loaded resonant length. This necessitates the use of a non-standard waveguide 
with an increased narrow wall dimension. 

• The effect of offsetting the strip horizontally was to increase the resonant 
frequency and to reduce the resonant conductance. The vertical offset of the 
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strip, however did not show any noticeable effects. The radiated power was 
found to ^0 thi ou^h a rnaximum with respect to the substrate thickness. 


6.2 Quadrature Methods for MoM Matrix Entries 

The analysis of the strip excited slot via. integral equation (solved using MoM), 
was exceptionally involved. The inclination of the strip along the .s direction and 
the presence of the factor in the waveguide Green’s function, compli- 

cates the MoM matrix quadratures. This problem was identified as one which 
calls for an elaborate symbolism and an algorithm was developed for imple- 
menting the same. The algorithm was based on two ideas: the recognition of a 
multiple integral as a binary tree, and the representation of trigonometric func- 
tion in exponential form. These two results were used to decouple the source 
and field points such that the quadrature could be expressed algorithmically. 


6.3 Green’s Function as a Distribution 

Green’s functions are distributions, whereas the scattered field derived using 
them are classical functions. The subtle issues involved, in using a Green’s 
function, if one of the source domain quadratures are dropped has been ex- 
plained, It was shown that, in the event of neglecting the x' quadrature, the 
E: field remains a distribution due to the presence of (5(a: - x'). The major re- 
sult of this study is that, for the Ez component to make sense, it is necessary 
to use an LSM represented Dirac Delta. An algorithm for implementing the 
Gram-Schmidt orthonormalization has been developed in this connection. The 
numerical subtraction of Dirac delta is also shown to be unstable in the limit 
of summing large number of modes. 
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6.4 Solution of LSE and LSM Eigenvalues 

The LSE and LSM eigenvalues are used extensively in the computations, since 
the waveguide Green’s function involves a series of LSE and Lb.M normal 
modes. A robust algorithm for solving these values, for an inhomogeneously 
filled waveguide has been presented. The algorithm makes use of a modi- 
fied Newton-Raphson-Bisection root finding. The procedure is based on set of 
simple results for the LSE and LSM curves based on the properties of tan 
functions. These provided good starting values and tight bounds for the root 
finding scheme. 


6.5 Theory and Practice of Method of Moments 

Method of Moment attempts to solve a first kind integral equation with com- 
plete disregard for the essential ill-posedness of the problem. The method does 
not give progressively better solutions as it claims to be, but the sequence of 
solutions initially improves and then goes unstable; a behavior typical of direct 
methods for first kind integral equations. A set of examples were used to illus- 
trate this point. Tikhonov Regularization forms the right basis for explaining 
the instability issues and for developing stable MoM algorithms. 


6.6 Measurement Techniques 

TRL calibration has been used for measuring the S parameters of the strip 
excited slot, since conventional waveguide standards could not conveniently set 
the measurement plane at the slot center. The design and implementation of the 
TRL calibration standards were dealt with. The test antenna was fabricated on 
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a machined wavc‘guide, to which a printed substrate was attached via. soldering. 
The effect of the* giound plane was a shift in the resonant frequency by 0.5%. 
The los.se.s in the sti ucture were found to be substantial if the strip is made 
resonant in the* absence of the slot. 


6.7 Scope for Further Work 

The stuciy carri(>d out in this thesis could be repeated for an L shaped strip, since 
the L-sirii> i.s .simpler to analyze. It is necessary to generate accurate design 
data, and develop an algorithm for implementing an array of strip excited slots. 
The practical is.sues of making a planar array are to be sorted out. Since the 
analysis can be used to compute (to a first order), both the external and internal 
LSE'io mode mutual coupling, this data also should be generated and validated. 

The study on /■.’.■(./. ) Green’s function as a distribution gave a nearly accept- 
able solution to the problem of subtracting a Dirac Delta numerically. Since 
this problem is of interest in inhomogeneously filled waveguides with an axial 
source, the method is to be improved, if the errors due to residual delta is to be 
removed completely. 

A numerically acceptable implementation of MoM should be developed on the 
lines of Tikhonov Regularization. Such an attempt will help the EM community 
consolidate the thumb rules used to get optimum solutions in a numerical 
implementation. It improves the authenticity of the computed results, which 
otherwise is difficult to accept, in the absence of specifying how the author 
reaches at the correct value. 
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